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ABSTRACT 


The  random  fluctuation  of  concentration  levels  in  a 
plume  is  studied  using  a  wind  tunnel  simulation  of  dispersion 
in  the  neutrally  stable  atmospheric  boundary  layer.  A  fast 
response  hot-film  probe  is  developed  and  used  to  measure 
high  frequency  fluctuations  of  helium  tracer  concentration  in 
the  wind  tunnel.  Plume  i ntermi t tency  is  well  described  by  a 
Gaussian  error  function  model  originally  developed  for  turbu¬ 
lence  i ntermi ttency  in  round  free  jets.  The  classic  reflected 
Gaussian  model  for  mean  concentrations  successfully  reproduces 
the  spatial  distribution  of  mean  concentration  in  the  elevated 
tracer  plume.  An  expression  for  the  mean  field  eddy  diffus- 
ivity  is  obtained  from  the  reflected  Gaussian  model,  and 
compares  favorably  with  experimentally-derived  values.  For 
non-uniform  winds,  the  vertical  eddy  diffusivity  has  a  complex 
variation  with  height  above  ground,  becoming  discontinuous 
near  the  source  height.  Measurements  show  that  local  pro- 
duct  ion  and  dissipation  of  concentration  variance  c'  are  not 

in  balance.  In  fact,  as  the  plume  moves  downwind,  production 

2 

becomes  negligible  compared  with  dissipation  of  c  within 

— 2 

the  plume.  A  new  dissipated  Gaussian  model  for  c'  is  derived 
from  the  variance  balance  equation  for  the  special  case  of 
negligible  production  and  fluctuations  which  decay  at  two- 
thirds  the  rate  of  those  in  isotropic  grid  turbulence.  This 


model  provides  a  good  description  of  the  spatial  distribution 
of  concentration  variance  within  the  plume.  Support  is  pro¬ 
vided  for  assumptions  made  in  earlier  work  on  concentration  fluctuations, 
in  that  the  rate  of  variance  dissipation  is  observed  to  be  proportional 
to  local  concentration  variance,  with  the  proportionality 
factor  being  the  reciprocal  of  a  decay  time  scale  which  is 
constant  across  the  plume  and  which  increases  linearly  with 
distance  from  the  source.  A  statistical  intermittent  log¬ 
normal  concentration  model  is  developed  and  successfully 
compared  with  measured  peak-to-mean  ratios  and  complete 
probability  density  functions  of  the  concentration  fluctu¬ 
ations;  the  model  parameters  are  concentration  mean,  variance 
and  i ntermi ttency .  The  highest  peak  concentrations  are 
observed  to  occur  about  two  thirds  of  the  distance  out  to  the 
point  of  maximum  mean  concentration. 
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CHAPTER  I 


THE  SIGNIFICANCE  OF  CONCENTRATION  FLUCTUATIONS 

1 . 0  Introduction 

The  most  serious  consequence  of  high  energy  consumption 
by  industrial  societies  is  the  rapid  depletion  of  fossil  fuels. 
As  the  high  quality  liquid  and  gaseous  hydrocarbon  fuels 
become  exhausted,  a  return  to  the  use  of  coal  becomes  more 
attractive.  However,  one  reason  for  the  reduced  use  of  coal 
was  that  it  is  a  dirty  fuel  which  releases  sulfur  and  parti¬ 
culate  matter  as  byproducts  of  combustion.  Industries  that 
propose  to  burn  coal  for  energy  must  satisfy  the  appropriate 
government  agencies  that  the  surrounding  environment  will  not 
be  adversely  affected  by  atmospheric  emissions  of  sulfur  and 
particulates.  Coal-fired  plants  attract  attention  because 
they  are  increasing  in  number  and  are  among  the  most  visible 
sources  of  atmospheric  emissions;  in  fact,  any  plant  using  a 
chimney  will  also  be  subject  to  the  same  rules. 

The  constraint  most  often  imposed  on  industry  by 
government  is  that  plant  operations  must  not  lead  to  ground 
level  concentrations  of  sulfur  dioxide  in  excess  of  a  specified 
level.  The  most  noticeable  problem  with  such  constraints  is 
that  even  the  most  modern,  clean-operating  plants  are  occasion¬ 
ally  in  violation  of  air  quality  standards,  usually  because 
meteorological  conditions  have  temporarily  prevented  adequate 
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atmospheric  dispersion  of  emissions.  In  fact  it  is  not 
practical  to  specify  a  constant  plant  emission  rate  which 
will  never  under  any  circumstances  lead  to  violations  of  exist¬ 
ing  air  quality  standards. 

Public  disapproval  of  industrial  emissions,  if  focussed 
on  such  unavoidable  violations  of  the  air  quality  standards, 
can  place  both  government  and  industry  in  a  difficult 
position.  For  this  reason  there  is  growing  pressure  for  more 
realistic  constraints  which  require  that  specified  ground 
level  concentrations  never  be  exceeded  more  than  a  smal 1 
fraction  of  the  time.  This  crucial  change  is  practicable 
only  if  there  is  adequate  understanding  of  the  factors  which 
influence  the  frequency  of  occurrence  of  a  specified  concentra¬ 
tion.  This  occurrence  frequency  is  called  the  probability 
density  function  (pdf)  of  concentration  fluctuations,  and  a 
study  of  its  properties  and  the  factors  influencing  its  form 
is  the  purpose  of  this  thesis. 

1 . 1  Concentration  Fluctuation  Probabilities 

It  is  well  known  that  air  pollution  levels  near 
industrial  plants  vary  with  time,  and  a  series  of  pollution 
measurements  will  most  likely  show  short-term  peak  values 
which  are  much  higher  than  the  mean.  These  peaks  may  produce 
the  illegal  pollutant  concentrations  referred  to  previously. 

A  measure  of  this  tendency  is  the  ' pea k-to-mean »  concentration 
ratio,  where  the  peak  concentration  Cp  is  defined  as  that  value 
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which  is  not  exceeded  more  than  a  fraction  P  of  the  time 

(or  equivalently,  that  the  specified  concentration  has  a 

probability  P  of  being  exceeded). 

Revision  of  short-term  air  quality  standards  can  only 

be  done  if  P  is  established  for  every  peak  concentration  of 

interest.  Long  term  statistics  of  full  scale  peak-to-mean 

concentration  ratios  [e.g.  Shruaikar  &  Patel  (1977)]  provide 

useful  information  on  the  relationship  between  P  and  c  . 

P 

Because  concentration  fluctuations  are  a  random  variable,  P 
can  be  computed  for  a  specified  peak  concentration  Cp  from 

rCP  r°° 


P(c  )  = 
P 


1- 


p ( c  )dc  = 


0 


P  (c  ) dc  , 


(1.1) 


where  p ( c )  is  the  probability  density  function  of  the  concen¬ 
tration  fluctuations.  Current  knowledge  about  the  form  of 
the  pdf  has  been  obtained  from  two  approaches: 


t  measurements  of  actual  concentration  time  series 
provide  specific  examples  of  the  pdf,  to  which 
empirical  equations  are  fitted. 

•  physical  arguments  are  used  to  derive  heuristic  models 
for  the  pdf,  which  are  then  compared  with  observed 
values. 


The  empirical  approach  is  guaranteed  to  provide  a  good  descrip¬ 
tion  of  observed  pdf  forms,  and  the  observations  can  provide 
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physical  insight  to  the  processes  causing  concentration 
fluctuations.  Because  the  input  of  physical  insight  is 
minimal,  the  method  is  weak  in  establishing  in  advance  how  the 
pdf  values  should  vary  as  meteorological  conditions  change. 

The  heuristic  approach  is  stronger  in  this  regard  because  the 
physical  processes  causing  concentration  fluctuations  are 
usually  linked  in  some  way  with  local  meteorology.  It  is 
difficult,  however,  to  include  all  the  relevant  physics  in  a 
heuristic  model  of  the  concentration  pdf. 

A  good  example  of  the  pdf  curve-fitting  approach  is 
that  of  Apt  (1976)  who  fitted  a  Weibull  three-parameter  distri¬ 
bution  to  two-week  averages  of  atmospheric  radioactivity  data 
taken  over  periods  of  one  year.  The  Weibull  curve  gives  a 
zero  probability  that  concentration  zeros  occur,  which  means 
that  the  observed  phenomenon  of  pollutant  concentrations 
which  are  intermittently  zero  cannot  be  incorporated  by  this 
model.  This  disadvantage  is  removed  in  the  three-parameter 
lognormal  curve  successfully  fitted  to  full  scale  data  by  Ott 
and  Mage  (1976). 

Examples  of  heuristic  pdf  models  include  the  exponential 
distribution  function  favored  by  Barry  (1974)  and  the  two- 
parameter  lognormal  model  proposed  by  Csanady  (1969,  1973  p. 
228).  Rather  than  determining  the  model  parameters  by  curve¬ 
fitting  to  observed  probability  distributions,  as  required  by 
empirical  models,  the  exponential  model  parameters  are  simply 
the  local  concentration  mean  and  i ntermi ttency ,  while  the 


■ 
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lognormal  model  parameters  are  local  concentration  mean  and 
variance.  If  the  concentration  mean,  variance  and  inter- 
mittency  are  known  for  a  given  meteorological  condition, 
then  the  heuristic  models  allow  a  direct  calculation  of  the 
concentration  pdf.  Typical  distribution  shapes  for  some 
empirical  and  heuristically -derived  pdf  models  are  sketched 
in  Fig.  1.1. 

The  various  probability  models  discussed  above  will 

be  described  in  detail  in  Chapter  VI,  and  a  new  heuristic 

three-parameter  lognormal  model  will  be  developed  and  tested. 

This  intermittent  lognormal  model  is  an  extension  of  the  two- 

parameter  lognormal  pdf  model  developed  by  Csanady  (1969,  1973 

p.  228),  and  incorporates  eddy  dilution  and  concentration  inter- 

mi  ttency  in  a  physically  reasonable  way.  The  model  input 

parameters  are  the  local  values  of  mean  concentration  c, 

2 

concentration  variance  c'  ,  and  concentration  i ntermi ttency  y. 


1 . 2  Mean  Concen trat i on 

The  spatial  distribution  of  mean  concentration  in 
an  elevated  plume  is  usually  considered  to  be  well  known.  The 
most  widely  used  model  for  this  distribution  is  the  reflected 
Gaussian, 


=  S'  ^)exp(-^){exp 


2lTVyaz 


2  a 


-(z-H) 

2a_2 


+  exp 


-  ( z+H ) 

77^ 


(1.2) 
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the  application  of  which  is  considered  in  Chapter  V.  This 
equation,  originally  an  empirical  extention  of  the  Gaussian 
solution  to  the  one-dimensional  Fickian  (constant  eddy 
diffusivity  K)  diffusion  equation 


dc  _ 
dt 


K 


.2- 
3  c 


3x 
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has  recently  been  derived  for  dispersion  in  bounded  flows 
by  Veigele  and  Head  (1978).  Vertical  and  lateral  profiles 
of  mean  concentration  in  an  elevated  plume  have  been  measured 
in  the  present  study,  and  in  Chapter  V  are  fitted  to  Eq. 
(1.2).  It  is  shown  that,  in  order  for  Eq.  (1.2)  to  satisfy 
the  diffusion  equation  under  non-uniform  wind  conditions, 
the  concentration  eddy  diffusivity  must  vary  through  the  plume 
in  a  specific  way.  The  results  verify  that  the  simple  reflec¬ 
ted  Gaussian  model  adequately  describes  the  spatial  distribu¬ 
tion  of  mean  concentration  in  elevated  plumes. 


1  . 3  Concentration  Variance 

Expressions  for  the  spatial  distribution  of  concentra¬ 
tion  variance  c'2  inside  elevated  plumes  are  rare,  due  mainly 
to  lack  of  appropriate  data.  The  meandering  plume  model  of 
Gifford  (1959)  requires  that  the  concentration  variance  be 
zero  at  a  fixed  distance  from  the  moving  plume  centerline, 
while  the  exponential  pdf  model  of  Barry  (1974)  requires  that  the 
standard  deviation  /c'2  he  directly  proportional  to  the  mean  concentration  c 
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within  the  plume.  The  concentration  variance  predicted  by 
these  models  is  due  not  to  atmospheric  turbulence  within  the 
plume,  but  rather  to  large  scale  bodily  motions  (meandering) 
of  the  plume  centerline.  The  models  cannot  account  for  the 
observations  of  concentration  variance  reported  by  Becker, 
Rosensweig  and  Gwozdz  (1966)  for  pipe  flows  in  which  plume 
meander  was  negligible.  The  observations  show  concentration 
fluctuation  intensities  i  =  J c  '  2 /c  in  excess  of  unity  at  the 
plume  centerline,  and  increasing  without  limit  toward  the 
plume  edges. 

To  explain  such  observations,  Csanady  (1967)  has 

derived  from  the  equation  of  continuity  a  conservation  law 

for  the  concentration  variance  in  turbulent  diffusion.  The 

components  of  the  conservation  equation  are  similar  to  those 

found  in  the  energy  balance  of  a  turbulent  mixing  layer  or 

— 2 

boundary  layer.  The  conservation  equation  is,  with  S  =  c'  , 


U 


as 

ax 


2 

as 


2 

as 


=  K(^-f  +  ^)  +  2K 


ay 


az 


(3£)2  +  (l£) 

la y1  v9z; 


s 

T  ,  ’ 


(1.3) 


-J  C. 


advection 


gradient 

flux 


production  dissipation 


where  U  is  the  mean  wind  speed,  K  the  constant  eddy  diffusivity 
of  concentration  mean  and  variance,  and  T, t h e  dissipation 
time  scale.  Csanady  obtains  a  closed  form  solution  for  the 
spatial  distribution  of  concentration  variance  by  assuming 
that  both  c  and  c'2  diffuse  with  constant  diffusivity  K  at 
equal  rates.  This  assumption  leads  to  axially  symmetric 
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profiles  of  c'  ,  and  restricts  the  solution  to  regions  of 
uniform  wind  velocity  far  from  solid  flow  boundaries.  The 
derivation,  solution  and  necessary  assumptions  for  Eq.  (1.3) 
are  discussed  in  more  detail  in  Chapter  V. 

An  approximate  solution  to  the  conservation  equation 
derived  by  Csanady  (1967)  can  be  obtained  by  assuming  that 
the  rate  of  producti on  of  conce ntrati on  fluctuations  is 
exactly  equal  to  the  rate  at  which  the  concentration  fluctua¬ 
tions  are  dissipated.  In  that  case  Eq.  (1.3)  reduces 
immediately  to  a  solution  for  the  concentration  variance. 


2KT 


(^-)2+  (H)2 


3y 


8z 


(1.4) 


where  c'  is  expressed  in  terms  of  the  mean  field.  This 
relation  has  been  used  by  Kewley  (1978)  in  a  study  of  fast 

chemical  reactions  in  a  turbulent  plume.  The  solution  Eq. 

2 

(1.4)  for  c‘  is  unrealistic  because  it  predicts  zero  variance 
on  the  plume  centerline.  Thus,  Kewley  recommends  an  empirical 
correction  to  the  solution,  as  sketched  on  the  next  page. 

The  region  of  most  interest  for  atmospheric  dispersion 
lies  close  to  the  ground,  where  plumes  are  influenced  by 
the  nearby  solid  boundary.  The  solution  to  Eq.  (1.3)  for 
c'2  obtained  by  Csanady  (1967)  is  therefore  an  approximation 
to  the  actual  distribution  of  c*  in  plumes  moving  in  this 
region  of  the  atmospheric  boundary  layer,  because  the  solid 
boundary  presented  by  the  ground  is  not  included  in  the  analysis. 
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Radial  Distance  Off  Plume  Centerline 

However,  without  additional  data  there  is  little  hope  of 
developing  a  more  exact  model.  For  this  reason,  the  present 
study  provides  new  measurements  of  vertical  and  lateral  pro¬ 
files  of  concentration  variance  in  plumes  emitted  near  the 
ground.  Using  these  measurements  as  a  basis.  Chapter  V  shows  that  local 
production  and  dissipation  of  concentration  fluctuations  are  not  in  balance. 

As  the  plume  moves  downwind,  production  becomes  negligible  compared  with 

2 

dissipation  of  c  .  Then,  for  plumes  in  uniform  winds  far  above 
the  ground,  concentration  variance  can  be  calculated  from  the 
simplified  equation 
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U  3x  = 

K(3^  +  ***) 

ay  az 

S 

"  Td  * 

(1.5) 

advection 

gradient  flux 

V '  ■ 

dissipation 

2 

which,  for  Td  =  a  /2K,  predicts  a  Gaussian  lateral  and  vertical 
distribution  for  S  =  c'  ,  identical  to  that  found  for  mean 
concentration  c.  This  observation  leads  in  Chapter  V  to  a 
new  dissipated  Gaussian  model  for  the  spatial  distribution 
of  concentration  variance  c‘  in  plumes  carried  close  to  the 
ground . 

1 . 4  Concentration  I n termi ttency 

Only  scanty  information  is  available  on  the  spatial 
distribution  of  concentration  i n termi ttency  in  elevated 
chimney  plumes.  Barry  (1970)  discusses  field  measurements 
of  plume  i n termi ttency ,  and  relates  these  to  long-term  wind 
direction  fluctuations.  This  approach  is  reasonable  when  the 
plume  centerline  is  seldom  above  the  receptor  location,  but 
does  not  provide  information  on  internal  plume  i n termi ttency 
caused  by  atmospheric  turbulence.  I n termi ttency  y  can  be 
described  by  two  components, 

Y  =  Y0  .  Yj  , 

where  y  is  caused  by  large  scale  meandering  movements  of  the 
0 

plume;  the  contribution  yT  is  due  to  smaller  scale  atmospheric 
turbulence  which  transports  uncontaminated  eddies  at  random 


across  the  receptor  location.  The  turbulent  fluctuations 
producing  y-j-  are  normally  assumed  to  have  periods  of  one 
hour  or  less  in  the  atmosphere. 

The  phenomenon  of  plume  meander  occurs  over  periods 
of  from  1  to  3  hours,  and  is  caused  by  wind  direction  changes 
at  the  high  frequency  end  of  the  fluctuation  range  included 

in  yQ  (the  low  frequency  components  of  yA  are  due  to  wind 

t)  y 

direction  changes  by  weather  systems).  As  mentioned  pre¬ 
viously,  Gifford  (1959)  developed  a  theory  which  predicts 
concentration  fluctuation  probabilities  due  to  meandering 
alone,  with  no  contribution  from  the  randomness  of  atmospheric 
turbulence.  The  turbulent  contribution  to  i n termi ttency  has 
not  previously  been  measured  in  elevated  plumes,  but  some 
guidance  is  available  from  the  work  of  Becker,  Hottel  and 
Williams  (  1  965),  who  examined  scalar  intermittency  in  round  free 
jets.  They  found  that  in  sel f -gene rated  jet  turbulence,  the 
radial  position  of  the  jet  boundaries  fluctuated  in  a  Gaussian 
manner,  so  the  intermittency  factor  could  be  well  described 
by  the  Gaussian  error  function.  New  measurements  of  plume 
intermittency  caused  by  boundary  layer  turbulence  have  been 
made  in  the  present  study,  and  Chapter  IV  shows  that  the 
Gaussian  error  function  accurately  describes  the  spatial 
distribution  of  concentration  intermittency  inside  elevated 
pi umes . 
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1 . 5  Measuring  Concentration  Fluctuations 

In  this  experimental  investigation,  concentration 
fluctuations  were  measured  in  a  laboratory  simulation  of  a 
neutrally  stable  (adiabatic)  atmospheric  boundary  layer. 
Chapter  III  will  discuss  the  rationale  for,  and  the  methods 
used  to  provide  this  simulation.  Briefly,  the  main  reason 
for  this  choice  is  that  the  inability  to  model  all  atmos¬ 
pheric  phenomena  simultaneously  in  the  laboratory  is  far 
outweighed  by  the  ability  to  provide  adjustable,  reproducible 
steady  flow  conditions  in  which  reasonable  averaging  times 
give  statistically  reliable  results. 

The  laboratory  measurement  of  high-frequency  scalar 
fluctuations  in  turbulent  flows  has  attracted  little  attention 
during  the  last  several  decades.  The  most  recent  reviews  of 
current  work  appear  to  be  those  of  Csanady  (1973,  Ch.  7)  and 
Hinze  (1  975  ,  Ch.  3 & 6 ) .  Scalar  properties  of  interest  have 
usually  been  the  fluid  temperature  or  species  concentration. 
The  present  lack  of  fluctuation  data  is  due  in  large  part  to 
the  experimental  difficulties  encountered  in  making  the 
measurements.  The  experimental  problems  generally  include: 
low  s i gnal -to-noi se  ratios  (often  less  than  unity),  detector 
sensitivity  to  extraneous  noise  sources  (environmental  or 
electronic),  insufficient  frequency  response,  inability  to 
obtain  dynamic  calibration,  cumbersome  or  expensive  design, 
and  toxic,  explosive  or  messy  tracer  materials. 
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In  spite  of  these  problems  which  are,  after  all,  not 
unique  to  the  measurement  of  turbulent  scalar  fluctuations, 
several  measuring  systems  have  been  developed  which  have 
given  satisfactory  results;  a  survey  of  the  available 
techniques  is  presented  in  Appendix  A.  These  devices  fall 
into  one  of  three  groups: 

•  light-scattering  probes,  in  which  the  amount  of 
light  scattered  into  a  photo-receptor  is  proportional 
to  the  number  of  tracer  particles  in  a  sample  volume. 

•  heated  element  probes,  in  which  the  voltage  required 
to  maintain  a  hot  wire  or  film  at  a  constant 
temperature  is  proportional  to  the  temperature  and 
composition  of  the  gas  flowing  over  the  element. 

•  conductivity  detectors,  in  which  the  electrical 
conductivity  of  an  aqueous  electrolytic  solution  is 
proportional  to  the  concentration  of  electrolyte 
flowing  over  the  sensor. 

An  assessment  of  the  potential  usefulness  of  various 
fluctuation  detectors  is  presented  in  the  review  by  Fackrell 
(1976).  The  advantages  and  disadvantages  of  the  various  probe 
types  are  compared  in  Table  1.1,  from  which  it  is  apparent 
that  no  single  measurement  technique  works  well  under  all 
circumstances.  The  development  and  testing  of  a  new  hot-film 
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fast  response  concentration  sensor  designed  specifically 
for  the  present  study  is  described  in  Chapter  II, 

1  . 6  Summary 

The  purpose  of  this  study  is  to  provide  new  informa¬ 
tion  on  the  behavior  of  the  probability  distribution  of 
concentration  fluctuations  inside  an  elevated  plume.  The 
intermittent  lognormal  model  developed  for  the  pdf  can  be 
used  only  if  the  spatial  distributions  of  mean  concentration 

_  2 

c,  concentration  variance  c‘  and  concentration  i n termi ttency 
y  are  known  inside  the  plume. 

A  new  fast  response  concentration  sensor  was  developed 
and  used  to  measure  these  parameters  in  a  wind  tunnel  simula¬ 
tion  of  the  atmospheric  boundary  layer.  An  existing  model 
for  the  spatial  distribution  of  mean  concentration  in  plumes 
is  verified,  and  a  model  for  i n termi ttency  in  jets  is  shown 
to  work  well  in  elevated  plumes.  In  addition,  a  new  'dissipated 
Gaussian'  model  is  shown  to  successfully  describe  the  spatial 
distribution  of  concentration  variance  inside  elevated  plumes. 
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TABLE  1 .1 

SCALAR  FLUCTUATION  DETECTORS 


DETECTOR 

TYPE 

ADVANTAGES 

DISADVANTAGES 

LASER  - 
DOPPLER 

LIGHT¬ 

SCATTERING 

PROBE 

[Yang  & 
Meroney 
0975)] 

-  good  for  dispersion  of  particulates 
and  gases. 

-  no  flow  distortion  of  sample  volume- 

-  provides  real-time  data  measurements. 

-  good  frequency  response,  adequate 
spatial  resolution- 

-  simultaneous  fluctuating  velocity 
measurements 

-  point,  sheet  or  volume  sampling 

-  dual  system  measures  spatial 
correlations 

-  robust  probe  design 

-  simple  probe  construction 

-  slow  molecular  diffusion  of  oil 
fog  tracer  (S^  4  x  10^)  limits 
accurate  observation  of  dissipative 
eddies. 

-  large  probe  precludes  use  in  small 
spaces. 

-  difficult  to  measure  overall  fre¬ 
quency  response. 

-  special  device  needed  for  tracer 
generation;  tracer  can  be  messy 
or  toxic. 

-  tracer  may  adhere  to  solid  boundaries. 

-  cannot  take  measurements  flush  to 
solid  boundaries. 

INFRA-RED 

LIGHT¬ 

SCATTERING 

PROBE 

[Motycka  & 
Leutheusser 
(1972)] 

-  good  for  dispersion  of  particulates 
and  gases. 

-  wider  applicability  due  to  small 
probe  size. 

-  provides  real-time  measurements. 

-  adequate  frequency  response. 

-  small  amount  of  flow  distortion  in 
sample  volume. 

-  robust  probe  design. 

-  simple  probe  construction. 

-  slow  molecular  diffusion  of  oil  fog 
tracer  (Sc  -  4  x  101*)  limits 
accurate  observation  of  dissipative 
eddies. 

-  relatively  poor  spatial  resolution. 

-  cannot  measure  velocity  fluctuations. 

-  difficult  to  measure  overall  frequency 
response. 

-  special  device  needed  for  tracer 
generation;  tracer  can  be  messy  or 
toxic. 

-  tracer  may  adhere  to  solid  boundaries. 

-  cannot  measure  flush  to  solid 
boundaries. 

DUAL  - 
ELEMENT  HOT 
WIRE/FILM 
PROBE 

[Way  &  Libby 
(1971)] 

-  good  for  disperion  of  gases 

-  Helium  tracer  (Sc  -  .25)  is  clean 
and  easy  to  use. 

-  wide  applicability  due  to  small 
probe  size. 

-  good  frequency  response  and  spatial 
resolution. 

-  accurate  observation  of  dissipative 
eddies. 

-  simultaneous  fluctuating  velocity 
measurements. 

-  negligible  flow  distortion  in  sample 
volume. 

-  can  measure  close  to  solid  boundaries 

-  requires  very  accurate  calibration, 
careful  use  and  complicated  off-line 
data  processing. 

-  optimum  design  difficult  to  establish. 

-  difficult  to  measure  overall  frequency 
response. 

-  delicate  probe  design  and  construction. 

FLUCTUATING 

TEMPERATURE 

PROBE 

[Fiedler 

(1974)] 

-  good  for  dispersion  of  gases. 

-  temperature  tracer  (Pr  =  .7)  is  clean 
and  easy  to  use. 

-  wide  applicability  due  to  small 
probe  size. 

-  good  frequency  response  and  spatial 
resolution. 

-  accurate  observation  of  dissipative 
eddies. 

-  negligible  flow  distortion  in  sample 
volume. 

-  simple  probe  construction. 

-  provides  real-time  measurements. 

-  can  measure  close  to  solid  boundaries. 

-  may  be  sensitive  to  velocity 
fluctuations. 

-  adequate  signal -to-noise  ratios 
may  be  difficult  to  obtain. 

-  solid  flow  boundaries  act  as 
heat  sinks. 

-  delicate  probe  design. 

-  difficult  to  measure  overall  fre¬ 
quency  response. 

SINGLE  - 
ELECTRODE 
CONDUCTIVITY 
PROBE 

[Gibson  & 
Schwartz 

0963)] 

-  good  for  disperion  of  liquids. 

-  salt  tracer  is  clean  and  easy  to  use. 

-  can  detect  species  concentration  or 
temperature  fluctuations. 

-  good  frequency  response  and  adequate 
spatial  resolution. 

-  small  amount  of  flow  distortion  in 
sample  volume. 

-  robust  probe  design. 

-  provides  real-time  measurements. 

-  slow  molecular  diffusion  of  salt 
tracer  (Sc  =  700)  limits  accurate 
observation  of  dissipative  eddies. 

-  delicate  probe  construction. 

-  difficult  to  measure  overall  fre¬ 
quency  response. 

-  limited  to  dispersion  in  liquids. 

-  needs  redesign  to  measure  flush  to 
solid  boundaries. 
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FIG.  1.1  Typical  Distribution  Shapes  for 

Some  Concentration  PDF  Models 


CHAPTER  II 


THE  HOT-FILM  CONCENTRATION  DETECTOR 
2 . 0  I ntroduct i on 

The  fast-response  sensors  currently  available  for 
detecting  high  frequency  scalar  fluctuations  in  laboratory  or 
industrial-scale  turbulent  flow  were  discussed  in  Chapter  I. 
None  of  the  sensors  are  perfect,  and  selection  of  an  optimum 
probe  type  is  at  best  a  compromise  because  sensor  limitations 
make  some  goals,  such  as  small  probe  size,  real-time  data 
analysis,  high  signal -to-noi se  ratio  and  good  spatial  resolu¬ 
tion,  difficult  to  obtain  from  one  system.  The  purpose  of 
this  chapter  is  to  describe  the  design  and  calibration  of  a 
new  helium  concentration  detector  which  uses  a  single  hot-film 
sensor  as  the  active  element.  This  new  design  provides  a 
useful  alternative  to  other  available  fast-response  scalar 
detectors . 

2 . 1  Probe  Design 

One  of  the  objectives  was  to  develop  a  fast  response 
concentration  detector  that  would  be  small  enough  for  use  in 
a  wind  tunnel  of  cross-section  30  cm  x  30  cm.  An  anemometer- 
driven  hot-film  was  selected  for  the  sensing  element  because 
it  combined  robustness  with  good  frequency  response.  In  order 
to  detect  concentration,  the  user  must  exploit  the  sensitivity 
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of  hot-film  temperature  to  the  thermal  conductance  of  the  gas 
flowing  past  the  element.  The  classic  theory  of  convection 
heat  transfer  in  Holman  (1968,  p.  134)  shows  that  heat  lost 
from  the  film  is  related  to  the  gas  properties  by 

h  ^  y'1/6C  1/3  k  V3  p  1/2  U  V2  . 

P 

Helium  was  chosen  as  the  tracer  material  because  its  properties 
are  much  different  from  those  of  air,  and  it  is  therefore 
easily  detected.  Also,  it  has  the  desireable  attributes  of 
being  clean,  non-toxic,  non-flammable  and  low  in  cost. 

The  details  of  the  probe  design  and  construction  are 
given  in  Appendix  B.  In  brief,  the  fast  response  concentration 
sensor  (FRCS)  is  an  aspirated  probe  in  which  the  velocity  is 
held  constant  by  means  of  a  choked  orifice  downstream  of  the 
hot  film  sensing  element.  This  configuration  is  similar  to  that  used  by 
Blackshear  and  Fingerson  (1962)  and  Kuretsky  (1967)  for  detecting  scalar 
fluctuations  in  a  turbulent  flow.  A  short  fiber  filter  in  the  probe 
inlet  reduces  sensitivity  to  wind  tunnel  turbulence.  The  aspiration 
velocity  is  not  adjustable,  so  the  probe  usually  samples  anisokinetically. 
This  means  that  external  flow  streamlines  must  either  draw  together  or 
move  apart  while  entering  the  probe,  depending  on  whether  the  aspiration 
velocity  is  faster  or  slower  than  the  local  tunnel  speed.  Isokinetic 
sampling  was  not  practical  because  it  requires  that  the  suction  speed  be 
adjusted  to  local  tunnel  speed  at  every  probe  location,  and  changing  the 
suction  speed  causes  the  probe  calibration  to  vary.  It  will  be  shown  in 
Chapter  III  that  corrections  to  the  measured  concentration  time  derivative 
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are  necessary  to  compensate  for  signal  dissipation  caused  by  the  probe 
inlet  filter. 

Following  the  format  of  Table  1.1,  some  attributes 
of  the  new  probe  are  given  in  Table  2.1,  and  an  axial 
probe  cross-section i s  shown  in  Fig.  2.1  along  with  a  close-up 
photograph  of  the  probe. 


TABLE  2.1 

CHARACTERISTICS  OF  NEW  FAST  RESPONSE  CONCENTRATION  DETECTOR 


DETECTOR 

TYPE 

ADVANTAGES 

DISADVANTAGES 

SINGLE 

HOT-FILM 

FLUCTUATING 

CONCENTRATION 

DETECTOR 

-  good  for  dispersion  of  gases. 

-  Helium  tracer  (Sc  -  .25)  is  clean 

and  easy  to  use. 

-  small  probe  size. 

-  good  frequency  response  and  spatial 
resolution. 

-  provides  real-time  measurements. 

-  overall  frequency  response  can  be 
measured. 

-  robust  probe  design. 

-  inexpensive  and  easy  to  construct. 

-  anisokinetic  sampling. 

-  unknown  flow  distortion  of  sample 
volume. 

-  care  required  to  eliminate  sensitivity 
to  velocity  and  pressure  fluctuations. 

-  probe  can  diffuse  or  distort  small 
eddies. 

-  needs  redesign  for  measurements 
close  to  solid  boundaries. 

2 . 2  Calibration  of  Steady-State  Response 

Before  the  concentration  fluctuation  detector  could  be 
used  in  the  laboratory,  its  steady  state  and  transient 
response  to  helium  concentration  had  to  be  determined.  Steady- 
state  calibration  of  the  concentration  detector  required  helium 
air  mixtures  of  accurately  known  composition,  in  a  reproducible 
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inexpensive  form.  The  mixtures  were  obtained  by  using  rotameters 

to  meter  the  flow  of  separate  streams  of  helium  and  air  into  a 

mixing  chamber.  The  helium  was  of  research  grade  and  was 

purchased  in  standard  high  pressure  bottles  from  a  commercial 

gas  supplier.  The  air  was  obtained  from  the  laboratory 

compressed  air  supply,  and  was  filtered  and  dried  prior  to  use. 

Five  glass-tube  rotameters  were  used  to  cover  a  flow  range  of 
3  3 

20  cm  /min.  to  20,000  cm  /min.  The  flow  meters  were 
calibrated  to  an  accuracy  of  ±  2%  of  reading  for  flows  of  air, 
helium  or  Refrigerant  R - 1 2  (molecular  weight  122).  Calibration 
at  low  flow  rates  used  a  soap-bubble  flowmeter,  and  at  higher 
rates  positive  displacement  of  water  in  an  inverted  graduated 
cylinder. 

The  Thermo-Sys terns  1 276-1 0 A  hot-film  sensor  was  driven 
by  a  DISA  55D05  battery-powered  constant  temperature  anemometer, 
chosen  because  of  its  low  noise  characteristics.  In  order  to 
ensure  operation  with  the  proper  hot-film  overheat  ratio  of 
1.5,  the  cold  probe  resistance  was  measured  with  a  high-resolu¬ 
tion  digital  ohmmeter.  The  corresponding  operating  resistance 
was  then  calculated  and  an  external  decade  resistance  box  was 
adjusted  to  properly  balance  the  circuit.  A  battery  voltage 
supply,  voltage  offset  box  and  nulling  volt-meter  were  connected 
to  the  output  of  the  anemometer  so  that  the  zero-helium  out¬ 
put  voltage  could  be  set  at  any  desired  value.  This  flexibility 
was  useful  for  later  signal  processing,  especially  when  measur¬ 
ing  voltages  which  fluctuated  about  zero.  The  Hewlett-Packard 
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5326B  DVM  could  not  correctly  average  a  signal  of  varying 
polarity,  so  a  stable,  adjustable  voltage  offset  capability 
was  necessary.  Low  noise  input  was  obtained  by  using  batteries 
to  provide  the  diode-stabilized  offset  voltage.  Circuit 
diagrams  for  the  voltage  offset  devices  are  presented  in 
Appendix  B. 

A  schematic  of  the  calibration  setup  is  shown  in  Fig. 

2.2.  Early  attempts  at  steady-state  calibration  were  compli¬ 
cated  by  drift  in  the  zero-helium  voltage  from  the  probe. 

The  drift  was  traced  to  small  amplitude,  low  frequency,  random 
temperature  fluctuations  in  the  calibration  gas  mixture  and 
in  the  laboratory.  Drift  of  this  sort  occurs  in  anemometer 
systems  where  there  is  no  temperature  compensating  arm  in  the 
bridge.  The  problem  was  overcome  by  passing  the  mixture  first 
through  copper  coils  immersed  in  water  at  room  temperature, 
and  then  blowing  the  mixture  over  the  probe  inside  the  wind 
tunnel.  The  tunnel  was  operated  at  low  speed  to  provide 
ventilation  past  the  probe  and  to  aid  in  stirring  the  air- 
conditioned  laboratory  air  thoroughly  for  temperature  uniformity. 
In  addition,  the  effects  of  zero  drift  were  corrected  for  by 
switching  the  helium  on  and  off  every  10  seconds  during 
calibration  at  a  fixed  helium  level,  and  subtracting  this  zero 
reading  from  adjacent  helium  readings.  Three  groups  of  ten 
"helium  on"  readings  were  averaged  to  obtain  the  resulting 
probe  output  voltage.  The  continuous  signal  sketched  (next  page) 
on  the  left  was  therefore  "chopped"  as  at  the  right,  so  that  the 
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zero  drift  caused  no  problem.  A  strip  chart  was  used  to 
monitor  the  probe  output;  quantitative  readings  were  obtained 
using  a  Hewlett-Packard  5623B  vo 1 tage-to-f requ ency  counting 
DVM  connected  to  a  line  printer. 


Helium  time  Helium  time 

ON  OFF 


The  calibration  curves  obtained  for  the  probe  showed 
a  linear  relationship  between  output  voltage  and  helium  con¬ 
centrations  of  up  to  1%  (10  ppth)  by  volume.  The  concentra¬ 
tions  observed  during  the  experimental  phase  of  the  work  were 
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usually  much  less  than  this  10  ppth  linear  limit.  A  typical 
calibration  curve  for  helium  is  shown  in  Fig.  2.3.  Six  hot- 
film  sensors  were  calibrated  during  the  experiments;  the 
calibration  constants  ranged  from  1.5  to  2.0  ppth/mv.  One 
sensor  was  also  calibrated  for  Refrigerant  R-12/air  mixtures 
so  the  effect  of  tracer  density  and  molecular  diffusivity 
could  be  studied.  The  R - 1 2  calibration  curve  obtained  is  also 
shown  in  Fig.  2.3.  Note  that  the  probe  output  voltage  decreases 
with  increasing  concentration  for  the  R-l 2  tracer,  and  in¬ 
creases  with  increasing  helium  tracer  concentration.  This 
occurs  because  the  heat  conductance  of  R - 1 2  is  lower  than  that 
of  air,  while  that  of  helium  is  higher.  The  curves  also  show 
that  the  probe  is  about  five  times  more  sensitive  to  helium 
than  to  R - 1 2  concentrations,  and  for  this  reason  helium  was 
used  exclusively  during  the  experiments. 

2 . 3  Calibration  of  Transient  Response 

The  calibration  procedure  described  in  the  previous 
section  showed  that  the  probe  output  voltage  varies  linearly 
with  the  steady-state  tracer  concentration.  The  next  step 
was  to  investigate  the  probe  response  to  a  rapidly  changing 
tracer  concentration.  Ideally,  one  would  wish  for  a  concen¬ 
tration  which  fluctuated  sinusoidally  between  zero  and  some 
known  value.  The  frequency  of  oscillation  would  then  be 
slowly  increased  until  the  probe  was  no  longer  capable  of 
fully  responding  to  the  input.  This  would  be  indicated  by  a 
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decrease  in  the  amplitude  of  the  fluctuating  probe  output 
voltage.  Unfortunately,  this  ideal  calibrating  concentration 
signal  is  impracticable  because  of  the  difficulties  in 
accurately  modulating  the  concentration.  A  more  practical 
calibrating  signal  is  a  step-change  in  concentration. 

Diffusion  acts  to  smear  the  boundary  between  volumes 
containing  different  tracer  concentrations,  so  that  sharp 
gradients  are  difficult  to  maintain.  Therefore,  the  key  to 
calibration  of  transients  is  to  expose  the  probe  to  high 
temporal,  not  spatial,  concentration  gradients.  This  was  done 
in  the  present  study  by  swinging  the  probe  on  a  pendulum 
rapidly  across  the  outlet  of  a  10  ppth  helium/air  jet.  It 
was  felt  that  the  combination  of  sharp  jet  boundary  and  moving 
probe  would  provide  the  required  concentration  gradient.  The 
alternative  of  moving  the  jet  was  rejected  because  of  the 
distorted  and  smeared  jet  boundaries  that  would  result  from 
the  forced  entrainment  of  still  air  into  the  upwind  face  of 
the  jet.  A  schematic  of  the  configuration  is  shown  in  Fig.  2.4. 

A  stand  was  used  to  anchor  a  simple  hinge  device,  into 
which  the  0.5  m  long  probe  support  stem  was  inserted.  This 
resulted  in  a  pendulum-like  configuration,  with  the  sensor 
swinging  through  a  vertical  plane  on  an  arc  of  0.5  m  radius. 

The  pendulum  pivot  friction  was  adjusted  so  that  at  the  bottom 
of  its  arc  the  probe  swung  through  the  jet  at  between  5  and  7 
m/s.  This  speed  was  close  to  the  probe  aspiration  velocity, 
which  was  in  turn  close  to  the  wind  speed  in  the  tunnel  during 


25 


the  experiments.  The  probe  output  voltage  was  monitored  on 
a  storage  oscilloscope  which  was  set  to  trigger  just  before 
the  probe  passed  into  the  helium  jet. 

A  typical  oscilloscope  trace  is  also  presented  in  Fig. 
2.4.  Note  that  the  probe  output  reaches  the  steady  value  abouti 
1  -  2  ms  after  it  first  begins  to  rise.  The  probe  velocity 
across  the  helium  jet  was  about  5  mm/ms.  This  means  the  probe 
moved  5  -  10  mm  while  the  output  voltage  was  rising.  It  is 
possible  that  turbulent  mixing  at  the  jet  boundary  smeared  the 
step  change  across  a  region  this  thick,  but  2  -  4  mm  is  more 
likely  for  the  gradient  region.  Probably  only  the  first  0.5 
ms  of  voltage  rise  was  affected  by  concentration  gradients  at 
the  jet  boundary.  The  3%  overshoot  of  the  final  steady-state 
voltage  is  caused  by  the  electronic  feedback  circuitry  in  the 
DISA  55D05  constant  temperature  anemometer.  This  small  amount 
of  electronically  induced  signal  overshoot  improves  the  probe 
frequency  response. 

2 . 4  Probe-Induced  Signal  Attenuation 

It  is  not  possible  to  separate  the  effect  of  finite 
jet  concentration  gradients  from  the  effect  of  finite  probe 
frequency  response.  Therefore,  a  conservative  estimate  of  the 
probe  frequency  response  is  obtained  by  assuming  that  the 
observed  concentration-time  curve  of  Fig.  2.4  is  due  entirely 
to  non-i ns tantaneou s  response  of  the  probe  to  a  step  change  in 
concentration.  The  frequency  response  is  determined  by  examining 
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the  behavior  of  an  electric  circuit  which  responds  as  in  Fig. 

2.4  to  a  step  change  in  voltage  input.  A  simple  R-C  circuit 
analog,  which  gives  a  frequency  bandwidth  of  800  Hz  (see  Fig. 
2.6),  is  inadequate  because  it  does  not  exhibit  the  initially 
sluggish  response  or  the  subsequent  voltage  overshoot  apparent 
in  Fig.  2.4.  Later  work  by  Wilson  and  Netterville  (1979) 
shows  this  behavior  to  be  more  characteristic  of  a  third-order 
low  pass  filter  circuit  which  includes  a  damped  harmonic 
oscillator  in  addition  to  the  R-C  components.  Such  a  circuit  is 
an  approximate  electric  analog  to  two  processes  oc curing  in 
the  probe-anemometer  system,  shown  schematical  ly  i n  Fig.  2.5: 

•  Inside  the  FRCS,  the  square  concentration  wave-front 
is  diffused  by  the  fiber  filter,  and  acquires  an  S- 
shaped  Gaussian  error  function  distribution  according 
to  the  solution  for  dispersion  in  tubes  formulated  by 

G.  I.  Taylor  (1953).  The  error  function  is  approximated 
by  the  exponential  response  of  an  R-C  circuit  element 
having  a  natural  frequency  =  1/RC  =  1/t^  =  2tt f n  . 

•  The  electronic  feedback  circuitry  of  the  DISA  55D05 

constant  temperature  anemometer  drives  the  probe  voltage 

as  a  damped  harmonic  oscillator,  which  causes  the 

initially  sluggish  probe  response  and  the  subsequent 

signal  overshoot.  The  harmonic  oscillator  has  a 

natural  frequency  uj  =  1/t«  and  operates  with  a 

n  ^ 

critical  damping  ratio  r,. 

It  has  been  shown  by  Wilson  and  Netterville  (1979) 
that  such  a  system  has  the  impulse  response 
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where 


Vout  .  1  _  (t1/t2)2  exP(-t/T1 )  + 

Vin  1“2c(t1/t2)+(t1/t2)2 

+  exp(-ct/x2)sin(/l-c2t/T2  -  $) 
{(1-^2)[1-2c(t1/t2)  +  (t-j/t2)2]}  /2 


(2.1  ) 


(2.2) 


The  time  constants  x-j  and  t 2  and  the  damping  ratio  c  were 
determined  from  a  least-squares  fit  of  Eq.  (2.1)  to  the 
experimental  response  curve  of  Fig.  2.4.  The  optimum  values 
were  =  0.42  ms,  x2  =  0.30  ms  and  c  =  0.49.  The  theoretical 
system  response  Eq.  (2.1)  has  been  plotted  in  Fig.  2.4  using 
these  values;  the  fit  between  theory  and  experiment  is  seen 
to  be  good.  Note  that  for  the  special  case  of  x^  =  x2  =  x 
and  c  =  1/2,  the  impulse  response  (2.1)  would  reduce  to  that 
of  a  third-order  Butterworth  filter. 

The  probe's  attenuation  of  a  steady  unit  sine  wave 
input  of  frequency  f  is  (Wilson  and  Netterville  (1979)) 


=  {[I-O  +  Zc^/t^F2]2  +  [F{2?+(1-F2)t1/t2J]2}  2, 
in  (2.3) 


with  F  =  f/f  ,  where  fn  =  1  / 2xrx 2  =  530  Hz.  The  attenuation 
curve  Eq.  (2.3)  has  been  plotted  in  Fig,  2.6  along  with  a 
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typical  experimental  concentration  fluctuation  spectrum. 

The  fast  response  sensor  attenuates  the  concentration  signal 
by  10%  at  300  Hz,  and  the  signal  is  attenuated  by  a  factor 
of  two  (3  db  down)  at  about  600  Hz.  Concentration  fluctuations 
above  450  Hz  were  not  observed  in  the  wind  tunnel;  the  energy- 
containing  eddies  which  produce  virtually  all  of  the  concen¬ 
tration  variance  occur  in  the  range  0  -  150  Hz,  and  are 
attenuated  by  a  negligible  amount.  When  the  tunnel  speed  was 

systematically  varied  (this  is  discussed  in  Chapter  III), 

2 

the  measured  concentration  variance  c'  proved  to  be  insensitive 

to  the  presence  or  absence  of  high  frequency  signal  components. 

Thus,  it  is  concluded  that  the  probe's  frequency  response  is 

2 

adequate  for  measuring  c'  . 

The  sensor  does  not  appear  to  work  as  well  when  measure¬ 
ments  of  concentration  dissipation  are  desired.  The  natural 
dissipation  of  concentration  fluctuations  occurs  mainly  at  the 
high  frequency  end  of  the  spectrum,  typically  peaking  at  fre- 
quencies  3  to  5  times  higher  than  c1  ,  where  Eq.  (2.3)  predicts 
the  most  probe-induced  signal  attenuation.  An  increasing 
wind  speed,  such  as  occurs  vertically  through  the  wind  tunnel 
boundary  layer,  shifts  the  concentration  spectrum  towards  higher 
frequencies,  thereby  exposing  more  of  the  eddies  producing  time- 
derivatives  to  probe  - induced  attenuation..  Chapter  III  will 
show  that  there  is  considerable  prohe-i nduced  dissipation  of 
the  high  frequency  fluctuations,  more  than  is  indicated  by  the 
theoretical  attenuation  of  Eq.  (2.3).  Specifically,  a  velocity- 
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dependant  signal  attenuation  appeared  in  measurements  of  con- 

—  o 

centration  time  derivative  variance  ,3c\  needed  for 

^  3 1 


estimating  natural  dissipation  rates.  The  question  of  success¬ 
fully  applying  the  fast  response  sensor  to  concentration 
measurements  in  the  wind  tunnel  is  discussed  in  more  detail 
in  the  next  chapter. 
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The  Fast  Response  Concentration  Detector 
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FIG.  2.1b  The  Fast  Response  Concentration  Detector 
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FIG.  2.2  Steady  State  Calibration  Equipment 


FIG.  2.3  Calibration  Curves  for 

Steady  Concentrations 
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FIG.  2.4 


Probe  Response  to  a 
Concentration  Step-Change 
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Concentration  Signal 
Attenuation  Eq.  (2.3) 
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FIG.  2.6  Probe-Induced  Attenuation 

of  Concentration  Spectrum 


CHAPTER  III 


WIND  TUNNEL  MODELLING  AND  CONCENTRATION 
MEASUREMENT  TECHNIQUES 


3 .  0  Introduction 

Atmospheric  dispersion  can  be  studied  by  making  full- 
scale  measurements  in  the  field,  or  alternatively,  the  study 
can  approximate  real  conditions  by  laboratory  simulation  of  the 
atmospheric  boundary  layer.  This  latter  approach  provides 
control  over  flow  parameters  such  as  wind  speed  and  direction, 
and  measurements  in  the  laboratory  can  be  repeated  until 
statistically  reliable  quantities  are  obtained,  without 
concern  for  the  vagaries  of  the  weather.  One  disadvantage  of 
laboratory  simulations  is  that  some  atmospheric  flow  properties, 
such  as  vertical  gradients  of  temperature  and  wind  direction, 
are  difficult  to  simulate.  Also,  the  constraining  roof  and 
side  walls  of  the  tunnel  can  cause  flow  features  which  are 
absent  in  most  atmospheric  flows.  An  example  of  this  is  the 
flow  acceleration  and  turbulence  inhomogeneity  in  the  downwind 
direction  which  was  observed  in  the  present  wind  tunnel  flow. 
However,  for  a  great  many  interesting  dispersion  problems  the 
advantages  of  adjustable,  reproducible  steady  flow  conditions 

far  outweigh  the  disadvantages  of  laboratory  simulation. 

?  4 

Because  stable  c'  values  required  about  10  integral  scales 
of  turbulence  to  be  sampled,  it  is  unlikely  that  the  present 
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wind  tunnel  concentration  measurements  could  have  been 
obtained  in  the  field  with  sufficient  accuracy  and  detail 
to  provide  successful  verification  of  the  theoretical  models 
developed  in  Chapters  IV,  V  and  VI. 

3 . 1  Simulation  of  the  Atmospheric  Boundary  Layer 

Atmospheric  boundary  layers  are  generated  by  the 
natural  interaction  between  rough  ground  and  the  wind.  For 
neutrally  stable  flows  of  sufficiently  high  Reynolds  number, 
all  naturally  grown  turbulent  boundary  layers  are  similar  in 
structure,  so  that  relatively  thin  laboratory  boundary  layers 
reflect  in  reduced  scale  the  structure  of  the  thick  atmos¬ 
pheric  boundary  layer.  The  difficulty  with  laboratory  flows 
is  that  naturally  developed  turbulent  boundary  layers  must 
pass  over  a  very  long  fetch  of  roughness  elements  before  they 
grow  thick  enough  for  model  testing. 

Artificial  stimulation  of  turbulent  mixing  can  be  used 
to  create  thick  turbulent  boundary  layers  within  the  confines 
of  short  wind  tunnel  test  sections.  Photographs  of  the  wind 
tunnel  are  shown  in  Figs.  3.1  and  3.2;  the  components  are 
shown  schematically  in  Fig.  3.3.  This  wind  tunnel  configuration 
generates  an  artificially  thickened  boundary  layer  with  vortex 
generators  (spires)  and  barrier  used  to  stimulate  upwind  mixing. 
The  method  is  similar  to  that  used  by  Counihan  (1969,  1970),  and  reduces 
the  required  boundary  layer  generation  length  from  about  80 
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boundary  layer  heights  to  about  8, 

The  technique  of  stimulating  rapid  boundary  layer 
growth  has  certain  practical  obstacles  that  must  be  overcome. 
One  of  the  most  serious  problems  was  poor  lateral  (crosswind) 
uniformity  of  the  flow  field.  The  non-uni formi ty  was  traced 
to  a  wake  instability  behind  the  spires  causing  the  main  air 
flow  to  become  attached  to  one  or  the  other  of  the  tunnel 
walls.  Uniform  flow  was  obtained  by  slightly  rotating  each 
spire  around  its  vertical  axis,  just  enough  to  prevent  the  wall 
flow  attachment.  To  do  this  effectively  required  a  quick 
measure  of  the  resulting  lateral  wind  profile;  this  was 
obtained  with  a  streamlined  pitot  tube  rake  which  connected 
to  an  inclined  multitube  manometer.  The  effect  of  slight 
changes  to  the  spire  angles  of  attack  could  thus  be  observed 
immediately,  making  the  trial  and  error  adjustment  of  the 
lateral  wind  profile  a  relatively  simple  task.  The  final 
profiles  of  wind  speed  and  turbulence  were  measured  with  DISA 
hot  wires;  they  are  shown  in  Fig.  3.4,  in  which  it  is  seen 
that  the  lateral  uniformity  is  with  in  ±  3%  throughout  the 
central  90%  of  the  test  section  width. 

The  simulated  atmospheric  boundary  layer  employed  for 
this  study  has  also  been  used  by  Wilson  (1977)  to  model 
boundary  layer  flow  over  two-dimensional  ridges.  By  comparing 
mean  wind  and  turbulence  parameters  with  full  scale  data 
reported  by  Couni han  (1975),  Wilson  found  that  the  model  scale 

factor  for  the  present  simulated  atmospheric  boundary 
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adapted  from  Wilson  (1977) 
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layer  is  approximately  4000:1.  Table  3.1  shows  the  various 
turbulence  parameters  measured  in  the  present  study. 

The  model  scale  factor  relating  the  simulated  atmos¬ 
pheric  boundary  layer  to  the  full  scale  is  not  as  important 
for  the  present  experiments  as  it  is  for  problems  dealing 
with  dispersion  around  obstacles.  The  plume  concentration 
measurements  described  in  later  chapters  are  used  to  derive 
models  of  the  spatial  and  statistical  distribution  of  concen¬ 
tration  fluctuations  within  an  elevated  plume.  Most  conclu¬ 
sions  are  drawn  from  examinations  of  the  internal  plume 
structure,  and  the  theoretical  predictions  should  apply  over 
a  wide  range  of  scale  factors.  The  single  source  height  can 
represent  both  tall  industrial  stacks  and  low  sources  emitting 
close  to  the  ground. 

3 . 2  The  Turbulent  Wind  Field 

Hot-wire  anemometers  (DISA  55D30)  were  used  to  measure 
the  mean  and  turbulent  wind  field  in  the  tunnel,  with  a  0.3  mm 
diameter  pitot-static  tube  and  inclined  manometer  providing 
occasional  wind  profiles  for  comparison.  A  pitot-static 
reference  tube  upstream  of  the  spires  was  used  daily  for 
setting  the  wind  tunnel  speed  to  a  constant,  repeatable  value. 
This  convenient  upwind  location  meant  that  the  reference  tube 
never  had  to  be  moved  during  the  experiments.  The  hot  wire 
probes  were  calibrated  for  mean  velocity  in  a  small  wind  tunnel 
using  dynamic  head  pressure  taps  as  a  standard.  Day-to-day 
calibration  adjustments  in  the  experimental  wind  tunnel  were 
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accomplished  by : 

1.  using  a  pitot-static  tube  to  measure  wind  speed 
in  the  test  section  at  a  point  20  cm  above  the 
tunnel  floor,  where  turbulence  was  low. 

2.  replacing  the  pitot-static  tube  with  the  hot  wire 
probe . 

3.  making  fine  adjustments  to  the  probe  overheat 
ratio  so  the  measured  velocity  agreed  with  that 
from  the  pitot-static  tube. 

Mean  velocity  was  obtained  bypassing  the  linearized 
hot  wire  anemometer  output  voltage  to  a  vol tage-to-f requency 
counting  digital  voltmeter  (Hewlett  Packard  5 3 2 6 B )  set  for 
100  sec  averages.  The  turbulent  velocity  components  were 
measured  by  first  filtering  the  linearized  signal  to  remove 
noise  components  above  20,000  Hz,  then  passing  the  signal  to 
an  electronic  signal  processor.  The  processing  included: 

1.  removal  of  low  frequency  mean  and  noise  components 
(<  0.1  Hz)  by  A.C.  capacitive  coupling  of  the 
processor  input. 

2.  signal  amplification,  addition  and  subtraction  to 
obtain  the  fluctuating  velocity  components  u', 

v 1  and  w  1  . 

3.  signal  differentiation  to  obtain  the  velocity 
time  derivative 
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4.  signal  multiplication  to  obtain  the  variance 

components  u'2,  v'2,  w'2,  (j^)2  and  the  Reynolds 
stresses  u'w',  u'v'  and  v'w'. 

The  processed  signals  were  then  averaged  for  100  seconds  on 
the  V-to-F  counting  DVM .  Single  DISA  55P11  hotwire  probes 
(1.25  mm  filament)  gave  good  spatial  resolution  for  detailed 
mean  wind  measurements.  DISA  55P61  X-wire  probes  with  a 
sample  volume  of  about  one  cubic  mm  were  used  for  rms  turbu¬ 
lence  and  Reynolds  stress  measurements. 

The  wind  tunnel  reference  speed  of  10  m/s  was  chosen 
high  enough  to  avoid  Reynolds  number  dependance  in  the  turbu¬ 
lence  and  concentration  measurements.  The  criteria  for 
selecting  this  speed  are  discussed  in  Section  3.5.  Vertical 
profiles  of  the  mean  and  turbulent  velocities  are  shown  in 
Fig.  3.5;  these  values  were  maintained  for  all  dispersion 
measurements.  The  vertical  profiles  are  given  for  x/H  =  12.5, 
midway  along  the  wind  tunnel  test  section.  Measurements  upstream 
or  downstream  of  this  position  exhibited  the  same  profile 
shapes,  but  varied  systematically  in  magnitude.  These  gradients 
in  the  flow  direction  are  shown  in  Fig.  3.6. 

There  are  two  separate  effects  which  produce  the 
observed  alongwind  gradients.  First,  boundary  layers  growing 
on  the  tunnel  sides  and  roof  gradually  reduced  the  effective 
tunnel  cross-sectional  area  so  that  the  mean  velocity  increased 
slightly  in  the  downwind  direction.  Also,  the  surface  roughness 
was  not  able  to  generate  enough  turbulence  to  balance  the 
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decaying  turbulence  from  the  spires  and  barrier.  The  result 
was  the  observed  gradual  decrease  in  turbulence  levels  down 
the  tunnel.  In  later  chapters  it  will  be  seen  that  the 
effect  of  the  axially  nonhomogeneous  flow  field  is  observed 
only  in  the  rate  of  plume  growth,  and  can  be  effectively 
removed  by  correction  techniques. 

The  profiles  of  root  mean  square  (rms)  turbulence 
components  in  Fig.  3.5  have  the  desired  vertical  uniformity 
through  the  plume  dispersion  layer  0^z/H^3,  and  the  boundary 
layer  comparison  in  Table  3.1  shows  very  good  agreement  with 
full  scale  data  for  a  wind  tunnel  scale  factor  of  4000:1. 

This  indicates  that  the  method  of  artificially  stimulating  the 
upstream  turbulent  mixing  process  can  produce  a  realistic 
laboratory  simulation  of  the  full  scale  atmospheric  boundary 
1 ayer . 

3 . 3  Measuring  the  Concentration  Field 

The  fast  response  concentration  sensor  described  in 
Chapter  II  was  used  for  measuring  the  concentration  fluctua¬ 
tions  reported  in  subsequent  chapters.  However,  the  slow  but 
random  zero  drift  from  the  hot  film  anemometer  made  accurate 
mean  concentration  measurements  with  the  fast  response  detector 
very  cumbersome.  The  mean  concentration  field  was  determined 
with  a  slow  response,  high  accuracy  heated  element  conduct¬ 
ivity  detector  developed  by  Wilson  (1977),  shown  schematically 
in  Fig.  3.3.  A  continuous  sample  was  drawn  from  the  tunnel 
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test  section  through  a  1  mm  OD  moveable  probe,  into  the 
detector,  where  it  passed  over  a  thermal  conductivity  bridge 
of  the  same  type  used  in  gas  chromatographs.  A  reference 
sample  drawn  simultaneously  from  upwind  of  the  spires  was 
passed  over  the  opposite  arm  of  the  bridge  to  compensate  for 
fluctuations  in  background  helium  concentration;  zero  drift 
was  monitored  by  systematically  turning  the  helium  tracer 
off  and  on.  The  hel i urn- i nduced  bridge  imbalance,  directly 
proportional  to  helium  concentration,  was  amplified  prior  to 
measurement  on  a  vol tage-to-f requency  counting  digital  volt¬ 
meter  set  for  100  sec  averages,  and  was  continuously  displayed 
on  a  chart  recorder.  The  concentration  readings  from  the 
system  were  reproducible  within  ±  2%  on  a  month-to-month 
basis,  a  performance  which  matched  the  accuracy  of  concentra¬ 
tion  fluctuations  determined  with  the  fast  response  sensor. 

The  signal  from  the  fluctuation  detector  was  processed 
by  analog  electronic  circuits  to  measure  concentration  variance, 
time  derivatives  and  i n termi ttency .  The  basic  sequence  of 
signal  processors  and  monitors  is  shown  schematically  in  Fig. 
3.7.  The  signal  processing  steps  were  to: 

1.  Apply  an  offset  voltage  to  remove  the  pure-air 
component  of  the  hot  wire  signal. 

2.  Use  a  low-pass  filter  to  remove  the  high-frequency 
noise  component  (>  1500  Hz)  from  the  signal. 

3.  Use  capacitive  AC  coupling  (basically  a  high-pass 
filter)  to  remove  the  very  low  frequency  zero  drift 
component  (<  0.1  Hz). 
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Amplify  the  signal  as  necessary  and  send  to  correlator 
and  spectrum  display  for  analysis  of  autocorrelations, 
probability  density  functions  and  spectra. 

Pass  the  signal  through  a  differentiator  circuit  to 
obtain  concentration  time  derivative  . 

O  L 


Square  the  signal  ( c '  or 


8c  1 
8 1 


as  required)  with  a 
,  2 


multiplier  circuit  to  obtain  c‘2  or  (f^— )  . 

a  L 

7.  Use  a  low  pass  filter  to  remove  the  high  frequency 
noise  (>  5000  Hz)  introduced  by  the  signal  processing 
circuits. 

8.  Use  a  vol tage-to-f r equency  counting  digital  voltmeter 

to  take  100  second  averages  of  c'  and  )  . 


3 . 4  Measuring  Concentration  I n termi ttency 


Near  the  fringes  of  the  plume, the  concentration  signal 
spent  an  appreciable  length  of  time  at  the  level  corresponding 
to  pure  air.  A  typical  signal  trace  is  sketched  on  the  next  page,  and 
shows  detectable  levels  of  helium  tracer  to  be  present  only 
about  80%  of  the  time.  This  time  fraction  is  called  the 
concentration  i ntermi ttency  y»  and  varies  from  near  unity  in 
the  plume  core  to  zero  in  the  uncontaminated  air  outside  the 
plume.  It  is  obvious  that  the  concentration  i ntermi ttency 
must  be  known  in  order  to  describe  the  statistical  properties 
of  the  concentration  fluctuations.  The  method  used  to  measure 
signal  i ntermi ttency  was  similar  in  concept  to  that  used  by 
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Becker,  Hottel  and  Williams  (1  965),  and  consisted  of  alter¬ 
nately  amplifying  the  concentration  signal  by  a  factor  of  10 
and  then  clipping  the  signal  at  0.8  volts.  After  three  such 
stages,  a  nonlinear  amplification  step  was  obtained  by  squaring 
the  signal.  This  improved  the  s i gnal - to -noi se  ratio,  because 
the  low  level  noise  components  were  amplified  less  than  the 
larger  signal  components.  After  two  more  amplification¬ 
clipping  stages,  a  low  amplitude  500  kHz  sine  wave  was  added 
to  the  signal  to  provide  regularly  spaced  time  steps  for  a 
frequency  counter.  Zero  drift  was  rough-trimmed  with  a  voltage 
offset  and  nulling  voltmeter  prior  to  the  first  amplification¬ 
clipping  stages,  and  was  fine-trimmed  with  a  voltage  offset 
and  oscilloscope  before  input  to  the  frequency  counter.  To 
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reduce  sensitivity  to  the  amplified  air  noise  over  the  probe 
and  the  electronic  system  noise,  the  sine  wave  amplitude  and 
the  frequency  counter  trigger  level  were  set  so  that  the 
counter  was  triggered  only  about  20%  of  the  time.  This 
resulted  in  a  count  of  about  100  kHz  with  the  helium  off.  The 
counting  procedure  is  shown  schematically  in  the  sketch  below. 


A  low  sine  wave  amplitude  is  used  because,  if  the  noise  was 
submerged  within  a  high  amplitude  sine  wave,  so  that  500  kHz 
was  counted  with  the  helium  off,  any  low  level  tracer  gas 
signals  would  also  trigger  the  counting  process  and  falsely 
indicate  that  no  tracer  gas  was  present.  The  portion  of  time 
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that  tracer  gas  is  present  is  defined  as  the  i nterm i t tency  y, 
calculated  from 

_  i  frequency  with  tracer  on  (-  100  kHz) 

T  frequency  with  tracer  off  (-  100  kHz) 

A  schematic  of  the  i  ntermi ttency  circuit  is  given  in  Fig.  3.8, 
and  results  of  the  i ntermi ttency  measurements  are  discussed 
in  Chapter  IV. 

3 . 5  Reynolds  Number  Effects  on  Mass  Diffusion 

The  Reynolds  number  of  the  flow  in  full-scale  atmospheric 

g 

boundary  layers  is  on  the  order  of  10  ,  whereas  the  Reynolds 

number  of  the  flow  in  the  wind  tunnel  simulation  was  only 

5 

about  10  .  Turbulent  dispersion  in  the  wind  tunnel  must 
therefore  be  independant  of  Reynolds  number  if  the  simulation 
is  to  provide  information  applicable  to  the  full  scale. 

The  wind  tunnel  speed  had  to  be  set  high  enough  to 
avoid  Reynolds  number  effects.  However,  to  keep  the  signal- 
to-noise  ratio  of  the  fast  response  concentration  sensor  as 
high  as  possible,  slow  wind  tunnel  speeds  giving  low  dilution 
and  low  noise  from  impact  pressure  fluctuations  were  desireable. 
Therefore,  the  tunnel  was  operated  at  the  lowest  wind  speed 
which  was  still  high  enough  for  Reynolds  number  i ndependance . 

For  a  tunnel  reference  speed  of  10  m/s,  signal-to-noise  ratios 
for  rms  concentration  were  usually  greater  than  5,  Typical 
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noise  levels  corresponded  to  an  rms  concentration  fluctuation 
of  0.2  ppth.  The  variation  of  mean  concentration  and 
Reynolds  stress  -u 'w‘  with  tunnel  speed  is  shown  in  Fig.  3.9, 
along  with  the  effect  of  changing  the  ratio  of  source  exit 
velocity  to  wind  speed  at  source  height.  The  ratio  c  U T/ Q 
was  independant  of  tunnel  speed  wherever  Reynolds  number 
effects  were  absent.  This  region  corresponded  to  tunnel 
speeds  greater  than  about  10  m/s,  and  this  speed  was  used 
for  all  the  experiments.  The  ratio  of  source  velocity  to 
local  tunnel  speed  was  fixed  at  unity  to  ensure  that  plume 
spread  characteristics  resulted  entirely  from  atmospheric  (as 
opposed  to  sel f-generated )  turbulence. 

3 . 6  Correction  for  Anisokinetic  Sampling 

Anisokinetic  sampling  occurs  when  the  fast  response 
sensor  aspirates  samples  at  a  velocity  different  from  the 
local  tunnel  speed  Uy.  The  wind  speed  increased  vertically 
through  the  tunnel  boundary  layer,  so  there  was  only  one  height 
at  which  the  sampling  occurred  isokinetically  (in  this  case,  at 
z  =  1  cm,  where  Uy  =  Up  =  7  m/s).  Fig.  3.10  illustrates  the 
ingested  streamline  pattern  for  the  case  where  Up  <  Uy. 

The  mass  of  an  eddy  cannot  change  while  moving  through 

the  probe,  so 

meddy  =  (pAut)T  =  (pAut)p  =  (pAut>s  =  constant  , 


(3.1  ) 
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where  p  is  density,  A  the  cross-sectional  area  of  the  ingested 
stream  tube,  U  the  eddy  velocity,  and  t  the  transit  time  for 
an  eddy  to  pass  a  point  in  space.  The  subscript  T  indicates 
conditions  in  the  tunnel,  P  indicates  conditions  just  inside 
the  probe  inlet,  and  S  indicates  conditions  at  the  hot-film 
sensor,  after  the  eddy  has  passed  through  the  filter  at  the 
probe  inlet. 

The  mass  flow  per  unit  time  must  be  the  same  at  any 
point  along  the  ingested  stream  tube,  so  that 


m  =  eddy- —  =  PjAtUt  =  ppA.Up  =  p$A  U 
z  t  r  a  n  s  i  t  1  f  ° 


constant . 
(3.2) 


Combining  Eqs.  (3.1)  and  (3.2)  gives  t<~  =  tp  =  t  y ,  which 
means  the  time  for  an  eddy  to  pass  a  fixed  point  in  space  is 
unchanged  by  an  isokinetic  sampling,  so  that  fluctuation 
frequencies  inside  the  fast  response  sensor  are  the  same  as 
in  the  tunnel. 


The  fluctuation  frequency  fy  within  the  tunnel  must  be  directly 
proportional  to  tunnel  speed;  a  series  of  concentration 
spectrum  observations  at  different  tunnel  speeds  showed  the 
expected  proportional  relationship  between  spectral  peaks  and 
tunnel  speed,  confirming  Eg.  (3.3). 


. 


'  . 
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Fig.  3.9  verified  that  for  sufficiently  high  tunnel 
speeds  the  mean  concentration  c  is  related  to  tunnel  speed  U  y 
and  emission  rate  Q  by 


c  U-j. 

— q-  =  A] (x,y,z)  , 


(3.4) 


where  A-j  is  a  function  only  of  x,  y  and  z,  independant  of  Q 
and  U-p  The  rms  concentration  should  also  vary  with  Q/Uy, 


U 


T  


Q  =  A2  ( x  >y  >  z ) . 


(3.5) 


2 

Measurements  of  c‘  were  made  at  various  values  of  Q  and  Uy; 
the  observed  linear  relationship  between  2  Uy  and  Q  shown 

in  Fig.  3.11  verifies  Eq.  (3.5).  Thus  no  an  isokinetic  sampling 

2 

correction  is  necessary  for  concentration  variance  c'  . 

The  dissipation  <f>c  of  concentration  fluctuations  can 
be  related  to  the  variance  of  the  concentration  time  derivative 
by 


6D 

U-r 


(— ) 
'  ^  +■  / 


3t 


tunnel 


(3.6) 


as  shown  in  Chapter  V  .  Measurements  of  dissipation  rate  were 

required,  so  the  dependance  of  time  derivative  variance  on  tunnel 
speed  (or  Reynolds  number)  was  investigated.  The  measurements 

plotted  in  Fig.  3.11  show  a  relation  between  9c1,  Q  and  Uy 

9 1 
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given  by 


(1^) 


at 


A3(x,y,z) 


measured 


(3.7) 


which  was  quite  startling  considering  that  variations  of  c 
and  Jc ' ^  given  by  (3.4)  and  (3.5)  were  inversely  proportional 
to  Uy  rather  than  its  root.  A  clue  to  the  unravelling  of  this 
mystery  is  found  in  Chapter  V  *  where  it  is  shown  that  concen¬ 
tration  time  derivatives  can  be  related  to  other  properties 
of  the  concentration  field  by  a  diffusion  equation,  so  that 


+ 


a_ 

az 


(Kzff) 


u— 1 

uaxJ  * 


(3.3) 


2 

where  all  the  K's  are  proportional  to  Uy,and  S  =  c'  .  Then, 
using  Eqs.  (3.4)  and  (3.5), 


o,  S[terms  of  0(UyQ2 ) ]  . 


(3.9) 


The  concentration  time  derivative  variance  is  proportional 
to  UyQ^,  so  that 
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1 


UTV2  Q 


N 


(If1) 


3t 


AA(x,y,z)  . 


tunnel 


(3.10) 


Combination  with  Eq.  (3.7)  gives  a  relation  between  actual  and 
observed  rms  time  derivatives. 


■  A5UT 

,r,  2 

( 3  ) 
v  at  ; 

J 

tunnel  ) 

measured 


(3.11) 


The  explanation  for  this  behavior  must  lie  in  the 
sensor  response  characteristics.  According  to  Eq.  (3.11), 
the  measured  rms  time  derivative  shows  an  unexpected  inverse 
dependance  on  tunnel  speed.  However,  from  Eq.  (3.3),  eddy 
timescales  are  not  distorted  by  anisokinetic  sampling.  Thus 
the  observed  behavior  implies  that  high  frequency  c'  fluctu¬ 
ations  (responsible  for  most  of  the  time  derivative  variance) 
lose  amplitude  while  passing  through  the  probe,  probably  due 
to  dissipation  of  small  eddies  within  the  filter  at  the  probe 
inlet.  The  observed  inverse  relationship  between  measured  rms 
time  derivative  and  tunnel  speed  means  that  the  probe's 
dissipative  effect  on  the  high  frequency  portion  of  the  concen¬ 
tration  spectrum  is  directly  proportional  to  the  fluctuation 
frequency,  as  sketched  on  the  next  page.  The  specific  relation  between 
probe-induced  dissipation  and  signal  frequency  (or  equivalently. 
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tunnel  speed)  probably  depends  on  the  shapes  of  the  concentra¬ 
tion  spectrum  and  the  attenuation  curve  for  the  probe;  a 

physical  reason  for  the  apparently  linear  relationship  is  not 
obvious. 


f  signal  frequency 


An  estimate  for  the  magnitude  of  the  proportional i ty 
constant  Ag  appearing  in  Eq.  (3.11)  can  be  obtained  from  an 
expression  derived  by  Hinze  (1975,  p.  290)  relating  the  concen¬ 
tration  and  velocity  fluctuation  microscales. 


(3.12) 


« . 


* 


56 


Here  Sc  =  v/D  is  the  Schmidt  number  (0.25  for  dilute  helium 
into  air),  and  A^  and  Ac  are  the  Taylor  microscales  of  along- 
wind  velocity  and  concentration  res  pec t i vely ,  calculated  from 


(3.13) 


and 


tunnel 


(3.14) 


These  expressions  are  combined  with  Eq .  (3.11)  to  give  a 
probe  'dissipation  velocity'  u^. 


-V2 


u 


A 


U-j-S  c 


N 


u 


measured 


(3.15) 


which  is  used  to  correct  the  measured  time  derivative  variance 
according  to 


( 


2 

tunnel 


(icLl) 

v3t  ' 


measured 


(3.16) 


Measurements  of  velocity  and  concentration  time  derivative 
variance  were  made  in  a  vertical  plane  through  the  plume 
centerline  at  x/H  =  12.5;  the  resulting  estimates  for  U.  are 
presented  in  Table  3.2. 


. 
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TABLE  3.2 

PROBE  DISSIPATION  VELOCITY  FOR  CONCENTRATION 
TIME  DERIVATIVE  VARIANCE 


Within  experimental  accuracy  the  fast  response  probe's  dissipation 
velocity  is  Ud  =  7  m/s  -  remarkably  close  to  the  probe's 
aspiration  velocity  Up  =  6.94  m/s,  although  this  may  be 
coincidence.  Thus,  the  concentration  time  derivatives 
measured  in  this  study  were  corrected  for  dissipation  in  the 
anisokinetic  fast  response  sensor  according  to  Eq.  (3.16) 
with  =  7  m/s. 

It  is  interesting  to  substitute  (3.16)  back  into  the 
relation  (3.6)  between  dissipation  and  time  derivative.  Then 


(ic\ 
Kdt  ' 


measured 


» 


which  means  that  the  measured  profiles  of  concentration  time 


58 


derivative  variance  are  directly  proportional  to  the  actual 
dissipation  rate  of  concentration  fluctuations.  As  a  result, 
the  anisokinetic  fast  response  sensor  can  be  used  to  measure 
local  concentration  dissipation  rates  without  the  need  for 
knowledge  of  the  local  wind  field. 


. 
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FIG.  3.1  The  Wind  Tunnel  Laboratory 
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FIG.  3.2  The  Tunnel  Test  Section 
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FIG.  3.3  Wind  Tunnel  System  for  Plume  Simulation 
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U  Wind  Speed  (m/s) 


FIG.  3.5  Vertical  Profiles  of  Mean  Wind  and 

Turbulence  on  Plume  Centerline 
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FIG.  3.7  Equipment  for  Measuring 

Concentration  Fluctuations 


FIG.  3.8  Equipment  for  Intermittency 

Measurements 
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FIG.  3.10  Anisokinetic  Sampling  by  Aspirated  Probe 
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3  Q  Emission  Rate  (cm3/min) 


FIG.  3.11  Dependance  of  Concentration  and 

Time  Derivative  Variances  on 

Emission  Rate  and  Wind  Speed 


CHAPTER  IV 


CONCENTRATION  INTERMITTENCY  AND  RATES  OF  PLUME  GROWTH 
4. 0  I ntroduc ti on 

One  object  of  this  study  is  to  construct  a  simple 
but  flexible  statistical  model  for  concentration  probabili¬ 
ties  in  elevated  plumes.  The  statistical  model  developed 
in  Chapter  VI  achieves  simplicity  through  pdf  shape  parameters 
written  in  terms  of  the  local  values  for  mean  concentration 

_  2 

c,  concentration  variance  c'  and  concentration  i n termi ttency 
y.  Model  flexibility  results  from  the  implicit  inclusion  of 
atmospheric  physical  processes  through  the  spatial  distribu- 
tions  of  these  parameters;  predicted  values  for  c,  c‘  and 
y  can  be  used  in  the  statistical  model  to  forecast  concentration 
p  r  o  b  a  b  i  1  i  t  i  e  s  . 

This  chapter  deals  with  the  downwind  spread  of  the 
mean  concentration  field,  and  relates  this  rate  of  plume 
growth  to  the  velocity  autocorrelation  function.  Then,  the 
spatial  distribution  of  concentration  i n termi ttency  in  elevated 
plumes  is  shown  to  follow  a  Gaussian  error  curve  similar  to 
that  for  turbulence  i n termi ttency  in  free  jets.  The  width 
of  plume  in termi ttency  profiles  appears  to  be  proportional 
to  the  width  of  the  mean  concentration  field. 
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4.1  Measurements  of  Plume  Spread 

Simple,  reliable  models  for  predicting  c'2  and  y  in 
chimney  plumes  have  not  yet  appeared  in  the  literature.  How¬ 
ever,  the  spatial  distribution  of  mean  concentration  c  in 
plumes  has  been  widely  studied,  and  a  number  of  predictive 
techniques  for  c  are  available.  The  simple  reflected  Gaussian 
model  introduced  in  Chapter  I  is  shown  in  Chapter  V  to 
provide  a  good  fit  to  mean  concentration  profiles  obtained 
in  the  present  study.  The  reflected  Gaussian  model  Eq.  (1.1) 
has  the  form 

C  =  (  2  tt  U  ^  a  a  )  e*P  )  {exp[^4P-  ]  +  exp[' ( ]}  , 

Hy  z  2a  1  Zo  Zo  ^ 

with  Q  the  source  emission  rate  and  UH  the  wind  speed  at 
source  height  H;  y  and  z  are  the  lateral  distance  from  plume 
centerline  and  height  above  ground,  respectively.  The  plume 
geometry  parameters  a  and  a  ,  functions  of  downwind  distance 
x,  are  the  lateral  and  vertical  standard  deviations  of  the 
plume  concentration  profiles.  They  describe  the  rate  at  which 
the  plume  grows  as  it  moves  downwind. 

The  literature  contains  literally  hundreds  of  reports 
on  full  scale  measurements  of  and  az  in  elevated  plumes. 

One  example  is  that  of  Singer  and  Smith  (1966),  which  has  been 
recommended  by  the  ASME  (1973).  Values  of  a and  a z  observed 
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in  wind  tunnel  simulations  of  plume  dispersion  must  agree 

with  these  full  scale  data  before  a  successful  simulation  can 

be  claimed.  Profiles  of  mean  concentration  measured  for  the 

present  study  were  compared  with  values  calculated  from  the 

reflected  Gaussian  Eq.  (1.1).  A  least-squares  error  technique 

was  used  to  obtain  values  for  a  and  a  z .  The  location  of 

plume  cross-sections  measured  in  the  wind  tunnel  is  shown  in 

Fig.  4.1;  the  best-fit  values  of  a  and  a  are  given  in 

y  z 

dimensionless  form  in  Table  4.1. 

TABLE  4.1 

DIMENSIONLESS  STANDARD  DEVIATION  FOR 
PROFILES  OF  MEAN  CONCENTRATION 


The  measurements  used  to  obtain  the  velocity  integral  length 

scale  A.  =  $UT_  will  be  discussed  in  Section  4.2. 

Lu  Eu 


< 
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Taylor's  (1922)  statistical  theory  of  turbulent 
diffusion  predicts  that  sigma  increases  linearly  with  x 
close  to  the  source,  and  approaches  a  square-root  dependance 
on  x  at  large  distances  from  the  source.  Large  distances  are 
defined  as  x  >>  A^u, where  =  UT^  is  the  Lagrangian  integral 
length  scale  of  velocity  fluctuations,  discussed  in  the 
next  section.  In  either  case,  sigma  should  depend  linearly 
on  the  (constant)  turbulence  intensities  i  =  / v ' ^ / U  and 

/  2 

i  =  /  w'  /U  respectively. 


The  generalized  x-dependant  form  a  =  ax13  includes  both 
theoretical  extremes  of  linear  and  square-root  variation  and, 
when  fitted  to  the  data  using  least-squares  error  techniques, 
gives  (with  a  and  x  in  cm). 


1£;o  .806  f  2  qqq\ 

a  =  .159  x  (r  =.999) 


(4.2) 


and 

a  =  .165  x'724  (r2  =  .979)  .  (4.3) 


The  observed  axial  variation  lies  between  the  two  theoretical 
extremes,  with  more  tendency  to  linearity  than  to  a  square- 
root  dependance.  However,  Chapter  III  described  how  the  mean 
flow  along  the  tunnel  centerline  accelerated  due  to  boundary 
layer  growth  on  the  tunnel  walls,  and  turbulence  decayed 
because  the  roughness  elements  could  not  maintain  the  high 
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levels  of  turbulence  generated  by  the  upwind  spires  and 
barrier.  As  a  result  the  turbulence  intensities  decreased 
with  x  (in  cm)  at  the  rate 

1  v  =  .  145  -  .003  x'63 

and  (4.4) 

i  =  .1  22  -  .003  x-62  . 
w 

This  non-homogeneous  turbulence  field  means  that  statistical 
properties  of  homogeneous  turbulence  are  just  approximations 
to  the  wind  tunnel  flow  properties.  For  homogeneous  turbu¬ 
lence,  the  rate  of  lateral  plume  spread  is  proportional  to 
the  (constant)  turbulence  intensity  according  to 


o=2 

y 


'2  ( t-T  )  R  (x )  dx  =  2v  1  2 


0 


(t-T)Rv(x)dT  , 


0 


(4.5) 


where  the  autocorrelation  coefficient  R  ( t )  is  defined  as 


R  (x)  e 

v  — ? 

..  i 


(4.6) 


An  expression  similar  to  Eg.  (4.5)  can  be  written  for 


o  . 
z 


H.:  •  *f  -  t,  e,'j  C*  ' ’  1  •»'$  ^0  C  3 
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A  first  approximation  to  the  value  of  a  in  non- 

T 

homogeneous  turbulence  is  calculated  by  replacing  v'  in 
Eq.  (4.6)  by  its  spatial  average  value  over  the  range 
0  5  x  ^  Ut.  This  approximation  is  obtained  by  noting  that 
the  fluctuating  velocity  is  strongly  sel f -correl a  ted  ( R  v  >  0 ) 
for  about  one  integral  length  scale,  a  short  distance  within 
which  the  turbulence  levels  remain  almost  constant  (i.e. 
locally  homogeneous  flow).  This  permits  the  velocity  variance 
to  be  defined  as  the  spatial  average  value,  which  can  be 
removed  from  the  integrand  of  Eq.  (4.5).  These  average 
turbulence  values  are,  from  Eq.  (4.4), 


i v ( x )  =  .1  45  -  .002  x 


.63 


and 


(4.7) 


i  (x)  =  .122  -  .002  x-62. 
wv  ' 


Then  combining  Eqs.  (4.2),  (4.3)  and  (4.7)  in  the  form 


a  =  const  •  i  x  ^  gives 


a  =  .95  i  x 

y  v 


.88 


,  .  0  - —  .  8 1 

a  =  1  . 1  3  i  x 
w  w 


and 


(4.8) 


1 3  'I 
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Using  the  model  scale  factor  of  4000:1,  these  expressions 
for  0^  and  a z  in  cm  can  be  transformed  to  meters  in  the  full 
scale  by 


and 


1  .46 


(4.9) 


a 


z 


2.25 


The  average  value  for  i  through  the  test  section  is 
( i  v ) =  0.129.  Similarly  the  average  value  for  is 

(i  )  =  0.107.  Using  these  average  values  in  the  above 
equations  gives  the  simulated  full  scale  power  laws  (in  meters) 

n  ,  n  .88 
o  =  0 . 1  9  x 

y 

and  (4.10) 

ct  =  0.24  x-81  , 
z  * 


which  are  compared  in  Fig.  4.2  with  Singer  and  Smith's  (1966) 
full  scale  one-hour  average  values  for  a and  o z  measured  at 


Brookhaven  National  Laboratory. 

Values  of  a  and  a  for  the  simulated  plume  were 

y 

derived  from  concentration  field  measurements  obtained  using 
100  sec  averages.  A  direct  application  of  the  4000:1  model 
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scale  factor  would  indicate  that  these  plume  spreading  rates 

correspond  to  approximately  4-day  averages  of  sigma  in  the 

full  scale,  rather  than  the  1-hour  averages  shown  in  Fig.  4.2. 

However,  measurements  showed  the  wind  tunnel  turbulence  spectra  to 

peak  around  30  Hz,  with  little  energy  below  about  1  Hz; 

this  lower  frequency  corresponds  to  approximately  1-hour 

fluctuations  in  the  full  scale.  Because  the  wind  tunnel  does 

not  simulate  the  slow  velocity  fluctuations  which  cause  full 

scale  sigma  values  to  continuously  increase  as  averaging  time 

is  extended  beyond  one  hour,  the  wind  tunnel  sigma  values  were 

essentially  independant  of  averaging  time  after  just  a  few 

seconds;  the  long  100  sec  averages  were  necessary  to  smooth 

out  random  fluctuations  in  the  measured  concentration  variance 
2 

c '  .  The  conclusion  to  be  drawn  from  Fig.  4.2  is  that,  out  to 
the  farthest  modelled  distance  of  2.5  km,  and  for  averaging 
times  up  to  one  hour  full  scale,  the  simulated  atmospheric 
boundary  layer  disperses  material  at  the  same  rate  as  a  neutral 
to  slightly  unstable  full  scale  boundary  layer. 

4 . 2  Velocity  and  Concentration  Autocorrel a t i ons 

The  rate  of  plume  growth  is  related  to  the  turbulent 
velocity  fluctuations  through  equations  (4.5)  and  (4.6).  It 
should  be  possible  to  predetermine  plume  spread  rates  solely 
from  a  time  series  of  the  lateral  and  vertical  wind  components. 
The  plume  spreading  rate  has  asymptotic  limits  which  depend 
on  the  length  of  time  diffusion  has  been  taking  place.  For 
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short  times  (t+0,  R(t)+1)  equation  (4.5)  reduces  to  the  linear 
(a^t)  form 


a i  =  u i  t  (4.11  ) 

and  for  long  diffusion  times  (  R(t)+0)  equation  (4.5) 

can  be  written  as 

ay2  =  TLv  (t-TLv)  ,  (4.12) 

where  ?Lv  is  the  Lagrangian  integral  time  scale  obtained  from 
integrating  the  autocorrelation  function  for  long  diffusion 
times, 


T 


Lv 


3 


°0 

0 


REv(x)dT  . 


(4.13) 


Following  Hay  and  Pasquill  (1961),  the  empirical  constant  3 
is  a  rescaling  factor  for  the  time  axis  which  approximates 
the  relation  between  Lagrangian  and  Eulerian  time  series  of 
turbulent  variables.  Full  scale  observed  values  of  3  vary 
with  turbulence  intensity  i  [Pasquill  (1975)],  going  from 
about  1  to  10  as  the  atmosphere  changes  from  strongly  unstable 
to  strongly  stable  (Nappo,  1979).  The  values  group  around 
3-4  for  neutrally  stable  flows  such  as  occur  in  the  wind  tunnel. 
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The  Lagrangian  integral  time  scale  of  velocity 
fluctuations  is  an  important  parameter  because  it  character¬ 
izes  the  persistence  of  the  largest  eddies  contained  in  the 
plume,  and  can  be  used  to  calculate  the  final  asymptotic 
rates  of  plume  spread.  A  similar  Lagrangian  integral  time 
scale  can  be  defined  for  the  scalar  concentration  fluctuations. 


TLc  “  6TEc  =  B 


.00 


0 


REc(x)di  , 


(4.14) 


where 
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Ec 


c  1  ( t  )c 1 ( t+x) 

T2 

c 


(4.15) 


Hinze  (1975,  p.  290)  shows  that  the  ratio  of  the  integral 
scale  of  vector  velocity  fluctuations  to  that  for  scalar 
concentration  fluctuations  should  be  proportional  to  the  square 
root  of  the  Schmidt  number  when  diffusion  times  are  large 
( t-*») , 


^  =  b/sT  .  (4.16) 

The  constant  of  proportionality  B  is  approximately  unity. 

The  autocorrelations  of  concentration  and  alongwind 
velocity  were  measured  in  the  wind  tunnel  on  the  plume  center- 
line  at  x/H  =  12.5,  and  the  integrations  (4.13)  (with  u  in 
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place  of  v)  and  (4.14)  were  performed  graphically;  several 
observed  autocorrelations  are  reproduced  in  Fig.  4.3.  The 
proportionality  constant  was  (with  Sc  =  0.25  for  dilute 
helium  into  air) 


B 


2  (z/H  =  0.2) 


and 


(4.17) 


B  -  5  (z/H  =  1 )  . 


These  values  for  B  are  quite  a  bit  larger  than  unity,  and 
are  believed  to  be  due  to  younger  'ages'  for  the  observed 
concentration  fluctuations  as  compared  to  velocity  fluctu¬ 
ations.  At  the  measurement  location  x/H  =  12.5,  the  concen¬ 
tration  fluctuations  had  been  diffusing  for  only  t  =  x  / U 
-  46  ms,  whereas  the  velocity  fluctuations  generated  at  the 
spires  and  barrier  had  been  diffusing  for  t  =  Xu/U  ~  140  ms, 
and  were  effectively  three  times  older.  The  relation  (4.16), 
which  links  the  integral  scales  of  velocity  and  concentration 
computed  from  Eqs.  (4.13)  and  (4.14),  is  strictly  valid  only 
for  large  diffusion  times  t  =  t  =  t  ->•  «>.  The  observed  large 
values  of  B,  possibly  resulting  from  low  values  of  T^c  because 
tc  was  only  one-third  of  ty ,  are  not  surprizing  in  light  of 
the  difference  in  diffusion  times  for  concentration  and  velocity 


fluctuations. 
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4.3  Plume  Intermittent 

Concentration  measurements  near  the  boundary  of  a 
turbulent  plume  include  time  periods  during  which  no  concen¬ 
trations  are  observed.  If  all  such  segments  of  the  time 
series  are  removed,  the  remaining  portions  account  for  a 
fraction  "y"  of  the  original  concentration  record.  The 
parameter  y,  the  fraction  of  time  during  which  the  plume 
concentrations  are  non-zero,  is  called  the  concentration 
i ntermi ttency  factor,  and  is  similar  in  concept  to  that  used 
for  turbulence. 

The  i ntermi ttency  factor  approaches  zero  at  the  edges 
of  a  plume,  where  uncontaminated  eddies  are  entrained  from 
the  surrounding  atmosphere  and  the  contaminated  eddies  spread 
outward.  On  the  plume  centerline  y  often,  but  not  always, 
approaches  unity.  In  this  chapter  and  Chapter  VI,  frequent 
reference  will  be  made  to  the  i ntermi ttency  of  concentration 
fluctuations,  as  simply  ' i ntermi ttency ' .  References  to  inter- 
mittency  of  turbulent  velocity  will  be  plainly  stated  as  such. 

Hinze  (1975,  Ch.  6)  reviews  previous  studies  of  the 
eddy  exchange  processes  which  cause  turbulence  i ntermi ttency. 
Measurements  of  the  instantaneous  boundary  of  turbulent  flows 
show  a  highly  contorted  surface  of  bulges  and  valleys  which 
appear  and  disappear  at  random.  Hinze  reports  that  Corrsin 
and  Kistler  (1954)  have  found  from  work  on  a  round  free  jet 
and  a  plane  wake  flow  that  the  statistical  variations  in  the 
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height  of  the  bulges  follow  a  Gaussian  distribution.  Then  the 
integrated  bulge  position  gives  an  expression  for  the  turbu¬ 
lence  intermittency  factor  in  terms  of  the  Gaussian  error 
function, 

Y  =  j  erfc 

where  r is  the  radial  position  at  which  the  turbulence  in¬ 
termittency  has  dropped  to  y  =  |,  and  is  the  standard 
deviation  of  the  turbulent  boundary  bulges  around  this  position. 
The  measurements  of  Corrsin  and  Kistler  indicate  that  both 
r\.  and  aD  increase  linearly  with  downwind  distance  x,  although 

'  2  D 

at  different  rates  and  apparently  from  different  victual  origins. 
More  recent  work  by  Becker,  Hottel  and  Williams  (1965,  1966) 
on  free  and  ducted  round  jets  shows  scalar  intermittency 
distributions  are  also  well  reproduced  by  Eq.  (4.18).  However, 
they  found  that  the  ratio  r^/c^  varied  so  slowly  with  x  that 
al 1  the  radial  i n termi ttency  profiles  could  be  represented  by 
a  single  average  value  for  r v2 / a B •  The  elevated  plume  inter¬ 
mittency  profiles  measured  as  part  of  the  present  study  also 
show  this  convenient  behavior. 

The  electronic  circuit  used  for  the  measurement  of 
vertical,  lateral  and  axial  profiles  of  intermittency  in  an 
elevated  plume  is  described  in  Chapter  III.  The  measured  inter¬ 
mittency  profiles  are  shown  in  Figs.  4.4  to  4,6  along  with 
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ri 
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fitted  Gaussian  error  function  curves  of  the  form  given  in 
Eq.  (4.18).  The  spatial  distribution  of  i n termi t t ency  in  the 
elevated  plume  was  well  reproduced  by  the  product  of  three 
complementary  error  functions. 


Y  =  ^erfc 


T2-  (~-l) 

erf  c 

y‘4  (*-  -i) 
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with  the  empirical  values 
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-  2/2~  , 


xy 


(4.19) 


yi/2 

- —  -  5.3/2“  ,  (4.20) 

yy 

Zi-H 

— —  -  4.6/2". 
azy 

Within  experimental  accuracy  these  ratios  were  independant 
of  the  downwind  distance  x,  a  behavior  observed  by  Becker, 
Hottel  and  Williams  (1965)  for  i ntermi ttency  within  an  axi- 
symmetric  free  jet.  Their  i nterm i ttency  profiles  were  best 
represented  by  the  shape  factor  ri ~  4/2"  . 
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The  variation  of  yv  and  z1 -H  with  x  is  tabulated 
in  Figs.  4.4  and  4.5.  It  is  reasonable  to  expect  that  the 
magnitudes  of  y  y2  an(i  zy2  are  related  to  the  atmospheric 
turbulence  intensities  iv  and  i^  as  well  as  downwind  distance. 
This  dependance  is  conveniently  included  by  relating  variations 
in  y  i/2  and  to  the  values  of  a  and  .  Two  possible 
relations  are  the  power  law  fits 


yi/2  =  3.2  c  ‘ 7  , 

Zi  -H  =  3.6  a  *  ^ 
'2  z 


and  the  linear  fits 


1.7  (a  +3d  )  , 
y  s’ 


zly  -H  =  1.9  (a  +3d  ) , 

/  2  Z  S 


(4.21  ) 


(4.22) 


where  d  is  the  source  diameter.  There  is  no  theoretical 
s 

guidance  for  choosing  between  the  linear  or  power  law  forms. 
Within  experimental  accuracy,  both  fit  the  date  equally  well. 
Both  contain  three  different  constants,  with  the  leading 
constants  in  the  ratio 


3.6/3. 2=1.9/1.7=1.12. 
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It  is  not  known  what  external  factors  affect  the  values  of 
the  constants,  or  the  ratios  between  them.  The  studies  reported 
by  Hinze  (1975,  Ch.  6)  are  of  little  help  in  this  regard 
because  they  refer  almost  entirely  to  velocity  intermi ttency  in  the 
sel f-genera ted  turbulence  of  wakes  and  jets  in  quiescent  or 
laminar  flows.  The  most  significant  difference  between  such 
flows  and  the  present  case  is  the  effect  of  boundary  layer 
turbulence  which  governs  the  growth  rates  a  and  a  (and  hence 

y i.  and  z i,  of  the  passive  tracer  plume. 

'2  '  2 

Becker,  Hottel  and  Williams  (1965)  present  equations 
for  calcexT&ting  the  concentration  mean  and  variance  of  the 
contaminated  eddies  alone,  based  on  knowledge  of  the  inter- 
mittency  factor  and  the  overall  concentration  mean  and  variance. 
Using  the  subscript  E  to  indicate  the  property  of  a  tracer- 
containing  eddy,  the  mean  eddy  concentration  is 

c£  =  c/y  ,  (4.23) 

During  a  fraction  y  of  the  time,  the  instantaneous  excursion 
from  c  is  cE  -  c,  and  during  ( 1 -y )  it  is  simply  -c.  Thus 

_  2 

c ' 2  =  y (cE~c  )  2  +  (1 -y)c  , 

into  which  is  substituted  cE  =  cE  +  c  E  to  obtain 
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c^2  =  [c2  -  (1  -y)c2]/y  .  (4.24) 

These  equations  are  useful  for  combustion  problems  and  flows 

incorporating  fast  chemical  reactions,  applications  in  which 

the  eddy  concentrations  are  important. 

Example  profiles  of  measured  concentration  mean  c 
2 

and  variance  c*  have  been  abstracted  from  Chapter  V  and  are 

shown  in  Fig.  4.7.  Also  shown  are  the  corresponding  profiles 
_  2 

of  c^  and  c£  obtained  from  Eqs.  (4.23)  and  (4.24).  As  the 
plume  boundaries  are  approached,  the  eddy  concentration  fluctu¬ 
ations  achieve  a  uniform  distribution  of  concentration  mean 
and  variance,  with  the  fluctuation  intensity  becoming 

1 CE  c^  =  4.3.  This  uniform  fringe  region,  in  which 

concentration  variance  is  due  almost  entirely  to  i ntermi ttency , 
is  similar  to  that  observed  by  Becker,  Hottel  and  Williams 
(1966),  who  measured  a  smaller  asymptotic  value  of  i  =  0.73. 
The  different  fluctuation  intensity  may  be  due  to  lower  turbu¬ 
lence  levels  within  their  axisymmetric  jet/quiescent  fluid 
configuration,  as  compared  with  the  laterally  uniform  boundary 
layer  turbulence  at  plume  level  in  the  present  study. 

Returning  our  attention  to  Eq.  (4.19)  and  Fig.  4.5, 
it  is  apparent  that  the  Gaussian  error  function  which  fits 
data  above  the  plume  centerline  does  not  work  so  well  near  the 
ground.  This  is  probably  because  the  plume  boundary  fluctuations 
near  the  ground  are  constrained,  so  their  motions  do  not  follow 
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a  Gaussian  distribution.  Near  the  ground,  observed  variation 
of  y  with  x  and  z  is  the  result  of  complex  interactions. 

More  detailed  experiments  are  needed  before  the  i ntermi ttency 
model  (4.19)  can  be  further  refined  by,  for  example,  the  use 
of  a  truncated  Gaussian  position  distribution  for  the  lower 
plume  boundary. 
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x  Full  Scale  Downwind  Distance  (m) 


FIG.  4.2  Variation  of  Plume  Spread 

with  Distance  from  Source, 
for  4000 : 1  Scale  Factor 
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FIG.  4.3  Measured  Autocorrelations  of  Velocity  and 

Concentration  on  Plume  Centerline 


y/yV2  Dimensionless  Lateral  Distance 
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FIG.  4.4  Lateral  Concentration  Intermittency 

Profiles  at  Plume  Height 


Dimenionless  Vertical  Distance 
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FIG.  4.5  Vertical  Concentration  Intermittency 

Profiles  Through  Plume  Centerline 


Along  Plume  Centerline 
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y/o-y  Lateral  Distance  Off  Plume  Centerline 


y/a-y  Lateral  Distance  Off  Plume  Centerline 

FIG.  4.7  Concentration  Mean  and  Variance 

in  Dilute  Helium  Eddies 
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CHAPTER  V 


DETERMINISTIC  PROPERTIES  OF  THE 
CONCENTRATION  FLUCTUATIONS 


5. 0  I  ntroduct i on 

A  dynamically  neutral  material  released  continuously 
into  the  turbulent  atmospheric  boundary  layer  is  transported 
by  random  gusts  and  eddies,  so  that  the  instantaneous  material 
concentration  c(x,t)  downwind  of  the  source  also  becomes  a 
random  variable.  A  full  statistical  description  of  the  random 
concentration  fluctuations  is  necessary  to  identify  possibly 
hazardous  peak  concentrations,  and  this  requires  information 
about  the  concentration  probability  density  function  (pdf). 

The  intermittent  lognormal  pdf  model  to  be  developed  in  Chapter 

VI  is  expressed  in  terms  of  concentration  i n termi ttency  y, 

_  2 
time-mean  concentration  c,  and  concentration  variance  c'  . 

The  i n termi ttency  ,  mean  concentration  and  concentration  variance 
are  statistical  rather  than  physical  quantities.  However,  the 
instantaneous  concentration  c,  even  though  a  random  variable, 
is  actually  transported  and  dispersed  physically  by  the  turbu¬ 
lent  winds,  and  must  be  constrained  instantaneously  by  the 
equation  of  species  mass  conservation, 

=  -V- (uc-DVc )  (5.1) 
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where  c  is  the  concentration,  u^  the  random  wind  vector  and  D  the  molecular 
diffusi vity.  This  equation  is  useful  for  dispersion  in  a  turbulent  wind 
field  because,  as  will  be  shown,  it  leads  to  a  deterministic  description 
for  the  spatial  distribution  of  concentration  mean  c  and  variance  c1^. 

This  chapter  compares  experimental  data  with  predic¬ 
tions  from  the  well-known  reflected  Gaussian  solution  for  the 
mean  concentration  profile.  Then,  a  new  deterministic  model 
for  concentration  variance  is  developed  which  compares  well 
with  wind  tunnel  data.  This  concentration  variance  is  needed 
as  input  to  the  complete  statistical  model  of  concentration 
fluctuations  developed  in  Chapter  VI.  Finally,  key  assumptions 
made  in  previous  work  on  concentration  fluctuations  are  tested 
against  new  wind  tunnel  data,  and  their  regions  of  applicability 
are  defined. 

5  . 1  Present  State  of  Knowledge 

One  recent  theoretical  work  on  concentration  fluctu¬ 
ations  is  that  of  Csanady  (1967,  1973)  who  developed  an  analytic 
expression  for  the  variance  field  downwind  of  a  continuous 
point  or  line  source  in  a  uniform,  unbounded  flow  containing 
isotropic  turbulence.  Both  references  provide  a  good  review 
of  earlier  research  which  led  to  Csanady's  model.  A  key 
element  provided  by  Batchelor,  Townsend  and  Howells  (1959) 
was  the  observation  that  concentration  variance  is  a  transport¬ 
able  quantity,  and  can  be  transferred  and  dissipated  in  the 
same  way  as  turbulence  kinetic  energy.  Hinze  (1975,  p.  279) 
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contributed  the  concept  of  a  concentration  microscale,  and 
showed  that  self-similar  correlation  coefficients  require 
fluctuation  decay  rates  that  vary  inversely  with  travel  time. 
This  type  of  decay  was  observed  experimentally  by  Gibson  and 
Schwarz  (  1  963  )  . 

Csanady's  solution,  which  is  based  on  the  mass  con¬ 
servation  Eq.  (5.1),  requires  self-similarity  of  variance 
profiles,  and  contains  the  following  assumptions: 

•  the  flux  of  both  mean  concentration  and  concentration 
variance  may  be  expressed  as  the  product  of  an  eddy 
diffusivity  times  gradient; 

•  the  eddy  diffusivities  in  the  vertical  and  crosswind 
directions  are  equal  and  constant  for  a  given  cross 
section  of  the  plume,  although  they  vary  with  distance 
from  the  source; 

•  the  dissipation  rate  of  concentration  variance  is  propor¬ 
tional  to  the  local  concentration  variance; 

•  the  reciprocal  of  this  proport ional i ty  factor  is  a  'decay 
time  scale'  which  is  constant  for  a  given  section  but 
which  increases  linearly  with  distance  from  the  source. 

New  wind  tunnel  experiments  are  used  in  this  chapter  to  test 
the  last  two  of  these  assumptions,  which  are  found  to  be  valid 
in  the  wind  tunnel's  simulated  atmospheric  boundary  layer. 

The  measured  decay  time  scale  T  ^  increases  with  x  at  a  rate 
which  indicates  that  concentration  fluctuations  in  the  tracer 
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plume  are  dissipated  more  slowly  than  fluctuations  in  isotropic 
grid  turbulence. 

The  work  of  Kewley  (1978),  mentioned  in  Chapter  I, 
is  an  example  of  the  recent  application  of  fluctuation  theory 
to  the  analysis  of  fast  chemical  reactions  in  a  turbulent 
plume.  Kewley  makes  the  same  assumptions  as  Csanady,  with 
the  added  restriction  that  the  decay  time  scale  of  fluctua¬ 
tions  be  independant  of  distance.  In  order  to  simplify  the 

2 

solution  for  c'  ,  the  advection  and  gradient  flux  terms 
appearing  in  the  variance  Eq.  (1.3)  are  neglected  on  the 
assumption  that  they  are  small  compared  with  the  local  rates 
of  production  and  dissipation  of  fluctuations.  Then  the  local 
production  and  dissipation  must  be  in  balarice,  and  for 


K  ?*Kz,  c'2  is  given  by 


-  2 


-  2 


c'2  ■  2VKy(fy>  +  Kz{ff)  ]  • 


(5.2) 


Kewley  acknowledges  that  (5.2)  can  lead  to  unrealistic  profiles, 
and  proposes  an  empirical  adjustment  (sketched  in  Chapter  I) 
that  produces  flat  variance  profiles  on  the  plume  centerline. 

The  present  experiments  show  that  production  and 
dissipation  are  not  in  balance,  and  in  fact  that  with 

/V 

increasing  distance  from  the  source,  production  of  c1 
becomes  negl igible  compared  with  dissipation.  It  then 
follows  from  the  variance  balance  equation  (1.3)  that  the 
advection  and  gradient  flux  terms  must  be  of  the  same  order 
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as  the  dissipation  term,  and  cannot  be  neglected  as  Kewley 
suggests  . 

In  summary,  existing  solutions  for  the  variance  field 
apply  only  to  concentration  fluctuations  in  a  homogeneous, 
isotropically  turbulent  flow  with  uniform  mean  velocity,  far 
from  solid  boundaries  such  as  the  ground.  No  published 
theoretical  solutions  are  available  for  the  variance  field 
under  more  general  conditions. 

5 . 2  The  Governing  Equations 

The  diffusion  equation  (5.1)  is  used  as  the  starting 
point  in  the  analysis  of  concentration  fluctuations.  Separating 
the  turbulent  velocity  vector  and  concentration  scalar  into 
mean  and  fluctuating  components  u_=u+u_'  and  c  =  c  +  c',  then  taking 
ensemble  averages  of  both  sides  of  Eq.  (5.1),  gives  a  balance 
equation  for  the  mean  concentration. 


u-Vc  =  -V • u 1 c 1  +  DV2c  .  (5.3) 

a  L  — 

Subtracting  this  averaged  equation  from  the  original  equation 
provides  a  relationship  for  the  fluctuating  part  of  the  c-field. 
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■)  =  DV  c '  .  (5.4) 


Multiplying  this  equation  by  2c1,  then  taking  means  and 
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simplifying  by  use  of  the  continuity  equation  v*i[  =  0  results 
in  a  balance  equation  for  the  concentration  variance  c  '  ^ , 
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+  u_.^—  =  -2u  '  c 
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transient  advection  production  spatial  flux  dissipation 


where  (j)  is  a  dissipation  function, 
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Eq.  (5.3)  for  the  mean  concentration  field  is  simplified  by 
neglecting  axial  diffusion  (the  'slender  plume'  approximation) 
and  molecular  diffusion.  Furthermore,  we  assume  that  the 
Reynolds  flux  -u  '  c '  is  proportional  to  the  mean  concentration 
gradient,  i . e . 


u  .  c 


=  K 


3c 


i  3x 


(5.7) 


where  the  proporti onal i ty  factor  Ki  (the  'eddy  di ffusi vi ty ' ) 
is  a  function  of  x,  y  and  z.  Then,  for  steady  state  (^  +  0) 
and  a  unidirectional  mean  wind  (u-j  =  U,  u  2  =  u  3  =  ’  Eq.  (5.3) 

can  be  written  in  the  familiar  form 
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where 

U  =  U(z)  , 

Ky  =  Ky(x’z)’ 

Kz  =  Kz(x,y,z), 

c  =  c  (x  ,y  ,  z  )  . 

The  advection  term  on  the  left-hand  side  of  (5.8)  is  balanced 
by  the  gradient  flux  terms  on  the  right-hand  side.  The  wind 
field  is  assumed  to  be  laterally  homogeneous,  so  that  K  is 
constant  in  the  crosswind  direction. 

If  the  same  simplifying  assumptions  are  applied  to 

o' 

the  balance  equation  for  concentration  variance  c‘  ,  the  result 

2 

is  (with  Sec1  for  notational  simplicity), 


U 


9S_ 

9x 


K'.A  +  |I(K;|f}  +  2Ku(|f)+  2K,(f§)- 


y3y2 


y  ay 


z  '  9  z 


6D / 9c  1 
2l9t 


U 


advection  gradient  flux  production  dissipation 

where  the  eddy  diffusivities  for  S  are  defined  by 


(5.9) 


,  ,2  „  ,  3S 

-UiC  =  Ki 3x 


1 


(5.10) 


The  dissipation  term 


■  2D^|f)  +  ( 


3C')  +  (ffV] 
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6D  /  9c 
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has  been  simplified  by  assuming  locally  isotropic  concentra¬ 
tion  fluctuations,  with  Taylor's  hypothesis  —  =  1  replacing 
the  spatial  derivative.  The  functional  dependencies  are 


D  =  constant  , 

U  =  U  ( z  )  , 

Ky  =  KyU.z), 

K;  =  k;(x,Z), 
c  =  c(x,y,z), 
c  1  2  =  S  =  S(x,y,z) , 

Kz  =  Kz(x,y,z) , 

Kz  =  K^(x,y,z), 

Sr1  ^  Sr'  ^ 

(ft-)  =  {br)  (x-y>2)  • 


Equation  (5.9),  rewritten  in  polar  coordinates,  has  been 
solved  by  Csanady  (  1  967  )  for  the  case  of  K  =Kz=K  =K  =K(x)  and 
<|>c  proportional  to  S,  where  the  proportionality  factor  T^ 
grows  linearly  with  downwind  distance  x.  Later  in  this 
chapter,  new  experimental  data  is  used  to  show  that  this 
proportionality  assumption  on  <f>  is  a  very  good  one,  even 
with  complicated  boundary  conditions  on  Eq.  (5.9)  such  as  occur 
for  plumes  with  ground  level  reflection  of  the  mean  concentration 


'  1 
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profile  and  surface  dissipation  of  the  concentration  fluctuations. 

Eq.  (5.9)  governs  the  spatial  distribution  of  concen- 

tration  variance  c  downwind  of  a  point  or  line  source  in  a 

turbulent  atmospheric  boundary  layer.  The  terms  of  Eq .  (5.9) 

can  be  interpreted  in  the  same  way  as  for  the  turbulence 

kinetic  energy  equation:  the  advection  of  c1^  on  the  left  is 

balanced  on  the  right  by  the  gradient  flux  terms,  plus  the 

difference  between  local  production  and  dissipation  rates  of 
2 

c'  •  If  if  were  not  for  this  last  difference  term,  the  govern¬ 
ing  equation  (5.9)  for  concentration  variance  would  have  the 
same  form  as  Eq.  (5.8)  which  governs  the  spatial  distribution 
of  mean  concentration.  For  the  case  of  similar  boundary 
conditions,  any  differences  between  the  spatial  distribution 
of  concentration  variance  and  that  of  mean  concentration 
must  result  from  the  inclusion  of  this  net  dissipation  term. 

The  present  measurements  will  show  that,  for  an  elevated 
point  source  emitting  a  passive  contaminant  continuously 

into  a  turbulent  atmospheric  boundary  layer,  production  of 
2 

c'  becomes  negligible  compared  with  dissipation  as  the  plume 

moves  downwind.  The  resulting  simplification  of  the  variance 

equation  (5.9),  combined  with  a  judicious  choice  (see  Sec.  5.7)  for  the 

x-gradient  of  decay  time  scale  Td,  produces  an  equation  having 

the  same  form  as  that  for  mean  concentration.  With  different 

boundary  conditions  on  S,  this  leads  to  a  theoretical 

expression  for  the  spatial  distribution  of  concentration  var- 

2 

iance  that  agrees  very  well  with  the  measured  profiles  of  c1  • 
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5 . 3  Wind  Tunnel  Measurements 

Vertical  and  lateral  profiles  of  mean  concentrati on  c, 

2 

concentration  variance  c'  and  time  derivative  variance 

were  obtained  at  various  distances  4.3  -  x/H  -  19.2 
downwind  of  an  elevated  'point*  source  in  a  simulated  atmos¬ 
pheric  boundary  layer.  Details  of  the  measurement  methods 
are  given  in  Chapters  II,  III  and  IV.  This  data  was  analysed 
in  several  ways,  and  led  to  the  formulation  of  a  theoretical 
model  for  the  spatial  distribution  of  concentration  variance. 

5.3.1  Mean  Concentration 

The  first  step  was  to  see  if  the  reflected  Gaussian 

model  for  c  provided  a  good  description  of  the  observations. 

The  equation 

c  =  U  -)exp(^?)(exp[~(z'2)  ]  +  exp[--(z  +  !p  ]} 

MhVz  2a„  2aZ  2af 

y  z  z 

(5.12) 

is  widely  used  to  calculate  mean  concentrations  downwind  of  an 
elevated  point  source  emitting  continuously  into  an  atmospheric 
boundary  layer.  In  this  equation,  Q  is  the  pollutant  emission 
rate,  Uu  is  the  wind  speed  at  source  height  H,  y  and  z  are  the 
lateral  distance  from  plume  centerline  and  height  above  ground, 
respectively,  and  o  and  cr^  are  x-dependant  standard  deviations 
of  the  (assumed  Gaussian)  lateral  and  vertical  plume  concentra¬ 
tion  profiles.  Veigele  and  Head  (1978)  have  solved  the  govern¬ 
ing  mean  field  equation  (5.8)  for  bounded  flows,  and  Eq.  (5.12) 
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is  the  result.  Prior  to  their  derivation,  justification  for 
the  widespread  use  of  Eq.  (5.12)  rested  solely  on  its  good 
performance  in  predicting  measured  values  of  cf,  plus  the 
knowledge  that  the  reflected  Gaussian  is  a_  solution  of  the 
mean  field  equation  for  the  simplified  boundary  conditions  of 
uniform  flow  over  an  impermeable  surface. 

A  least-squares  error  fit  of  Eq.  (5.12)  to  normalized 

wind  tunnel  measurements  of  c  yielded  the  values  of  a  and  a 

y  z 

(discussed  in  Chapter  IV)  which  appear  in  brackets  in  Table 
5.1.  These  best-fit  values  have  been  used  to  generate  the 
theoretical  solution  (5.12)  for  c  at  the  downwind  locations 
where  measurements  were  taken,  and  the  theory  and  observations 
are  compared  in  Figs.  5.1  through  5.7.  Fig.  5.1  shows  that 
the  observed  centerline  decrease  in  mean  concentration  is  well 
reproduced  by  the  reflected  Gaussian  model.  The  observed 
lateral  distribution  of  c  is  also  well  represented  by  the 
Gaussian  model,  as  shown  by  Fig.  5.2.  Theoretical  vertical 
profiles  of  mean  concentration  are  compared  with  observations 
in  Figs.  5.3  to  5.7.  The  agreement  is  quite  good  over  the 
vertical  profile  as  a  whole,  but  near  the  ground,  differences 
as  large  as  40%  are  seen.  The  reflected  Gaussian  theory 
appears  to  first  overestimate  and  then  underestimate  the  ground 
level  concentrations  as  the  downwind  distance  increases.  Note, 
however,  that  the  prediction  error  is  not  nearly  as  large  as 
the  deviations  of  200  to  300  percent  often  quoted  in  the 
literature  (e.g.  Weil  and  Jepsen  (1977))  as  the  best  performance 
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to  be  expected  of  the  Gaussian  model.  These  large  error 
estimates  occur  in  full-scale  field  studies,  where  the  input 
parameters  Q,  U^,  and  cr^  are  often  not  accurately  known, 
and  factors  such  as  slow  wind  direction  shifts  play  an 
important  role.  Wind  tunnel  experiments,  because  of  their 
more  accurately  measured  and  steadier  conditions,  provide  a 
good  indication  of  what  to  expect  from  the  Gaussian  model 
under  ideal  field  conditions:  the  relative  error  in  mean 
concentration  predictions  is  less  than  10%  of  the  local  con¬ 
centration  except  at  the  plume  edges  and  near  the  ground, 
where  errors  of  20-40%  are  sometimes  observed. 

5.3.2  Concentration  Variance 

One  important  goal  of  the  wind  tunnel  experiments 

was  to  measure  the  spatial  distribution  of  concentration 

variance  in  an  elevated  plume.  The  observed  axial,  lateral 

2 

and  vertical  variations  of  c'  ,  the  root- mean-square  (rms) 
concentration,  are  plotted  in  Figs.  5.1  through  5.7.  The 
most  noticeable  difference  between  the  rms  and  mean  con¬ 
centration  fields  is  that  the  rms  concentration  decreases 
much  faster  near  the  ground  than  does  the  mean  concentration. 
This  is  a  very  important  observation,  because  it  supports  the 
idea  that  not  only  the  vector  fluctuations  of  turbulence, 
but  also  the  seal ar  fluctuations  undergo  a  net  dissipation 
near  solid  boundaries.  The  removal  of  scalar  fluctuations 
near  a  solid  surface  is  physically  reasonable  because  in  the 
thin,  highly  sheared  surface  layer,  any  concentration  gradients 
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tend  to  be  smoothed  out  by  molecular  diffusion.  In  the 
absence  of  mean  concentration  gradients,  there  is  no  produc¬ 
tion  of  new  fluctuations  to  replace  those  which  are  dissipated, 
s o  n e t  dissipation  is  large.  However,  the  extent  to  which 
this  surface-induced  dissipation  takes  place  is  difficult  to 
measure.  For  example,  all  full-scale  concentration  field 
studies  show  time-varying  concentrations  at  heights  only  one 
or  two  meters  above  the  ground,  which  indicates  that  surface 
dissipation  of  fluctuations  cannot  be  total.  On  the  other 
hand,  concentration  i n termi ttency  is  an  important  component 

of  observed  full-scale  concentration  time-series,  and  Eq.  (4.24)  shows 

that,  in  the  presence  of  intermi ttency ,  a  minimum  concentration  var- 

iance  of  c'  =  y  (l-y)c^  will  always  be  observed,  even  when  fluctuations 
2 

c'  £  in  the  non-zero  concentration  are  reduced  to  zero. 

With  the  fast  response  sensor  used  in  this  study,  it 
was  not  possible  to  extend  the  measurements  closer  than 
0.5  cm  (.17  source  heights)  to  the  surface,  so  the  dissipative 
surface  layer  could  not  be  observed  directly.  The  presence 
of  high  dissipation  at  the  surface  was  deduced  from  the  observed 
decrease  of  scalar  fluctuations  as  the  surface  was  approached. 

In  the  theoretical  analysis,  the  simplest  assumption  of  com¬ 
plete  surface  dissipation  (zero  fluctuations)  was  used.  The 
conclusions  drawn  from  the  subsequent  data  analysis  are  not 
sensitive  to  this  assumption. 

The  effect  of  surface-induced  dissipation  can  be 
thought  of  as  opposite  to  the  solid  boundary  reflection  exper¬ 
ienced  by  the  mean  concentrations.  The  following  sketch 
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illustrates  the  concept  of  surface-induced  dissipation;  it 

is  similar  to  the  "absorption"  boundary  condition  used  by 

Csanady  (1973,  p.40)  to  describe  deposition  of  aerosols  on  a 

vertical  wall.  An  obvious  question  derives  from  the  observed 
2 

c'  profiles: 

2 

How  well  can  the  observed  distribution  of  c‘ 
be  reproduced  by  a  Gaussian  plume  model  with 
dissipation  rather  than  reflection  at  the  wall? 
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The  rms  concentration  (instead  of  c'^)  was  selected  for 
plotting  in  Figs.  5.1  through  5.7  in  order  to  have  the  same 
physical  units  for  the  mean  and  fluctuating  concentration 
fields.  The  plots  do  not  show  that  lateral  distributions 
of  c1  appear  to  be  as  Gaussian  as  the  mean  concentration. 
Combining  this  fact  with  the  observations  from  Chapter  III 

2  p 

showing  c  proportional  to  (Q/Ur  (which  was  assumed  by 
Csanady  (1967)),  the  following  form  is  suggested  for  the 
variance  field: 


s  ( _ ^ 

So4ttU  ,,a'aT 


2  -v2 

)  exp  (~^— 


H  y  z 
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•)  (exp 


- ( z-H ) 

2a'  2 
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-(z+H) 

2a'  2 
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(5.13) 

where  sq  is  a  'source  variance'  parameter  whose  spatial  var¬ 
iation  must  be  determined  from  experiment.  The  parameter  sq 
is  analogous  to  the  normalized  concentration  source  variance 
g ( 0 )  which  appears  in  Csanady's  (1967)  solution  for  concen¬ 
tration  fluctuations.  The  last  exponential  term  in  Eq.  (5.13) 
is  the  surface  dissipation  component;  the  leading  negative  sign  serves 
to  differentiate  it  from  the  positive  reflection  term  in  Eq. 
(5.12)  for  mean  concentration.  The  dissipation  term  is  the 
result  of  a  negative  image  source  underground,  and  in  effect 
creates  a  distributed  surface  'sink'  which  exactly  cancels 
all  fluctuations  reaching  the  surface.  This  surface  sink 
represents  the  dissipation  of  fluctuations  taking  place  in  a 
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very  thin  layer  near  the  surface.  Surface  dissipation  of  fluctuations  is 
consistent  with  the  reflected  Gaussian  prediction  of  —■  +  0  as  z  -*  0  so 

o  Z  * 

that  the  surface  smooths  out  concentration  gradients,  at  least  in  the 
vertical . 

Note  that  total  dissipation  is  assumed  in  Eq.  (5.13),  so  that 
concentration  fluctuations  are  forced  to  zero  at  the  ground.  This  assump¬ 
tion  is  probably  too  restrictive,  but  the  observations  with  which  the 
dissipated  Gaussian  model  is  compared  do  not  provide  sufficient  detail  to 
distinguish  between  total  and  partial  dissipation.  A  way  to  incorporate 
partial  dissipation  into  the  model  will  be  suggested  in  the  next  section. 

A  least-squares  error  fit  of  this  dissipated  Gaussian  model  to 
the  normalized  wind  tunnel  measurements  of  c'2  yielded  the  values  of  a 

,  y 

and  az  given  in  Table  5.1  (the  corresponding  values  of  a and  a z  for  the 
mean  field  are  in  brackets). 


TABLE  5.1 

STANDARD  DEVIATIONS  FOR  CONCENTRATION  PROFILES 
OF  (MEAN)  AND  VARIANCE 


X 

a' 

(a) 

a* 

(a  ) 

(cm) 

y 

(cm)  y 

Z  (cm) 

z 

13.0 

1  .22 

(1.23) 

0.92 

(1.00) 

19.0 

1  .67 

(1.72) 

1  .39 

(1.41) 

27.5 

2.29 

(2.32) 

1  .98 

(1.97) 

37.5 

2.89 

(2.97) 

2.34 

(2.35) 

57.5 

4.21 

(4.08) 

2.70 

(2.91  ) 

The  first  point  to  note  is  that  Eq.  (5.13)  is  in  fact  able  to 
reproduce  the  observed  distribution  of  concentration  variance 
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and,  as  seen  from  Figs.  5.1  through  5.7,  it  fits  the  c '  data 
at  least  as  well  as  the  reflected  Gaussian  model  (5.12)  fits 
the  c  measurements.  The  second  observation  is  that,  within 
experimental  accuracy,  the  lateral  and  vertical  standard 

p 

deviations  appearing  in  the  dissipated  Gaussian  model  for  c‘ 
are  the  same  as  those  in  the  reflected  Gaussian  model  for  c. 
This  is  remarkably  fortunate,  because  it  means  that  the  large 
body  of  literature  dealing  with  the  spreading  rates  of  mean 
concentration  fields  can  also  be  applied  to  calculations  of 

the  fluctuating  concentration  field. 

After  the  dissipated  Gaussian  model  (5.13)  was 

2 

proposed  and  tested  against  the  c'  data,  it  was  discovered 
that  the  model  could  be  derived  theoretically  from  the 
balance  equation  (5.9).  The  theoretical  derivation,  presented 
in  Section  5.7,  is  based  in  part  on  the  assumption  that 
K^  =  Ky,  K^  =  Kz,  from  which  it  follows  that  ay  =  ay»  °z  =  az  must  be 
satisfied. 

The  'source  variance'  parameter  sq  was  determined  by 
dividing  Eq.  (5.13)  by  the  square  of  Eq.  (5.12)  and  rearranging 
to  obtain 


•  exp  (^^2 
2a 

y 


(exp[-(z-H)2/2az2]  +  exp [- ( z+H ) 2/ 2az 2] } 2 
{exp[-(z-H)2/2az2]  -  exp[-(z+H)2/2a^2] } 


> 


(5.14) 
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where  a ^  and  a ^  -  az  based  on  the  data  in  Table  5.1. 

Here  ic  is  the  intensity  of  concentration  fluctuations, 
defined  by 


(5.15) 


c 

Measured  values  of  i  at  128  receptor  points  in  the  wind  tunnel  were  used 
to  determine  sQ  from  Eq.  (5.14);  60%  of  the  data  came  from 

the  five  vertical  profiles,  and  the  remaining  40%  from  the 
five  corresponding  lateral  profiles.  The  calculated  values 
of  sq  showed  some  scatter  (^  15%)  but  no  significant  spatial 
dependence  trends  were  obseryed.  The  average  value  of  sQ  was 


s  =  2 . 6  ±  0.4  • 
o 

The  results  of  computing  s  using  observed  values  of  i  were 

W  C 

very  gratifying:  the  value  of  sq  was  found  to  be  independant 
of  receptor  position  within  the  plume.  This  supports  the 
interpretation  of  sq  as  a  measure  of  the  Source  variance* 
strength  which,  for  a  passive  tracer  released  at  the  mean  wind 
speed,  is  most  likely  proportional  to  the  level  of  atmospheric 
turbulence  near  the  source.  Recent  work  by  Robins  and 
Fackrell  (1979)  suggests  that  the  magnitude  of  sq  may  depend 
on  the  ratio  of  source  diameter  d$  to  the  integral  scale  of 
turbulence  A  near  the  source.  Further  experimental  work  is 
needed  to  quantify  the  dependance  of  sq  on  atmospheric  tur¬ 
bulence  and  source  characteristics. 
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5.3.3  Concentration  Fluctuation  Intensity 

The  most  useful  application  of  the  dissipated  Gaussian 
model  is  to  generate  the  values  of  concentration  fluctuation 
intensity  ic  required  as  input  by  the  statistical  concentration 
model  developed  in  Chapter  VI.  This  statistical  model  uses 
local  values  of  concentration  i ntermi ttency  factor  y  and  in¬ 
tensity  of  fluctuations  i  to  predict  the  peak-to-mean  concen¬ 
tration  ratios  at  any  point  downwind  of  an  air  pollution 
source.  The  fluctuation  intensity  has  been  previously  defined 

as  .  _  /  ,  2  * — 

1(;  =  /  C  /C. 

Fig.  5.1  and  Figs.  5.8  through  5.13  illustrate  the  ability 

2 

of  the  dissipated  Gaussian  model  for  c'  to  generate  accurate 

values  of  i  when  combined  with  the  reflected  Gaussian  model 
c 

for  c.  The  agreement  between  Eq .  (5.16),  obtained  by  dividing 
Eq.  (5.13)  by  the  square  of  Eq.  (5.12), 

o  2  {exp[-(z-H)2/2cr  2]  -  exp[- (z  +  H ) 2/2a  2] } 

i  =  s  exp(  ^  *) - 5 - 5 - 9 - 5 — 7  » 

c  0  2 a  (exp[-(z-H)  /2a  ]  +  exp[-(z+H)  /2 a  ]} 

y  z  z 

(5.16) 

and  experiment  is  excellent.  In  fact,  the  agreement  between 
(5.16)  and  experiment  is  better  than  the  separate  agreement 

of  either  of  its  components  c1  and  c.  This  may  be  due  in 

“7? 

part  to  systematic  theory-induced  errors  in  calculating  c 
with  (5.13)  being  partially  cancelled  by  similar  systematic 
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errors  in  the  calculation  of  c  using  (5.12).  Systematic 
experimental  error  is  unlikely  because  measurements  of  c'* 2 
were  obtained  with  an  instrument  different  from  that  used  for 
mean  concentration  c. 

2 

One  aspect  of  the  suggested  Gaussian  model  for  c‘  is 
that  it  predicts  a  complete  absence  of  concentration  fluctua¬ 
tions  at  ground  level.  For  purposes  of  hazard  assessment,  it  is 
useless  to  estimate  surface  fluctuation  intensity  by  placing  z=0  in 
Eq.  (5.16).  This  may  not  be  a  serious  deficiency, 
because  it  is  unclear  if  any  receptors  of  practical  interest 
occur  exactly  at  z  =  0.  Fig.  5.9  shows  extremely  high 
gradients  of  i  near  the  ground,  which  indicates  that  receptors 
just  above  ground  level  can  be  exposed  to  significant  concen¬ 

tration  fluctuations.  However,  it  may  be  that  an  empirical 
dissipation  coefficient  "d"  should  be  introduced  in  the 
dissipated  Gaussian  model,  i.e. 


2  2 

S  =  S  (——2 - )  exp(^-^-){exp[-(z-H)2/2a  2]  -  d  •  exp[-(z+H)2/2cr  2]}  > 

o  '2rrU  a  az  2a  2  z  z 

J  y 

(5.17  } 

where  d  is  a  number  less  than  unity.  Additional  wind  tunnel 
measurements  of  c  1  2  at  and  near  ground  level  are  needed  to 


determine  d. 
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5 . 4  Mean  Field  Diffusivity 

It  is  possible  to  derive  values  for  and  Kz  from 
the  experimental  data.  In  addition,  one  can  substitute  the 

reflected  Gaussian  model  into  the  mean  field  equation  (5.8). 
This  allows  determination  of  the  vertical  concentration 
diffusivity  Kz  which,  if  used  in  the  mean  field  equation 
(5.8),  would  give  the  same  concentration  profiles  as  the 
reflected  Gaussian  equation  (5.12).  The  form  of  this  eddy 
diffusivity  profile  is  of  interest,  because  it  may  help  to 
define  the  link  between  'K-Theory'  and  'Gaussian  Model' 
approaches  to  the  solution  of  plume  dispersion  problems 
under  non-u ni form  wind  conditions. 

In  calculating  plume  dispersion,  a  K-theory  solution 
for  the  mean  concentration  field  is  usually  obtained  by 
numerically  integrating  Eq.  (5.8)  for  a  given  eddy  diffusivity 
field.  A  general  specification  of  the  diffusivity  field  is 
not  possible,  because  and  Kz  are  not  true  properties  of 
the  turbulent  fluid.  They  depend  on  the  complete  source/flow/ 
boundary  configuration,  so  their  values  are  quite  problem- 
specific.  Because  of  the  resulting  lack  of  data  on  species 
diffusivity,  a  common  approximation  is  to  replace  the  unknown 
species  eddy  diffusivity  Kz  with  the  momentum  diffusivity  Km 
determined  from  the  turbulent  velocity  field.  This  approach 
ignores  the  fact  that  Kz  depends  on  the  source/boundary  config¬ 
uration,  a  behavior  which  appears  explicitly  in  the  'reflected 
Gaussian'  eddy  diffusivity  derived  in  this  section. 


' 
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Values  for  the  vertical  and  lateral  eddy  diffusivity 
were  derived  from  the  present  measurements  as  follows:  First, 
interpolation  and  numerical  differentiation  of  data  were  used 
to  obtain  values  for  U ,  c,  —  and  ^  at  2  mm  height  increments 
for  each  of  the  five  vertical  profiles  through  the  plume  center 
line.  Then,  to  simplify  the  data  analysis,  and  because  measure 
ments  of  lateral  flux  -v'c'  were  not  available,  the  1  a tera 1 
eddy  diffusivities  and  derivatives  of  c  were  calculated  using 
a  Gaussian  distribution  for  the  lateral  profiles,  so  that 


* 


3c 

3y 


•v 


(5.18) 


92c  = 

77 


(^7  -  1) 


*  on  plume  centerline  . 


The  vertical  eddy  diffusivities  were  derived  from  experimental 
measurements  by  numerical  integration  of  the  terms  in  Eq. 

(5.8)  to  obtain 


K. 


1 


3c/  3x 


<ulf  - 


K. 


32c 

3y2 


)  dz 


(5.19) 


o 


The  experimentally  derived  value  of  vertical  Reynolds  flux  is 
provided  by  (5.19)  in  combination  with  (5.7), 
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-w '  c  1 


r  Z 

(U  f£  -  K  ^f)  dz  . 

3x  y  ay  2 

0 


(5.20) 


Examples  of  the  eddy  diffusi  vi  ties  and  Reynolds  fluxes  derived 
from  experiment  are  shown  in  Figs.  5.14  through  5.18. 

Using  the  reflected  Gaussian  model  for  mean  concentra¬ 
tion,  the  expressions  required  for  insertion  in  Eqs.  (5.19) 
and  (5.20)  are,  for  the  power  laws  a =  axP  and  gz  =  bxq  shown 
to  be  valid  in  Chapter  IV,  numbers  (5.12)  and  (5.18)  plus  the 
alongwind  gradient  for  a  reflected  Gaussian, 


1)  + 


-{• 


[(z-H)2-az2]exp[-(z-H)2/2az2] 


a  2  exp[  -  ( z-H  )  /2a_  ]  +  exp[-(z  +  H)  /  2  cr  2  ] 


[ ( z  +  H  )  2  -  az2]exp[-(z  +  H)2/2az2] 
exp[-(z-H)2/2a  2]  +  exp [ - ( z  +  H ) 2/ 2 a  2] 


}] 


(5.21) 


and  the  vertical  gradient  of  the  reflected  Gaussian, 

-  -(z-H)exp[-(z-H)2/2az2] 

3z  az2  exp[-(z-H)2/2az2]  +  exp[ - (z+H  ) 2/ 2az2] 

(z+H)exp[-(z+H)2/2az2] 

}  . 


exp[-(z-H)2/2az2]  +  exp[- ( z+H ) 2/ 2 az2] 


(5.22) 


1 18 


Substituting  these  expressions  into  Eq.  (5.19)  and  rearranging 
gives  the  'reflected  Gaussian'  eddy  diffusivity  of  mean  concen¬ 
tration 


K. 


UH  da 


-{• 


M  +  1 


dx  -(z-H)exp[-(z-H)2/2az2]  -  (z+H ) exp [ - ( z+H ) 2/ 2 az2 ] 


} 


(5.23) 

A  similar  operation  on  Eq.  (5.20)  gives  the  vertical  Reynolds 
mass  flux  for  a  reflected  Gaussian  concentration  field. 


-w '  c  ' 


Q  da  /dx  o  o 

- (I.  +  ljexp(-y  /2a  ) 

2*Vz 


(5.24) 


where  the  integrals  1-j  and  I2  are  given  by 


z 

I-|  =  jj— exp[-(z-H)2/2az2][(z-H)2/az2-  1  -  (1 -y2/a^2)|^(l -^p)]dz 

H 

■  0 

(5.25) 

and 


■z 

{f-exp[-(z+H)2/2a  2][(z  +  H)2/oz2 

UH 
•  0 


i  -  (i-y2/oy2%0-jji)]dz  • 


(5.26  ) 


For  uniform  wind  conditions  U  =  UR  =  UQ,  I-,  and  Ig  can  be 
integrated  in  closed  form  to  yield 
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and 


K 


z 


U  da  2 
0  z 

2  dx 


(5.  27) 


-w  c 


c 


Z  , 


(5.  28) 


where  z  =  0  on  the  plume  centerline.  Note  that  (5.27), 
derived  originally  by  Batchelor  (1949),  is  valid  for  a 
bounded  or  unbounded  flow  in  a  uniform  wind. 

The  theoretical  expressions  (5.23)  for  Kz  and  (5.24) 
for  Reynolds  flux  -w'c'  are  compared  with  experimental  values 
in  Figs.  5.14  through  5.18.  The  experimental  values  were 
obtained  by  numerical  integration  of  Eqs.  (5.19)  and  (5.20), 
with  experimental  measurements  used  to  compute  the  integrand 
terms.  The  agreement  between  the  Gaussi an-deri ved  Kz  theory 
and  the  experiments  is  qualitatively  quite  good,  with  the 
theoretical  eddy  diffusivity  showing  the  observed  transition 
from  ^-shaped  to  -^-shaped  curves,  but  requiring  a  slightly 

longer  travel  time  before  the  centerline  discontinuity  changes 
character  to  match  observations  farther  downwind.  The  good 
fit  between  experimental  and  theoretical  profiles  of  mean 
concentration  (Figs.  5.1  to  5.7)  is  evidence  that  the  solution 
of  equation  (5.8)  is  relatively  insensitive  to  the  values 
used  for  Kz-  The  mathematical  reason  for  this  insensitivity 
is  that  the  mean  concentration  depends  on  the  spatial  i ntegra 1 


of  eddy  diffusivity,  i.e.  c  at  any  receptor  depends  on  the 
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entire  upstream  Kz  field,  and  the  effect  of  details  in  the  Kz 
field  is  obscured.  Thus  the  distance  at  which  the  eddy 
diffusivity  'flips'  from  ^-shaped  to  —^-shaped  curves  is  not 
of  major  importance.  This  insensitivity  to  Kz  profiles  helps 
to  explain  why  reasonable  c  profiles  can  be  obtained  even  for 
d i f f u s i v i t i e s  which  are  independant  of  source  height  (Kumar, 

1  978)  . 

A  surprizing  aspect  of  the  diffusivity  curves  in  Figs.  5.14  to 
5.17  is  that  both  experimental  and  theoretical  profiles  show  a  discontinuity 
near  the  source  height.  It  is  probable  that  the  discontinuity  is  nothing 
more  than  an  artifact  produced  by  accumulated  experimental  and  numerical 
errors.  However,  this  discontinuity  vanishes  from  the  theoretical  profile 
(5.23)  when  the  wind  field  becomes  uniform,  because  then  the  z-gradient 
denominator  in  Eq.  (5.23)  goes  to  zero  at  the  same  point  that  the  numera¬ 
tor  becomes  zero;  at  that  point,  the  ratio  0/0  has  the  value  unity  and,  as 
expected,  Kz  takes  the  constant  value  given  in  Eq.  (5.27).  The  discontin¬ 
uous  theoretical  solution  for  Kz  occurs  only  under  non-uniform  wind  con¬ 
ditions,  in  which  case  the  vertical  concentration  gradient  (denominator  of 
K  )  inconveniently  goes  to  zero  while  the  Reynolds  flux  (numerator  of  Kz) 
is  still  different  from  zero.  This  theoretical  non-uniform  wind  phenomenon 
is  an  alternative,  but  less  plausible,  explanation  for  the  discontinuity 
found  in  the  experimental  eddy  diffusivities. 

The  theoretical  eddy  diffusivity  equation  (5.23)  has 
limited  practical  application,  since  it  provides  no  more  infor¬ 
mation  about  the  mean  field  than  is  already  contained  in  the 
more  convenient  reflected  Gaussian  model  for  c.  The  theory's 


. 
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main  feature  is  that  it  illustrates  how  the  eddy  diffusivity 
of  concentration  in  non-uniform  winds  should  depend  on  the 
source  location  within  the  wind  field  as  well  as  the  effective 


diffusion  length  scales  and  az«  The  momentum  diffusivity 
is  often  used  in  place  of  the  scalar  eddy  diffusivity,  but 
depends  on  the  integral  turbulence  length  scale  a  ,  which  was 
shown  in  Chapter  IV  to  be  larger  than  the  observed  for 
scalar  concentrations.  The  difference  is  due  to  unequal  'ages' 
for  velocity  and  scalar  fluctuations;  most  of  the  vertical 
plume  spread  is  complete  before  the  scalar  fluctuations  are 
old  enough  to  be  influenced  by  the  low  frequency  turbulence 
contributing  to  A  .  This  fact  is  not  widely  appreciated  and 
is  often  ignored,  especially  when  the  Reynolds  analogy  is 
used  to  equate  scalar  diffusivity  with  momentum  diffusivity. 
Because  the  momentum  diffusivity  does  not  depend  on  any  source 
location  characteristics,  whereas  K  z  explicitly  does,  the 
analogy  is  valid  only  in  the  asymptotic  stage  (a  ^  x, 

K  independent  of  x)  of  plume  dispersion  in  uniform  winds. 

5 . 5  Variance  Production  and  Dissipation 

The  variance  equation  (5.9)  is  similar  to  the  mean 
field  equation  (5.3)  in  that  advection  and  gradient  flux  terms 
are  common  to  both  equations.  However,  Eq.  (5.9)  contains  two 
additional  terms  which  describe  the  local  production  and 
dissipation  of  concentration  fluctuations.  The  importance  of 
these  extra  terms  is  examined  in  this  section. 
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Csanady's  (1967)  solution  for  c'  produces  lateral 
profiles  which  have  a  saddle  at  the  plume  centerline  if  the 
ratio  of  dissipation  <|>  to  diffusion  SK  /a  is  large.  High 

«y  7 

diffusion  and  low  dissipation  smooths  out  the  of f -center  1 i ne 
peaks  which  occur  at  the  production  maxima,  and  results  in 
single-peaked  distributions  much  like  those  observed  in  the 
present  experiments.  The  relative  importance  of  production 
and  dissipation  in  the  simulated  plume  was  established  by 
calculating  their  values  from  measurements  of  c  and  (y^— )  • 
Interpolation  and  numerical  differentiation  of  data  from  the 
76  vertical  profile  receptor  locations  were  used  to  derive 
values  for  U,  and  (It-)  2  at  2  mm  height  increments  for  each 

o  Z  at' 

of  the  five  vertical  profiles.  The  experimental  mean  field 

diffusivities  derived  in  Section  5.4  were  used  for  Kz. 

Section  5.2  has  shown  that  the  eddy  dissipation  function 

<h  can  be  related  to  the  concentration  time  derivative  variance 
c 

if  local  isotropy  is  assumed  and  Taylor's  hypothesis  is  used. 


The  result  was 


6D  ,  9c  1  x 

T  r*  9  '  9  +  / 


u 


9t 


(5.29) 


In  this  case  D  is  the  molecular  diffusivity  of  helium  in  air, 
equal  to  0.65  cm2/sec.  An  averaging  time  of  100  seconds  was 
used  to  obtain  ^  3%  repeatability  for  measurements  of  the  time 
derivative  variance;  the  observed  values  were  corrected  for 
probe- induced  dissipation  as  described  in  Chapter  III. 
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The  variance  production  term  is,  from  Section  5.2, 


2K  (|£-) 
y  3y' 


+  2Kz  <H>2- 


In  the  vertical  (x,z)  plane  through  the  plume  centerline,  the 
lateral  derivative  is  zero,  so  that  production  in  the 
vertical  centerline  plane  is  calculated  from 


Pv  (y=o) 


(5.30) 


The  variance  dissipation  <J) c  and  variance  production  P  v 
in  this  vertical  plane,  plus  net  dissipation  <J>  ^  =  $  c  -  P  v , 
have  been  calculated  from  the  measurements,  and  are  plotted 
against  height  in  Figs.  5.19  through  5.23.  To  a  large  extent, 
the  curves  behave  as  expected: 


§  Production  is  low  at  the  plume  centerline  and  plume 
fringes,  where  concentration  gradients  are  small.  The 
largest  production  occurs  where  mean  concentration 
gradients  are  maximum,  in  a  region  defined  by  a  conical 
surface  increasing  in  diameter  downwind,  with  apex  at 
the  source.  At  larger  distances,  where  the  plume  has 
impinged  on  the  ground,  this  surface  becomes  s emi -con i c a  1 . 
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•  Dissipation  is  smallest  at  the  plume  fringes  and  greatest 
on  the  plume  centerline.  The  off-centerline  spatial 
distribution  of  variance  dissipation  is  similar  to  the 
distribution  of  variance  itself.  In  fact,  Csanady  (1967) 
proposes  that  the  dissipation  of  concentration  fluctua¬ 
tions  is  proportional  to  the  fluctuation  intensity,  since 
the  dissipation  mechanism  should  not  depend  on  the 
amplitude  of  fluctuations.  The  expected  proportional 
relationship  is  verified  in  the  next  section  of  this 
chapter. 

By  far  the  most  interesting  behavior  exhibited  by  the 

measured  production  and  dissipation  profiles  in  Figs.  5.19 

through  5.23  is  that,  as  the  plume  moves  downwind,  production 

decreases  and  eventually  becomes  negligible  compared  with  the 

2~ 

dissipation  of  c'  .  Then,  net  dissipation  can  be  approximated 
by  the  dissipation  term  itself,  as  suggested  by  Fig.  5.23, 
where  the  triangles  ( <J>  )  and  the  dotted  curve  (<j>^)  are  seen  to 
lie  close  together.  The  downwind  decrease  of  production  is  a 
natural  consequence  of  plume  growth,  during  which  the  plume 
material  is  spread  more  and  more  thinly  over  the  plume  cross- 
section,  thereby  reducing  local  concentration  gradients.  This 
assumption  of  negligible  variance  production  will  play  a  major 

role  in  the  theoretical  derivation  of  the  dissipated  Gaussian 

2 

model  for  c‘  . 

One  last  point  should  be  made  regarding  the  production 
curves,  particularly  the  one  shown  in  Fig.  5.23.  Evaluation 


. 
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of  production  using  Eq.  (5.30)  requires  the  experimental  eddy 
diffusivity  calculated  from  Eq.  (5.19).  This  in  turn  requires 
the  x-jjradi ent  of  concentration,  which  is  used  in  the  advection 

3  C 

term  U—  .  This  derivative  was  calculated  from  constant-x  curves  of 
ax'  ^  through  x-deri vati ves  of  c(x)  curves  fitted  through  constant-z 

values  of  c(x,z)  measured  at  the  five  downwind  locations.  Because  the  data 
of  Fig.  5.23  was  obtained  at  the  last  downwind  location,  the  x  - 
derivative  there  is  an  extrapolated  value,  and  is  subject  to 
more  error  than  x-deri vati ves  at  the  interior  points.  This 
error  propagates  through  the  eddy  di f fu si vi ti es  to  appear  as 
error  in  the  production  curve.  The  error  bar  on  production 
(and  hence  on  net  dissipation)  shown  in  Fig.  5.23  corresponds 
to  an  error  of  ±  25 %  in  the  value  of  ||  at  this  last  downwind 
location. 

5.6  Concentration  Microscale  and  Dec ay -Time  Scale 

The  results  of  this  chapter  permit  the  testing  of 

two  assumptions  made  by  Csanady  (1967)  in  formulating  his 

model  for  concentration  fluctuations.  The  first  and  most 

important  assumption  was  that  the  dissipation  rate  <f>  of 

concentration  fluctuations  is  directly  proportional  to  the 

2 

local  concentration  variance  c'  .  The  second  assumption  was 
that  the  proportionality  factor  can  be  written  as  1 / T  ^ ,  where 
the  dissipation  time  scale  T^  increases  linearly  with  travel 
time  t  (i.e.  using  Taylor's  hypothesis,  a  linear  increase  with 
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distance  x  from  the  source).  These  assumptions  were  essential 
to  Csanady  '  s  solution  because  they  led  to  sel f -s imi 1 ar i ty  of 
variance  profiles,  and  permitted  the  unknown  dissipation 
function  <j>  to  be  expressed  in  terms  of  the  unknown  concen¬ 
tration  variance  S,  in  the  form 


t  +t  ,x 
o  d 


(5.31  ) 


Csanady  states  that  the  constants  t  and  t^  depend  on  the 

relative  'age1  difference  between  the  atmospheric  turbulence 

and  the  concentration  fluctuations.  In  decaying  grid  turbulence, 

vector  and  scalar  fluctuations  originate  at  the  same  point 

in  space,  and  therefore  have  the  same  'age';  for  this  special 

case,  Gibson  and  Schwarz  (1963)  have  shown  that  t  =0  and 

o 

Td  =  2/3 1 .  Csanady  supposes  that  for  an  individual  plume 
cross-section,  the  'age'  of  S  is  constant,  so  that  tQ  and  t  d 
are  independant  of  y  and  z. 

These  assumptions  are  tested  here  in  the  following  way: 
first,  equations  (5.29)  and  (5.31  )  are  combined  and  rearranged 
so  that  Td  is  expressed  in  terms  of  the  ratio  of  concentration 
variance  to  time  derivative  variance, 


Hinze  (1975)  has  applied  the  concept  of  a  Taylor  microscale  to 


' 


■ 


■ 


concentration  fluctuations,  and  shows  that  the  dissipation 
rate  may  be  written  as 
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(5.33) 


c 


Thus  the  Taylor  concentration  microscale  and  the  decay-time 
scale  are  related  by 


(5.34) 


The  time  derivative  variance  was  measured  in  the  wind 


tunnel  (see  Ch.  Ill  for  details)  along  with  S  and  U ,  and 
vertical  profiles  of  the  derived  values  of  concentration 
microscale  are  shown  in  Fig.  5.24.  There  is  some  small  variation 
with  height  above  surface,  but  it  is  not  systematic  and  is 
probably  due  to  experimental  error.  Scatter  in  the  lateral 
profiles  (not  shown)  was  about  the  same.  For  each  of  the  five 
plume  cross-sections  a  laterally  and  vertically  averaged  value 
of  T.  was  calculated.  The  results  are  presented  in  Table  5.2. 


The  values  of  Td  (in  ms)  have  a  best-fit  linear  dependance  on 
x  (in  cm)  of 


T,  =  38.5  +  1.65  x 


d 


(5.35 ) 


with  a  coefficient  of  determination  r 


2 


0,995. 
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TABLE  5.2 

VARIATION  OF  DECAY  TIME  SCALE  WITH  DOWNWIND  DISTANCE 


This  experimental  data  is  in  agreement  with  Csanady's  suppositions 
that : 

•  the  dissipation  rate  of  concentration  fluctuations  downwind 
of  an  elevated  point  source  in  an  atmospheric  boundary 
layer  is  proportional  to  the  local  concentration  variance; 

•  the  proportionality  factor  is  independant  of  y  and  z,  and 
varies  inversely  with  T  ^ ,  a  decay  time  scale  which 
increases  with  downwind  distance. 

The  values  of  t  and  t  .  must  depend  on  the  level  of 

0  ^ 

atmospheric  turbulence,  although  other  factors  may  also  be 

important;  the  question  of  how  to  determine  a  priori  the  values 

of  t  and  t  is  still  unanswered.  However,  for  the  simulated 
o  d 

plume  the  need  to  know  these  constants  is  not  important. 


\  1  bn  I 
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because  the  dissipated  Gaussian  model,  which  provides  a 
satisfactory  description  of  the  concentration  variance  field, 
is  shown  in  the  next  section  to  imply  a  specific  dependance 


of  Td  on  x. 

5 . 7  Theoretical  Derivation  of  the  Dissipated  Gaussian  Model 

2 

There  is  an  excellent  agreement  between  c'  calculated 

from  the  dissipated  Gaussian  model  proposed  in  Section  5.3.2 

2 

and  the  values  of  c1  measured  in  the  simulated  plume.  Further¬ 
more,  the  mathematical  form  of  the  model  is  very  similar  to 
the  classic  reflected  Gaussian  model  for  mean  concentration 
c".  This  suggests  that  the  diffusion  equation  (5.8)  for  c  and 
the  variance  balance  equation  (5.9)  may  have  somewhat  similar 
solutions. 

The  analysis  which  follows  uses  the  approach  of  Csanady 
(1967,  1973).  First,  the  diffusion  and  variance  equations  are 
rewritten  for  cylindrical  coordinates,  implying  a  plume 
symmetric  about  its  centerline,  far  above  the  ground.  A 
similarity  solution  to  the  equations  is  sought,  in  which  the 
wind  field  is  uniform  and  the  x-dependant  plume  growth  is 
incorporated  in  the  empirical  variation  of  the  plume  standard 
deviation  a.  Thus,  the  mean  concentration  is  assumed  to  have 

the  form 

c  =  cQ  f (n) »  (5.36) 

where  the  dimensionless  radial  coordinate  is 


n  =  r/a 


(5.37) 
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and  the  x-dependant  centerline  variation  of  c  is 


C 

o 


Q 


2tt U  o 
o 


oo 

/0nf  (  n  )  dn 


(5.38) 


For  radially  constant  K  =  K  =  K  = 

y  z 

equation  (5.8)  can  be  written  (with 


U  o 
o 


diffusion 


I  _ 


da  .  . 

di>  the  , 

,  fa  =  ^-4) 


if 

dn 


as 


dn 


f"  +  (-  +  n)  f1  +  2f  =  o, 
n 


(5.39) 


which  is  a  homogeneous  ordinary  linear  second  order  differen¬ 
tial  equation  with  variable  coefficients.  The  boundary 
conditions  of  axial  symmetry  and  decay  to  zero  at  large 
distances  are: 


f(°°)  =  0, 

f  '  (0)  =  0. 


(5.40) 


It  is  easy  to  show  by  substitution  that  one  solution  of  (5.39) 
satisfying  all  the  boundary  conditions  is  the  simple  Gaussian 
distribution 


2 

f (n )  =  exp  ( -  ) » 


(5.41) 


for  which  j^nf(n)dn  =  1.  This  gives  the  well  known  expression 
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for  mean  concentration  in  an  axisymmetric  plume. 


C  -  (-9 


7)  eXP  ( - 7) 

2ttU  _o  2 cr 


(5.42) 


0 


These  steps  are  now  repeated  for  the  variance  field, 
with  the  similarity  assumption 

S  =  sQCo2  g(n)  (5.43) 

and  the  dissipation  (J)  given  in  terms  of  S  by  use  of  Eq. 

(5.31).  Eq.  (5.43)  for  S  is  different  from  the  form  used  by 

Csanady  (1967,  1973)  only  in  the  use  of  a  source  variance 

2 

strength  constant  s  .  The  balance  equation  (5.9)  for  c' 
can  now  be  written  as 


-4g  -  ng'  =  g"  +  ^  9'  +  2f'2  -  ag. 

- - y - <  ' - » - ‘  ' - * - '  ' - * - ’ 

advection  gradient  flux  production  dissipation 

Introducing  the  solution  (5.41)  for  f  and  collecting  terms 
yields  the  form  given  by  Csanady  (1967,  1973), 

O  9 

q"  +  (-  +  q)g'  +  ( 4-a) g  =  -2n  exp(-n  ),  (5.44) 

3  n 

which  i s  a  nonhomogeneous  ordinary  linear  second  order  differ¬ 
ential  equation  with  variable  coefficients.  The  parameter  a. 
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given  by 


2 

a  =  KT  *  (5.45) 

d 

results  from  dissipation  of  c  '  2 ,  and  is  assumed  constant 
throughout  the  plume.  The  nonhomogeneous  term  on  the  right 
hand  side  of  (5.44)  is  due  to  production  of  c'2. 

The  boundary  conditions  on  (5.44)  are,  in  this  case, 
the  same  as  for  the  mean  field:  a  decay  to  zero  at  large 
distances,  and  an  axisymmetric  plume, 


g(°°)  =  0, 

9 ' ( 0)  =  0. 


(5.46) 


Csanady  (1967)  has  obtained  a  closed  form  solution  to  (5.44) 
subject  to  boundary  conditions  (5.46).  What  is  sought  here  is 
an  approximate  solution  to  (5.44)  that  can  be  generalized  to 
the  more  complicated  case  of  an  asymmetric  plume  near  a  solid 
boundary . 

A  considerable  simplification  of  (5.44)  results  if 
the  nonhomogeneous  production  term  is  omitted.  The  figures  of 
Section  5.5  show  that  production  does  in  fact  become  negligible 
as  the  simulated  plume  moves  downwind.  For  the  assumption 
of  zero  production,  (5.44)  becomes  homogeneous. 
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9"  +  (^  +  n ) g  '  +  (4-a)g  =  0  (5.47) 

which  is  substantially  closer  to  the  form  of  the  mean  equation 
(5.39).  Only  one  small  step  is  needed  to  make  the  two 
equations  identical,  and  that  is  to  specify 

a  =  2. 

Then 

9"  +  +  n)g‘  +  2g  =  0,  (5.48) 

one  solution  of  which  is 

2  2 

g(n)  =  exp ( -  £-)  =  exp ( -  ?—?).  (5.49) 

2a^ 

In  cartesian  coordinates,  the  variance  becomes 

2  2  2 

S  =  s n ^ exP  (-^j)  exp  ,  (5.50) 

0  2ttU  2 cr  2a 

o 

which  is  very  close  in  form  to  the  dissipated  Gaussian  model 
proposed  for  S  earlier  in  this  chapter.  As  before,  sq  is  a 
source  variance  strength  parameter  to  be  determined  from 
exper i ment . 

By  analogy  with  the  known  mean  field  solution  for 


an  asymmetric  plume,  Eq.  (5.50)  can  be  generalized  for 

a  f  o  , 
y  z 
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2 


2 


2 


)  exp  (~y  o)  exp  ( 
2a 

y 


-z 


z 


2 '  ’ 


or,  shifting  the  z-coordinate  by  a  constant  amount,  so  the 
plume  centerline  is  at  z  =  H, 


2  2  2 
)  exp  exp  ].  (5.51) 

2oy  2oz 


S  =  s  ( - 3 - 

o  '  2ttU  a  a 
o  y  z 


Note  that  the  linearity  of  Eq.  (5.48)  permits  the 


solution  of  a  complicated  problem  to  be  written  as  a  super¬ 
position  of  solutions  to  less  difficult  problems.  Eq.  (5.48) 
remains  linear  when  written  in  cartesian  coordinates,  so  that 
the  process  of  superposition  lends  itself  to  the  difficult 
problem  of  solving  (5.48)  (in  cartesian  coordinates)  subject 
to  the  boundary  condition  of  complete  dissipation  of  fluctua¬ 
tions  at  a  solid  boundary  located  at  z  =  0,  a  distance  H  below 
the  plume  centerline: 


g ( o )  =  o. 


A  solution  to  this  problem  is  obtained  by  introducing  a 
negative  image  source  a  distance  H  on  the  opposite  side  of  the 
solid  boundary.  This  procedure  is  the  same  as  that  used  to 
derive  the  reflected  Gaussian  model  for  mean  concentration. 

The  S-field  is  then  given  by  a  superposition  of  the 
fields  from  both  the  positive  (real)  and  negative  (image) 


sources  , 
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(5.52) 


The  dissipated  Gaussian  model  ,  originally  proposed  for  S  on  the 
basis  of  empirical  arguments,  has  now  been  derived  as  a 
solution  of  the  variance  balance  equation  (5.9)  for  the 
special  case  of 

•  Negligible  variance  production, 

•  Dissipation  parameter  a  =  2. 

•  Constant  diffusivities  K  =  K 

y  y 

da 

1/1  —  1/  —  I  I  rr  Z 


The  physical  significance  of  setting  the  dissipation 
parameter  a  =  2  must  be  determined.  As  was  mentioned  in 
Section  5.6,  Gibson  and  Schwarz  (1963)  have  shown  that  the 
decay  time  scale  of  fluctuations  in  isotropic  grid  turbulence 
is  related  to  distance  from  the  grid  by  the  expression 

Td=§.  (5.53) 

2 

Csanady  (1973)  has  related  the  dissipation  parameter  a  =  a  /KTd 
to  downwind  distance  by  considering  the  asymptotic  stage  of 
diffusion  in  which  a  ^  ^  x,  and  obtains 


da 


=  U  a 


o  y  dx 


and 


■  f 
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2 ( x-x  ) 
_ _ <r 

aU 


(5.54) 


Obviously,  a  =  3  corresponds  to  dissipation  of  scalar  fluctu¬ 
ations  in  decaying  isotropic  grid  turbulence  described  by 
Eq.  (5.53).  Values  of  a  <  3  can  be  expected  in  turbulent 
flows  which  decay  more  slowly  than  grid-type  turbulence. 

Thus,  a  =  2  means  that  the  decay  time  scale  is  exactly  equal 
to  the  travel  time  of  fluctuations  from  their  effective  origin 
X  =  x0. 

Eq.  (5.54)  can  be  written  as 


T 


d 


a  +  bx  , 
s 


(5.55) 


where  x  is  distance  from  the  plume  source.  Combining  (5.54) 
s 

and  (5.55)  , 


from  which 


(5.56) 


(5.57) 


relative  to  the  plume  source.  Comparing  (5.55)  and  (5.57), 
the  dissipation  parameter  a  is  given  by 


a 


2_ 
Ub  * 


(5.58) 


A. 


■ 
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Measured  values  of  Td  (in  ms)  obtained  during  the 
present  study  fit  the  linear  relation 

Td  =  38.5  +  1 . 65  xs  .  (5.59) 

Taking  a  typical  velocity  in  the  tunnel  as  U  =  800  cm/sec, 
(5.57)  and  (5.59)  may  be  combined  to  give 

(5.60) 

This  shows  that  helium  fluctuations  in  the  wind  tunnel  were 
dissipated  at  only  half  the  rate  that  would  occur  in  isotropic 
grid  turbulence,  for  which  a  =  3 .  More  importantly,  the  wind 
tunnel  tracer  fluctuations  died  away  25%  slower  than  the 
Gaussian  solution  (a  =  2)  requires.  However,  the  favorable 
comparision  of  measured  S-profiles  (Figs.  5.1  through  5.7) 
and  measured  profiles  of  the  fluctuation  intensity  i  =  ■Js'/'c 
(Figs.  5.8  through  5.13)  with  the  Gaussian  model  prediction 
shows  not  much  sensitivity  to  variations  in  a.  The  profile 
shapes  are  accurately  predicted,  and  the  axial  variation  (Fig. 
5.1)  is  also  well  reproduced.  This  remains  true  even  close  to 
the  source,  where  the  assumption  of  negligible  production  is 
certainly  violated  (Fig.  5.19).  The  conclusion  drawn  is  that 
the  dissipated  Gaussian  model,  a  true  solution  of  the  variance 
equation  only  for  the  special  case  of  zero  production  and  decay 
time  scale  T^  equal  to  effective  travel  time,  is  in  general  a 


- 
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good  engineering  approximation  to  the  spatial  distribution  of 
concentration  variance  in  plumes. 

5 . 8  Summary 

The  most  important  results  from  the  measurements  and 
analysis  presented  in  this  chapter  are  summarized  here: 

•  For  steady  flow  situations  where  the  parameters  Q, 

UH*  ay  and  az  are  accurately  known,  the  reflected 
Gaussian  model  (5.12)  for  mean  concentration  fits 
observations  with  relative  error  less  than  10%  of  the 
local  concentration  except  at  the  plume  edges  and  near 
the  ground,  where  errors  of  20-40%  can  be  observed; 

•  In  flows  for  which  the  reflected  Gaussian  model  for 

mean  concentration  c  can  be  used,  the  concentration 
2 

variance  c‘  can  be  calculated  from  the  dissipated 
Gaussian  model  given  by  Eq.  (5.13).  This  model  in¬ 
cludes  a  source  variance  term  sQ  which  is  constant 
through  the  plume,  and  which  probably  depends  on 
atmospheric  turbulence  level  at  source  height; 

2 

•  The  local  production  and  dissipation  rates  of  c'  in 
a  plume  are  not  equal.  In  fact,  as  the  plume  moves 
downwind,  production  becomes  negligible  in  comparison 
with  the  variance  dissipation.  This  leads  to  a  balance 
equation  for  c 1 ^  which  has  the  same  form  as  the  diffusion 


equation  for  c; 
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O 

•  The  dissipated  Gaussian  model  for  c‘  is  a  true  solution 
of  the  variance  equation  for  the  special  case  of  zero 
in-plume  variance  production  and  fluctuations  which 
decay  at  two-thirds  the  rate  of  those  in  isotropic 
grid  turbulence; 

•  Measurements  of  the  local  dissipation  rate  of  concen¬ 
tration  fluctuations  support  Csanady's  (1967)  supposi¬ 
tion  that  this  rate  is  proportional  to  the  local  concen¬ 
tration  variance; 

•  The  measurements  confirm  Csanady's  (1967)  assumption 
that  the  reciprocal  of  this  proportionality  factor  is  a 
decay  time  scale  T^  which  is  constant  at  a  given  section 
and  which  increases  linearly  with  distance  from  the 
source. 
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FIG.  5.1  Centerline  Variation  of  Mean  and 
RMS  Concentration 


Centerline  Concentration  Fluctuation  Intensity 


z/H  Dimensionless  Height  Above  Surface 
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FIG.  5.3  Vertical  Variation  of  Mean  and  RMS 

Concentration  at  x/H=4.3 


z/H  Dimensionless  Height  Above  Surface 
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Normalized  Mean  and  RMS  Fluctuating  Concentration 


FIG.  5.4  Vertical  Variation  of  Mean  and  RMS 

Concentration  at  x/H=  6.3 


z/H  Dimensionless  Height  Above  Surface 
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Normalized  Mean  and  RMS  Fluctuating  Concentration 


FIG. 5.5  Vertical  Variation  of  Mean  and  RMS 

Concentration  at  x/H=9.2 


z/H  Dimensionless  Height  Above  Surface 
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Normalized  Mean  and  RMS  Fluctuating  Concentration 


FIG.  5.6  Vertical  Variation  of  Mean  and  RMS 

Concentration  at  x/H=I2.5 


z/H  Dimensionless  Height  Above  Surface 
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Normalized  Mean  and  RMS  Fluctuating  Concentration 


FIG.  5.7  Vertical  Variation  of  Mean  and  RMS 

Concentration  at  x/H  =  I9.2 
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FIG.  5.8  Lateral  Distribution  of  Concentration 

Fluctuation  Intensity  at  Plume  Height 


Downwind  Distance 
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FIG.  5.9  Vertical  Distribution  of  Concentration 

Fluctuation  Intensity  at  x/H=4.3 


Downwind  Distance 
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FIG.  5.10  Vertical  Distribution  of  Concentration 

Fluctuation  Intensity  at  x/H=6.3 


\ 


z/H  Dimensionless  Height  Above  Surface 
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Normalized  Concentration  Fluctuation  Intensity 
■cH 


FIG.  5.11  Vertical  Distribution  of  Concentration 
Fluctuation  Intensity  at  x/H=9.2 


rv 


z/H  Dimensionless  Height  Above  Surface 
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~r^—  Normalized  Concentration  Fluctuation  Intensity 
»cH 

FIG.  5.12  Vertical  Distribution  of  Concentration 

Fluctuation  Intensity  at  x/H=I2.5 


z/H  Dimensionless  Height  Above  Surface 
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7-— —  Normalized  Concentration  Fluctuation  Intensity 
'cH 


FIG.  5.13  Vertical  Distribution  of  Concentration 

Fluctuation  Intensity  at  x/H  =  I9.2 
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z/H  Dimensionless  Height  Above  Surface 
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FIG.  5.19  Vertical  Profile  of  Variance  Production  and 
Dissipation  on  Plume  Centerline  at  x/Hs4.3 


yi! 


z/H  Dimensionless  Height  Above  Surface 
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FIG.  5. 20  Vertical  Profile  of  Variance  Production  and 
Dissipation  on  Plume  Centerline  at  x/H=6.3 


z/H  Dimensionless  Height  Above  Surface 
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FIG.  5.21  Vertical  Profile  of  Variance  Production  and 

Dissipation  on  Plume  Centerline  at  x/H  =  9.2 


z/H  Dimensionless  Height  Above  Surface 
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FIG.  5.2  2  Vertical  Profile  of  Variance  Production  and 

Dissipation  on  Plume  Centerline  at  x/H=I2.5 


Production  A  —  Dissipation 
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FIG.  5.24  Vertical  Profile  of  Concentration  Microscale 
at  Several  Downwind  Locations 


CHAPTER  VI 


STATISTICAL  PROPERTIES  OF  THE 
CONCENTRATION  FLUCTUATIONS 


6.0  Introduction 

Most  practical  applications  of  atmospheric  dispersion 
theory  focus  on  identifying  the  mean  concentration  expected 
at  a  given  point  after  a  specified  averaging  time.  However, 
it  is  well  known  that  the  concentration  at  a  receptor  varies 
with  time,  and  a  concentration  history  which  results  in 
moderate  mean  concentrations  can  include  short-term  peak  values 
which  are  much  higher  than  the  mean.  A  measure  of  this 
behavior  is  the  'peak-to-mean'  concentration  ratio,  where  the 
'peak'  concentration  is  defined  as  that  value  which  is  not 
exceeded  more  than  a  fraction  P  of  the  time.  Knowledge  of  the 
peak  concentration  at  a  receptor  is  necessary  when  dealing 
with  pollutants  which  can  have  damaging  effects  at  high  concen¬ 
trations  over  short  times.  Examples  are  the  poisonous  gases 
hydrogen  sulphide  and  carbon  monoxide,  the  explosive  gases 
methane  and  hydrogen,  and  gases  such  as  sulphur  dioxide  which 
are  often  vented  to  atmosphere  and  to  which  some  vegetation  is 
sensitive. 

Calculation  of  peak  concentrations  requires  that  the 
probability  of  occurrence  be  known  for  the  entire  range  of 
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concentrations.  The  most  compact  way  to  present  such  data  is 
with  the  probability  density  function,  or  pdf.  In  this  chapter, 
an  intermittent  lognormal  pdf  is  proposed  as  a  model  for  con¬ 
centration  fluctuations  downwind  of  an  elevated  point  source 
in  the  atmospheric  boundary  layer.  The  model  input  parameters 
are  mean  concentration  c,  concentration  variance  c'2  and  inter- 
mittency  y,  and  can  be  computed  from  mathematical  models  or 
measured  in  the  field  or  laboratory  without  the  need  for  curve 
fitting.  The  pdf  model  has  been  compared  with  wind  tunnel 
measurements,  and  predicts  peak-to-mean  concentration  ratios 

within  20%  of  observed  values  so  long  as  the  dilution  parameter 
.  2  . 

D / i c  is  greater  than  400.  Peak  concentrations  at  the  ground 
occur  about  two  thirds  of  the  distance  out  to  the  point  of 
maximum  mean  concentration. 

6  .  1  Eddy  Dilution 

In  experimental  studies  of  atmospheric  diffusion  one 
often  observes  periods  of  zero  concentration  directly  downwind 
of  the  source  of  emissions.  This  means  that  the  sampled  eddy 
contains  no  portion  of  the  diffusing  cloud.  At  the  other  extreme 
is  observation  of  an  eddy  which  consists  of  pure  contaminant. 

In  most  real  situations  the  atmosphere  provides  dilutant  eddies 
in  such  numbers  that  the  contaminant  eddies  are  soon  attenuated 
by  the  dilution  process.  In  a  dispersing  plume  we  find  pure 
contaminant  eddies  very  rarely,  pure  dilutant  eddies  more 
frequently  and,  most  frequent  of  all,  the  intermediate  eddies 


. 


whose  concentration  fluctuations  are  the  subject  of  this 
chapter . 

Pure  contaminant  eddies  are  restricted  mainly  to  a 
core  region  at  the  stack  exit,  similar  in  shape  to  the  non- 
turbulent  core  region  of  a  laminar  jet  in  a  turbulent  cross- 
flow.  These  eddies  are  diluted  to  the  i n termed i ate  form 
so  quickly  that,  except  for  points  a  few  stack 
diameters  from  the  source,  they  can  be  ignored  as  a  class. 

Of  the  remaining  two  eddy  types,  the  frequency  of  occurrence 
of  the  pure  dilutant  eddies  can  be  described  with  a  single 
parameter  y,  the  concentration  i nt erm i ttency  factor  discussed 
in  Chapter  IV.  The  non-zero  fraction  of  the  concentration 
readings  (from  the  intermediate  eddies)  have  concentrations 
of  0  <  c  <  c  .  The  probability  of  observing  a  concentra- 

tion  which  lies  in  the  range  c  to  c+Ac  is  characterized  by 
p(c),  the  probability  density  function,  or  pdf,  of  eddy 
concentration.  The  probability  density  function  Pj(c)  for  the 
intermediate  eddies  has  the  property  that,  if  it  is  integrated 
over  all  concentrations, it  must  give  a  probability  of  y.  The 
delta-function  pdf  pD ( c )  for  the  pure  dilutant  eddies  must 
integrate  to  1 -y ,  so  that  the  total  pdf  p(c)  will  integrate  to 
unity  as  it  should. 

The  detailed  form  of  the  pdf  for  concentration  fluctu¬ 
ations  is  the  subject  of  much  investigation.  Csanady  (1973,  p. 
225)  has  presented  a  simple  physical  model  to  justify  a  lognormal 
distribution  for  the  pdf  of  non-zero  concentrations.  His 
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heuristic  argument  begins  with  the  assumption  that  eddy  con¬ 
centrations  are  reduced  by  a  random  series  of  independant 
diluting  steps,  so  that  after  N  dilutions  an  eddy's  concentra¬ 
tion  can  be  written  as  csource  multiplied  by  N  positive  random 
numbers  less  than  unity.  The  logarithms  of  a  multiplicative 
series  are  additive,  so  the  Central  Limit  Theorem  leads  to  a 
normal  distribution  for  the  logarithms  of  eddy  concentration 
if  the  number  N  of  diluting  steps  is  large.  Thus  the  concen¬ 
trations  themselves  are  not  normally  distributed,  but  have  what 
is  called  a  'lognormal'  distribution  given  by 


p(c) 


fZ ¥  a^c 


-In  (c/c  ) 

exp[ - 2-2-]  , 

2a, 


(6.1) 


where  cQ  is  the  median  concentration  of  the  lognormal  distribu¬ 
tion  and  a ^  is  the  logarithmic  standard  deviation.  Csanady 

(1973,  p.  229)  shows  that  these  two  parameters  are  related  to 

—  .  2 

concentration  mean  c  and  variance  c'  by 

°L  =  /ln(l  +  1CZ)  (6.2) 


and 


(6.3) 


where  fluctuation  intensity  1£  has  been  previously  defined  as 
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ic  =  c.  This  two-parameter  model  requires  the  mean 

—  '  P 

concentration  c  and  concentration  variance  c  1  as  input,  and 
does  not  account  for  the  intermittency  of  concentration  fluc¬ 
tuations.  The  two-parameter  lognormal  pdf  model  is  used  later 
in  this  chapter  as  the  starting  point  for  a  more  general  three- 
parameter  pdf  which  includes  intermittency. 

A  two-parameter  model  which  does  incorporate  inter¬ 
mittency  is  the  exponential  distribution  discussed  by  Barry 
(1974).  The  pdf  parameters  are  mean  concentration  c  and 
intermittency  y,  and  appear  in  the  cumulative  distribution 
function 


(P(c) 


fc 

p(c)dc  =  1  -  yexp[— . 


Barry  assumed  intermittency  to  result  from  two  effects,  so  that 

Y  =  Y0-  Yturb  , 

where  y  is  due  to  large  scale  fluctuations  in  horizontal  wind 
direction  and  Yturb  is  due  t0  atmospheric  turbulence.  For  wind 
tunnel  experiments  the  intermittency  must  be  generated  solely 
by  turbulence,  i.e.  y  =  Yturb>  since  y6  -  1  results  from  the 
physical  constraint  of  the  tunnel  walls.  However,  actual  field 
measurements  of  y  are  n°t  so  constrained.  The  result  is  that 
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short  data  averaging  times  (t  ~  hours)  produce  y-y.  ,  .and 

a  v  t  u  r  b 5 

longer  averaging  times  (tav>  days),  Y  -  Y0 •  Barry  obtained  the 
best  fits  to  P(c)  data  from  multiple  receptor  sites  with  long 
averaging  times,  when  the  model  assumption  that  individual 
concentration  observations  were  independant  was  most  likely 
to  be  valid. 

The  Weibull  three-parameter  distribution  has  been  used 
by  Apt  (1976)  to  fit  two-week  averages  of  atmospheric  radio¬ 
activity  data  taken  over  periods  of  one  year.  The  model  assumes 
that  y  =  1 »  and  has  the  form 

P(c)  =  1  -  exp[-(c-6)3/a]  , 

where  a,  3  and  6  are  empirically  adjusted  parameters.  The 
parameter  6  can  be  interpreted  either  as  the  background  concen¬ 
tration  or  the  lowest  observed  concentration,  in  which  case  its 
value  can  be  obtained  without  curvefitting.  However,  if  6  is 
used  as  a  fitting  parameter  (as  done  by  Apt)  then  it  must  be 
interpreted  as  the  lowest  concentration  likely  to  be  observed. 
The  shape  parameter  3  is  entirely  adjustable.  If  3=1  then  c 
has  an  exponential  distribution  as  described  by  Barry  (1974). 

If  3=2,  the  Rayleigh  distribution  results.  For  3-3.25  the 
concentration  distribution  would  be  close  to  Gaussian.  Usually 
the  scale  parameter  a  is  also  adjustable.  However,  for  the 
exponential  case  of  3=1  and  <5  =  0,  the  parameter  would  be  a=c. 
Such  special  cases  provide  the  only  opportunity  to  calculate 


. 
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a  without  curvefitting. 

A  'censored'  three-parameter  lognormal  distribution 

has  been  proposed  by  Ott  and  Mage  (1976)  as  a  model  which  can 

incorporate  a  specified  probability  of  zero  concentrations. 

Although  the  authors  do  not  say  so,  the  model  input  parameters 

can  be  cast  in  the  form  of  the  mean  concentration  c,  the  con- 

2 

centration  variance  c‘  and  the  i ntermi ttency  factor  y.  The 
calculation  technique  would  be  first  to  use  the  observed 

_  2 

values  of  c  and  c'  to  compute  the  two-parameter  lognormal 
pdf  (for  which  y=l).  Then  the  pdf  curve  is  displaced  to  the 
left  (toward  'negative'  concentrations)  until  the  area  under 
the  portion  to  the  left  of  the  origin  is  equal  to  1-y.  This 
probability  mass  is  then  lumped  as  a  delta  function  at  the 
origin,  so  that  negative  concentrations  do  not  occur.  As  a 
result  of  this  procedure,  the  most  probable  concentration, 
which  occurs  at  the  highest  point  on  the  pdf  curve,  is  reduced 
by  i ntermi ttency ,  whereas  the  probability  of  its  occurrence 
is  unchanged.  This  is  not  realistic,  because  the  more  inter- 
mittency  is  observed,  the  less  probable  is  observation  of  a 
specified  non-zero  concentration.  Also,  because  the  zero  and 
non-zero  populations  are  essentially  independant,  the  magnitude 
of  the  most  probable  non-zero  concentration  should  be  unaffected 
by  the  presence  of  i ntermi ttency.  Mage  and  Ott  did  not  use  the 
observed  mean  and  variance  as  input  parameters,  but  iteratively 
varied  their  parameters  (related  to  c  and  c  )  so  as  to  minimize 
the  least  squares  error  when  fitting  the  observed  cdf. 


. 
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With  very  few  exceptions  the  accepted  procedure  for 
verifying  any  model  for  the  pdf  of  concentration  fluctuations 
has  been  to  adjust  the  model  input  parameters  to  obtain  the 
least  squares  error  fit  to  the  observed  cumulative  distribu¬ 
tion  function.  This  method  has  several  disadvantages.  First, 
fitting  the  observed  cdf  is  not  as  sensitive  an  indicator  of 
an  incorrect  pdf  model  as  is  matching  the  observed  pdf  itself. 
Secondly,  the  curvefitting  approach  hides  errors  in  the  pdf 
model  that  might  show  up  if  the  model  parameters  were  calcu¬ 
lated  directly  from  observable  variables  like  c,  c'2  and  y. 
Thirdly,  curvefitting  is  a  cumbersome  way  to  determine  the  pdf 
parameters  in  situations  where  real-time  data  is  being  pro¬ 
cessed,  such  as  the  real-time  air  quality  prediction  systems 
used  by  some  industries  for  supplementary  emission  control. 
Whenever  the  atmospheric  or  source  conditions  change,  the 
fitting  process  must  be  repeated  with  a  new  observed  time  series 
which  can  not  begin  until  steady  conditions  are  established. 

The  purpose  of  this  chapter  is  to  develop  and  verify  a 
probability  density  model  for  concentration  fluctuations  which 
incorporates  concentration  i n termi ttency  in  a  physically  consis¬ 
tent  manner  and  whose  parameters  can  be  computed  from  models 
or  measured  in  the  field  or  laboratory  on  a  real-time  basis. 

6. 2  The  Intermittent  Lognormal  Distribution 

The  basic  two-parameter  lognormal  distribution  develop¬ 
ed  by  Csanady  (1973)  is  attractive  as  a  model  for  concentration 


. 
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fluctuations  because: 

•  the  model  is  based  on  a  plausible  heuristic  argument; 

•  the  parameters  required  are  concentration  mean  c  and 


2 


and  can  be  obtained  from  mathematical 


variance  c 


models  or  from  measurements  in  the  field  or  laboratory. 

The  disadvantage  of  the  two-parameter  lognormal  model  is  that 
there  is  no  provision  for  concentration  i ntermi ttency  (a  non¬ 
zero  probability  of  zero  concentrations). 

If  the  eddy  population  consists  of  two  groups, 

1.  contaminated  eddies  which  are  diluted  in  a  lognormal 
manner  as  described  by  Csanady  (1973,  Ch.  7),  and 

2.  uncontaminated  eddies  whose  observed  effect  is  signal 
i n term i ttency  , 

then  the  probability  of  observing  an  uncontaminated  eddy  is 
1-Y,  and  the  probability  of  observing  a  contaminated  eddy  is 
simply  y.  It  follows  that  the  probability  of  observing  an 
eddy  with  a  specified  degree  of  contamination  must  be  y  times 
the  lognormal  probability  function  proposed  by  Csanady  for 
non-zero  concentrations.  Then  the  probability  density  function 
p(c)  for  intermittent  concentration  fluctuations  must  have  the 


form 


p(c)  =  (l-y)6(c) 


+ 


9rt*  * 
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where  6(c)  is  a  delta  function  located  at  the  origin  c  =  0. 

The  cumulative  probability  distribution  is  obtained  by  integrating 


(6.4), 


P(c)  = 


(6.5) 


0 


which  gives  [P(°°)  =  1  and  (P(0)  =  1 -y  as  required.  An  express¬ 
ion  similar  to  Eq.  (6.5)  has  also  been  derived  by  Aitchison 
and  Brown  (  1  957  ,  p.  95)  in  the  context  of  economics  data  analysis  to 
account  for  zero  observations  in  an  otherwise  lognormal  population. 

In  order  to  plot  equation  (6.4)  for  comparison  with 
experiment,  all  the  parameters  in  p(c)  must  be  known.  These 
parameters  are  the  concentration  i ntermi ttency  y,  and  the 
median  concentration  cQ  and  logarithmic  standard  deviation 
of  the  concentration  fluctuations.  The  usual  method  would  be 
to  adjust  y,  cQ  and  a L  until  Eq.  (6.4)  for  p(c)  matches  the 
experimental  data.  However,  the  three  adjustable  parameters 
provide  so  much  flexibility  to  the  shape  of  p(c)  that  simply 
being  able  to  fit  Eq.  (6.4)  to  data  does  not  establish  the 
validity  of  the  intermittent  lognormal  assumption.  To  test 
the  lognormal  hypothesis  it  is  essential  to  compute  these 
three  parameters  from  directly  observable  variables.  The 


concentration  mean  and  variance  are  related  to  and  cQ  by 
taking  moments  of  the  pdf  about  the  origin  c  -  0.  Thus 
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f  00 

c  p(c  )dc  =  YC0exp(-a.2/2)  (6.6) 

0 


and 

[CO 

c'2  =  (c-c)2  p ( c  )  d  c  =  YC02exp(aL2)[exp(a  2)-y]  , 

j  0 

(6.7) 

which  can  be  combined  to  give  the  logarithmic  standard  deviation 

°L  =  /ln[y  (1  +  ic2)]  (6.8) 

and  the  median  concentration 

co  -  ~  C  (6.9) 

Y  -j  Y  ( 1  +  1  £  ) 

of  the  1 ognormal ly-distri buted  concentrations. 

If  the  i n termi ttency  y  is  set  to  unity,  these  expressions 

for  a.  and  c  reduce  to  the  non-i ntermi ttent  Eqs.  (6.2)  and 
L  o 

(6.3)  derived  by  Csanady  (1973).  It  is  also  possible  to  sub¬ 
stitute  into  Eqs.  (6.6)  and  (6.7)  the  expressions  (4.23)  and 
(4.24)  developed  by  Becker,  Hottel  and  Williams  (1965)  for  the 
mean  and  variance  of  non-zero  concentrations  in  an  intermittent 
concentration  field.  The  result  is  once  again  a  reduction  to 
the  non-i ntermi ttent  Eqs.  (6.2)  and  (6.3)  for  and  cQ.  This 
is  particularly  interesting  because  it  means  that,  in  spite  of 
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the  explicit  appearance  of  i ntermi ttency  in  Eqs.  (6.8)  and 

(6.9) ,  the  values  of  a L  and  cq  are  in  fact  independant  of  y. 

Thus  the  explicit  i ntermi ttency  in  Eqs.  (6.8)  for  aL  and  (6.9) 
for  cq  is  exactly  cancelled  by  the  implicit  effect  i  ntermi  ttency 
has  on  the  values  of  c  and  c'2.  It  follows  that  if  only  the 
non-zero  concentrations  were  used  to  calculate  c  and  c'2,  a 
and  cq  could  be  computed  from  the  non- i nterm i tte nt  expressions 
(6.2)  and  (6.3).  However,  it  is  more  common  that  c  and  c'2 
are  formed  over  periods  which  include  instances  of  zero  concen¬ 
tration,  in  which  case  the  intermittent  expressions  (6.8)  and 

(6.9)  must  be  used  for  o.  and  c  . 

L  o 

By  directly  computing  the  parameters  cQ,  a ^  and  y,  we 
avoid  the  uncertainties  which  arise  from  curvefitting.  If  the 
pdf  of  equation  (6.4)  is  computed  using  experimental  measure¬ 
ments  and  is  found  to  be  a  good  fit  to  the  observed  pdf,  the 
intermittent  lognormality  of  concentration  fluctuations  would 
appear  to  be  verified.  Another  advantage  of  physically  meaning¬ 
ful  pdf  parameters  is  that,  for  simple  geometries  and  boundary 

_  2 

conditions,  theories  exist  which  describe  the  c,  c'  and  y 
fields  in  plumes,  and  can  provide  all  input  data  necessary  for 
a  complete  probabilistic  description  of  the  concentration  field. 
Chapter  V  discusses  the  use  of  a  diffusion  equation  solution 
for  the  c*2  field,  and  presents  a  new  theoretical  model  for  the 
transport  and  diffusion  of  c'2  which  is  similar  to  the  classic 
Gaussian  plume  model  for  mean  concentrations. 


_ 
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6.2.1  The  Effect  of  Gaussian  Noise  on  Measured  PDFs 

Concentration  observations  in  the  field  or  laboratory 
are  subject  to  measurement  errors  which  obscure  some  of  the 
information  contained  in  the  observations.  Often  the  concen¬ 
tration  signal  is  mixed  with  noise  from  the  signal  processing 
electronics,  or  the  sampling  process  is  subject  to  small 
randomly  distributed  errors.  For  the  present  study  the  noise 
level  in  the  fast  response  detection  system  was  measured  by 
simply  turning  off  the  helium  source.  In  the  field  the  back¬ 
ground  noise  can  be  measured  whenever  the  wind  is  not  blowing 
contaminants  from  the  source  to  the  receptor.  If  the  data  and 
noise  are  statistically  independant,  signal  subtraction  allows 
easy  measurement  of  the  concentration  mean  and  variance.  How¬ 
ever,  a  parameter  such  as  the  probability  density  of  the  con¬ 
centration  fluctuations  will  be  inseparably  mixed  with  the 
noise,  and  signal  processing  can  not  extract  it.  The  observed 
signal  pdf  is  a  composite  of  the  unknown  concentration  pdf 
and  the  noise  pdf.  Thus  the  observed  signal  is  given  by 

s=c+n,  (6.10) 

where  s  =  s  +  s1  is  the  instantaneous  observed  signal, 
c  =  c  +  c‘  is  the  instantaneous  concentration  and  n  =  n  +  n' 
is  the  noise  component.  If  n  is  independant  of  c,  then 


■  . 
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c  =  s  “  n  (6.11) 

and 

c‘2  -  s*2  -  (6.12) 

_  2 

where  s  and  s  are  measured  when  contaminant  blows  over  the 
receptor,  and  n  and  n‘  when  it  does  not.  In  this  way  the  pdf 
of  n  can  also  be  directly  observed. 

If  p  -|  ( n )  is  the  known  noise  pdf  and  p2(c)  is  the 
unknown  concentration  pdf,  then  the  observed  signal  pdf  p(s) 
is  given  by  the  convolution  integral  (Spiegel  (1975,  p.  47)), 


,00 


p(s)  = 


r  oo 


P1  (n)p2(s-n)dn  = 


P  2  ( c  )  P i (s-c)dc  .  (6.13) 


_  oo 


0 


In  the  experiments,  the  noise  pdf  observed  with  the  helium 
source  turned  off  appeared  to  have  a  normal  distribution,  with 
mean  n  corresponding  to  zero  ppth  helium  and  standard  deviation 
a  equivalent  to  0.2  ppth  helium.  Using  a  Gaussian  pdf  for  the 
noise,  and  assuming  the  intermittent  lognormal  model  for  con¬ 
centration  fluctuations,  Eq.  (6.13)  becomes 


p(s)  = 


oo 


0 


— ! -  exp['(s~C2nj-]U1-T)g(c)  + 

°n/TiF  2an 

-ln2(c/c  ) 

■  Y - -exp[ - 2 - * 


(6.14) 


/2  7T  CJ ,  C 


2a, 
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which  is  partly  integrated  to  yield 


p(s)  = 


Y-  e x p [ -^s ~ ^  ]  + 


°n/  2,7 


Y 


2  a 


n 


2ira  a, 

n  L 


oo 


1  fs  c  ^ n  (°/ CQ) 

lexP[-il---c:n) - °]dc 

n  2a  c  2  a, 

0  n  L 


(6.15) 


The  integral  remaining  in  this  expression  for  p(s)  must  be 
evaluated  numerically. 

The  presence  of  noise  in  the  concentration  signal 
means  that  the  intermittent  lognormal  model  (6.4)  cannot  be 
tested  directly.  However,  Eq.  (6.15)  for  the  observed  signal 
is  based  on  the  intermittent  lognormal  assumption;  if  an 
observed  concentration  pdf  known  to  contain  Gaussian  noise  is 
well  predicted  by  this  theoretical  pdf,  then  the  assumption  of 
intermittent  lognormality  for  concentration  fluctuations  will 
be  validated. 


6 . 3  Comparison  of  Theory  with  Experiment 

A  number  of  pdfs  were  observed  in  the  wind  tunnel  at 

_  2 

the  same  time  that  measurements  were  made  for  n,  a^,  c,  c1  and 
y.  The  data  was  used  to  calculate  the  theoretical  pdf  p(s) 
from  the  convolution  integral  according  to  equation  (6.15), 
which  was  compared  with  the  pdf  obtained  by  digitizing  the 
concentration  signals  with  an  HP  3721A  correlator  and  electron¬ 
ically  sorting  them  into  100  equally  spaced  bins.  An  analog 
signal  multiplier  and  vol tage-to-f requency  counter  were  used 
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to  obtain  n,  ,  c  and  c'  .  I ntermi ttency  y  was  measured  by 
triggering  a  frequency  counter  whenever  the  signal  was  within 
a  few  microvolts  of  zero.  Chapter  III  describes  the  measure¬ 
ment  techniques  in  detail. 

A  measure  of  'goodness  of  fit'  is  the  coefficient  of 
2 

determination  r  ^.p,  which  indicates  the  amount  of  observed  con¬ 
centration  pdf  curve  'shape'  accounted  for  by  the  intermittent 

2 

lognormal  model  with  Gaussian  noise.  If  r  ^  >  0.8,  one  may 

conclude  that  the  theory  successfully  reproduces  a  significant 

portion  of  the  observed  variation  in  the  measured  pdf.  A  value 
2 

of  r p ^ .p  =  1  means  that  all  factors  affecting  the  observations 
are  properly  included  in  the  theory.  The  coefficient  of  deter¬ 
mination  is  calculated  from 


[NEx.y.  -  Zx.Zy.] 


pdf  [NEx?  -  (Ex1)2]2[NZyJ  -  (Ey.)2]2 


(6.16) 


where  x .  is  the  value  of  the  observed  pdf  P 0 b s  ^ s i ^  a  ^  concen_ 
tration  s.5  and  y.  is  the  corresponding  value  of  the  theoretical 
pdf  p ( s • )  calculated  from  Eq.  (6.15).  Summation  is  made  over 

all  values  of  s .  for  which  either  P  0  5  s  (  s  -j )  or  P(sj)  1  s  non- 

zero,  with  i  =100  corresponding  to  the  total  number  of  data 

max 

points  generated  for  an  observed  pdf. 

The  comparison  of  Eq.  (6.15)  to  data  yields  values  of 

r2Jjr  which  are  tabulated  against  distance  from  source,  height 
pdf 


above  ground  and  lateral  distance  from  plume  centerline  in 
Table  6.1  From  this  data  it  is  apparent  that: 
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•  as  one  moves  out  along  the  plume  centerline,  there  is 

a  definite  improvement  in  the  ability  of  the  intermittent 
lognormal  model  to  describe  the  observed  concentration 
fluctuations,  and 

•  there  is  a  tendency  for  the  theory  to  fit  observations 
better  at  the  plume  edges  than  on  the  plume  centerline. 

A  direct  graphical  comparison  of  the  observed  and 
theoretical  pdfs  is  made  in  Figs.  6.1  through  6.6;  indepen- 
ant  measurements  of  c,  c‘  and  y  were  used  to  compute  the 
theoretical  pdf  shape  parameters  and  cq .  The  figures  show 
that,  in  general, 

•  for  the  locations  where  measurements  were  made 

(1  <  x/ALu  <  4),  the  most  probable  concentrations  are 
about  half  as  large  and  occur  50%  more  frequently  than 
those  predicted  by  the  intermittent  lognormal  theory. 

•  small  concentrations  (<  c/4)  occur  more  frequently 
than  accounted  for  by  the  theory, 

•  midsize  concentrations  (>  c/4  and  <  2c)  occur  less  fre¬ 
quently  than  expected  from  the  theory,  and 
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t  the  frequency  of  occurrence  of  large  concentrations 

(>  2c)  is  well  predicted  by  the  intermittent  lognormal 
model  . 

It  must  be  emphasized  that  these  results  apply  only  for  down¬ 
wind  distances  at  which  measurements  were  made.  Some  other 
points  to  note  are: 

•  a  small  error  in  measured  i nt ermi t tency  can  produce  a 
large  error  in  the  calculated  pdf  of  low  concentrations; 

t  if  the  intermi ttency  factor  is  close  to  unity,  the  inter¬ 
mittent  lognormal  model  is  a  good  fit  to  concentration 
fluctuations  at  ground  level. 

6 . 4  The  Approach  to  Lognormality 

The  most  important  result  of  the  comparisons  in  Figs. 
6.1  to  6.6  is  the  uncovering  of  systematic  differences  between 
theory  and  observations.  The  flow  passing  through  the  sensor 
contains  more  low  concentrations  than  can  be  explained  by  the 
intermittent  lognormal  theory.  This  could  be  evidence  of 
probe  induced  flow  distortion  causing  rapid  diffusion  inside 
the  probe.  However,  probe  induced  distortion  of  the  pdf  does 
not  explain  the  strong  improvement  in  coefficient  of  determin¬ 
ation  with  distance  along  the  plume  centerline.  The  observed 
improvement  is  consistent  with  the  fact  that  the  opportunity 
for  plume  eddy  dilutions  increases  with  travel  time.  The  log¬ 
normal  assumption  requires  that  the  number  of  eddy  dilutions 
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be  large,  and  therefore  the  validity  of  the  theory  (and  the 

2 

magnitude  of  r^^)  s h ou  1  d  improve  with  downwind  distance,  as 
was  observed. 

Another  possible  explanation  for  the  difference  between 

theory  and  experiment  is  that  the  independence  of  dilution 

steps  assumed  by  the  lognormal  model  may  not  prevail  until 

the  diffusion  time  becomes  large  compared  to  the  Lagrangian 

time  scale  T^.  From  Chapter  IV,  the  observed  Euler i an  time 

scale  TEc  of  concentration  fluctuations  was  about  1.5  ms,  so 

that  for  ?Ec  =  3T^c  where  8  -  4,  the  Lagrangian  time  scale  was 

about  6  ms.  With  an  eddy  convection  velocity  just  over  8  m/s, 

the  Lagrangian  length  scale  A^c  was  5  cm.  The  requirement  that 

diffusion  time  be  large  compared  to  T^c  is  therefore  satisfied 

for  downwind  distances  greater  than  about  10  A^c  =  50  cm 

(x/Alu  -  3).  This  means  that  most  of  the  pdf  measurements  were 

made  relatively  close  to  the  source,  where  the  diffusion  time 

was  not  much  larger  than  the  Lagrangian  integral  time  scale  of 

eddy  diffusion.  This  would  help  to  explain  the  observed 

increase  of  r  with  distance  along  the  plume  centerline. 

pdf 

It  is  useful  to  have  a  parameter  which  indicates  when 
the  intermittent  lognormal  distribution  can  be  expected  to  fit 
a  set  of  concentration  fluctuations.  One  possibility  is  to 
consider  the  degree  of  plume  dilution  D.  The  Gaussian  relation 
for  pollutant  concentration  along  the  centerline  of  an  axi- 
symmetric  plume  is 

=  —9-=  >  (6.17) 


£ 

c 
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or,  in  terms  of  a  dilution  factor  D, 


(6.18) 


where  D  is  proportional  to 


and  is  computed  from 


Q/As  cs 

D  =  -0-1  •  r5-  •  (6.19) 

c 


The  coefficient  of  determination  r  ,r  between  the 

pdf 

observed  pdfs  and  the  intermittent  lognormal  theory  has  been 
plotted  against  dilution  factor  D  in  Figure  6.7.  There  is  a 
strong  trend  for  large  dilution  factors  to  be  associated  with  a 
high  degree  of  correl at i on  between  theory  and  experiment.  This 
reflects  the  fact  that  the  approach  to  lognormality  requires 
many  dilution  steps,  in  which  case  the  dilution  factor  must  be 
1 arge . 


It  is  important  to  note  that  care  must  be  taken  in 
choosing  an  eddy  convection  velocity  U  =  x/t  for  calculating  D 
from  Eq.  (6.19),  especially  when  flow  patterns  are  complex,  such 
as  in  the  neighborhood  of  hills  and  buildings.  Fortunately, 
the  eddy  travel  time  t  is  frequently  underestimated  in  complex 
flows.  This  can  produce  an  underestimation  of  D,  and  possibly 
a  decision  not  to  use  the  intermittent  lognormal  theory  when 
in  fact  it  could  be  used.  An  overestimation  of  D  can  have  the 
more  serious  consequence  of  providing  encouragement  to  use  the 
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lognormal  theory  in  regions  where  it  poorly  fits  the  concen¬ 


tration  pdf. 

The  coefficient  of  determination  r2^  gives  only  a 

qualitative  measure  of  the  accuracy  to  be  expected  from  the 

intermittent  lognormal  model  for  concentration  fluctuations. 

In  practice  the  coefficient  of  determination  r2  alone  tells 

pd  t 
o 

very  little.  For  example,  what  rpdf  (or  equivalently,  D)  is 

needed  to  give  calculated  peak-to-means  that  are  within  20% 

of  the  observed  values?  A  practical  answer  to  this  question 

is  given  in  the  next  section,  where  it  is  shown  that  a 

2 

lower  limit  to  the  ratio  D/i  must  be  exceeded  in  order  to 

0 

ensure  accurate  peak-to-mean  calculations. 


6 . 5  Peak-to-Mean  Concentration  Ratios 

The  'peak'  concentration  at  a  receptor  is  defined  as 
that  value  which  is  exceeded  less  than  a  fraction  P  of  the 
time.  The  probability  density  function  for  concentration 
fluctuations  may  be  used  to  calculate  the  peak-to-mean  concen¬ 
tration  ratios  at  various  receptor  locations.  If  cp  is  such  a 
peak,  then,  from  Eq.  (6.5), 


tP(cp)  =  (1  -  P)  =  1 


1  +  erf [ 


1  n (c  / c 
P 

/T  aL 


(6.20) 


Solving  for  c  /c  gives  the  peak-to-mean 

P 


ratio 


■  e 
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zP~  =  ^-exp{2^[2 JT erf'1  (l-^.)-aL]}  ,  (6.21) 

where  cL  is  computed  from  Eq.  (6.8).  Csanady  (1973,  p.  232) 
has  developed  similar  equations  for  the  case  where  y  =  1. 

The  peak-to-mean  ratios  which  would  be  exceeded  10%  of 
the  time  (P  =  .10)  along  the  plume  centerline  were  calculated 
using  Eq.  (6.21),  and  then  compared  with  the  values  obtained 
by  numerically  integrating  the  observed  pdfs. 

The  ratio  of  theoretical  to  observed  pea k-to-means  in¬ 
creased  towards  unity  as  dilution  factor  D  became  larger,  but 
decreased  with  distance  off  the  plume  centerline.  This  is 
evidence  that  some  other  parameter  besides  0  plays  a  role  in 
the  approach  to  lognormal i ty .  Inspection  of  Eqs.  (6.8)  and 
(6.21)  suggests  that  the  square  of  concentration  fluctuation 

intensity,  i  ,  has  a  suitably  strong  influence  on  computed 

0 

? 

peak-to-means ,  and  Ch.  V  shows  that  i  increases  rapidly  with 
distance  off  plume  centerline.  These  observations  prompted  the 
plot  of  theoretical  to  observed  peak-to-mean  ratio  against  the 
composite  dilution  parameter  D / i c  >  presented  in  Fig.  6.8.  The 
data  falls  approximately  on  a  straight  line  against  this  para¬ 
meter,  and  indicates  that  calculated  peak-to-mean  ratios  should 
be  accurate  to  within  20%  at  any  downwind  locations  for  which 
q  /  i  is  greater  than  400.  For  elevated  point  sources  in  the 

atmospheric  boundary  layer,  the  reflected  Gaussian  plume  model 
is  a  good  choice  for  estimating  D;  the  value  of  i c  c  a  n  he 
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obtained  from  measurements,  or  from  the  dissipated  Gaussian 

2 

model  proposed  for  c'  in  Chapter  V. 

The  intermittent  lognormal  pdf  fits  the  data  fairly 

well  close  to  the  ground  (see  Fig.  6.6),  and  Eq.  (6.21)  with 

aL  derived  from  measurements  has  been  used  to  calculate  the 

peak-to-mean  ratios  exceeded  10%  of  the  time  (i.e.  P  =  .10) 

along  the  ground  under  the  plume  centerline.  Fig.  6.9  shows 

that  the  peak-to-mean  concentration  ratio  for  P  =  .1  is  maximum 

at  a  distance  of  0.6  x  ,  where  x  is  the  distance  to  the 

in  a  x  max 

point  of  maximum  mean  concentration.  However,  the  largest 
actual  peak  concentrations  occur  farther  out,  at  a  distance  of 
xpeak  =  xmax  ^  p  =  •  1  •  This  distance  becomes  less  as  P 

is  reduced.  For  example,  if  P  =  .01  then  the  largest  short¬ 
term  concentrations  are  about  nine  times  larger  than  c  and 
occur  near  0.5  x( 
a  simple  rule  of  thumb  is  to  assume  that  peak  concentrations 
occur  at  two  thirds  of  the  distance  out  to  the  point  of  maximum 


max 


Since  Cpea^  changes  slowly  near  x  a^. 


mean  concentration,  or 


XPe-a-!l  =  |  (6.22) 

X  3 

max 

for  values  of  P  between  0.01  and  0.1. 

The  source  height  was  not  varied  during  these  experi¬ 
ments,  so  no  data  is  available  to  check  if  relation  (6.22) 
applies  to  a  range  of  source  heights.  It  is  possible  to  vary 
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2 

the  source  height  in  the  Gaussian  model  for  cl  presented  in 
Chapter  V.  The  results  could  then  be  used  in  Eq .  (6.21)  to 
establish  the  dependance  of  x  ,  on  source  height. 
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Intermittent  Lognormal 
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Experiment  (3  Runs) 


194 


195 


UOUBU.IIUJ810CI  JO  JU8PIJJ80Q  ipdzJ 


FIG.  6.7  Fit  of  Intermittent  Lognormal  PDF  Theory  vs.  Dilution  Factor 
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FIG.  6.8  Ratio  of  Calculated  to  Observed  Peak-to-Means 
(for  P=.l)  vs.  Dilution  Parameter  D/if 
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FIG.  6.9  Location  of  Highest  Mean  and  Fluctuating  Concentrations 
Near  Ground  Level  for  4000:l  Model  Scale  Factor 


CHAPTER  VII 


SUMMARY  AND  CONCLUSIONS 
7 . 0  Summary  of  Results 

Knowl edge  of  the  statistical  and  spatial  distributions 
of  fluctuating  concentration  in  elevated  plumes  is  necessary 
before  the  frequency  of  pollution  incidents  can  be  quantita¬ 
tively  established  for  the  purpose  of  hazard  assessment. 

This  study  uses  a  wind  tunnel  simulation  of  dispersion  in  the 
atmospheric  boundary  layer  to  provide  new  information  on 
concentration  fluctuations  in  plumes. 

An  aspirated  single  hot-film  sensor,  driven  by  a 
constant  temperature  anemometer,  was  developed  to  observe 
high-frequency  fluctuations  of  helium  tracer  concentration  in 
an  atmospheric  boundary  layer  wind  tunnel.  This  fast-response 
probe  produced  about  1  mv  output  signal  for  a  2  ppth  helium- 
air  mixture,  and  could  follow  helium  fluctuations  up  to  530 
Hz,  at  which  point  the  probe  output  was  3  db  down.  The  noise 
level  in  the  probe  signal  corresponded  to  an  rms  concentration 
fluctuation  of  0.2  ppth  about  a  mean  of  zero  ppth.  Undistorted 
measurements  of  concentration  variance,  probability  density 
function,  autocorrelations  and  spectra  were  obtained  with  the 
aspirated  probe;  measured  values  of  concentration  time 
derivative  required  correction  for  probe- induced  dissipation 
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of  the  high-frequency  signal  components.  Mean  concentration 
was  measured  with  an  existing  high  accuracy  slow  response 
detector,  since  zero-drift  problems  required  removal  of  the 
mean  component  from  the  fast  response  probe  signal  by  capaci¬ 
tive  AC  coupling  of  the  probe  output. 

Standard  modeling  techniques  were  used  to  stimulate 
the  rapid  growth  of  a  thick  turbulent  boundary  layer  within 
the  confines  of  a  short  wind  tunnel.  A  passive,  non-bo uy ant 
plume  of  helium  tracer  gas  was  released  above  the  tunnel  floor 
at  the  entrance  to  the  test  section.  The  wind  tunnel  boundary 
layer  simulated  a  neutrally  stratified  atmospheric  boundary 
layer  with  a  model  scale  factor  of  4000:1,  which  was  verified 
by  comparison  of  measured  tracer  plume  spreading  rates  with 
those  reported  for  the  full  scale.  With  this  model  scale 
factor,  the  helium  tracer  represented  a  chimney  plume  which 
had  leveled  off  at  a  height  of  120  m. 

Vertical,  lateral  and  alongwind  profiles  of  concentra¬ 
tion  i  n termi t tency  were  measured  in  the  tracer  plume.  The 
profile  shapes  were  well  described  by  the  Gaussian  error 
function  model  originally  suggested  by  Corrsin  and  Kistler  (1 954) 
for  turbulence  intermittency  in  round  free  jets.  This  inter- 
mittency  model  worked  well  everywhere  in  the  plume  except  near 
the  ground,  which  acted  to  constrain  the  Gaussian  plume  boundary 

movements  upon  which  the  model  is  based. 

Vertical  and  lateral  profiles  of  mean  tracer  concen¬ 
tration  were  fitted  to  the  classic  reflected  Gaussian  equation 
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often  used  to  describe  the  mean  concentration  field.  This 
model  reproduced  the  present  measurements  with  relative  error 
less  than  10%  of  the  local  concentration  except  at  the  plume 
edges  and  near  the  ground,  where  errors  of  20-40%  could  be 
observed . 

Profiles  of  concentration  variance  were  measured  in 
the  tracer  plume,  and  were  well  described  by  a  new  'dissipated 
Gaussian'  model  developed  during  this  study.  In  addition  to 
providing  a  good  fit  to  data,  the  model  is  easy  to  use  because 
its  shape  parameters  are  the  familiar  plume  spreads  a and  a z 
which  appear  in  the  reflected  Gaussian  model  for  mean  concen¬ 
tration.  However,  the  dissipated  Gaussian  model  for  fluc¬ 
tuations  has  two  disadvantages.  First,  it  contains  a  'source 

variance'  term  s  about  which  not  much  is  known;  s  is 

o  o 

constant  through  the  plume,  but  probably  varies  with  atmos¬ 
pheric  turbulence  levels.  Second,  the  model  assumes  comp  1 e te 
dissipation  of  concentration  flucutations  at  the  ground.  Total 
dissipation  is  unlikely  to  occur  in  practice,  so  it  may  be 
necessary  to  introduce  an  empirical  'dissipation  coefficient 
in  the  model.  The  present  measurements  lack  sufficient  detail 

to  resolve  this  question. 

The  reflected  Gaussian  model  has  been  used  in  the 
diffusion  equation  for  c  to  derive  an  expression  for  the 
vertical  eddy  diffusivity  under  non-uniform  wind  conditions. 
This  eddy  diffusivity,  which  relates  Reynolds  flux  -wlc'  to 
the  vertical  gradient  of  mean  concentration,  was  compared  with 


’  a rtr-*  w  m  $  *4  b  »  '*»*  ****  *  ' '  q  '’ ? D 1  1 
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the  vertical  eddy  diffusivity  obtained  by  using  measured 
data  in  the  diffusion  equation.  The  experimentally-derived 
vertical  profiles  of  eddy  diffusivity  were  complicated,  having 
a  discontinuity  near  the  source  height.  Nevertheless,  for 
non-uniform  winds  the  Gaussian  mean  concentration  model  was 
able  to  successfully  predict  these  diffusivity  shapes. 

An  exami nat i on  of  terms  in  the  variance  conservation 

equation  led  to  the  observation  that  measured  values  of 

variance  production  and  dissipation  were  not  in  local  balance. 

As  the  plume  moved  downwind,  production  became  negligible 

compared  with  the  local  dissipation  of  concentration  variance 
2 

c1  .  This  in  turn  led  to  an  approximate  balance  equation  for 
concentration  variance  which  has  the  same  form  as  the  diffusion 
equation  for  mean  concentration.  The  dissipated  Gaussian 
variance  model  was  shown  to  be  a  solution  to  this  approximate 
conservation  equation. 

The  measured  variance  dissipation  term  was  found  to  be  pro¬ 
portional  to  the  local  concentration  variance,  which  confirmed 
Csanady's  (1967)  supposition  that  this  should  be  so.  In 
addition,  the  measurements  supported  Csanady's  (1967)  assumpt¬ 
ion  that  this  proportionality  factor  is  the  reciprocal  of  a 
decay-time  scale  which  is  constant  across  a  given  plume  section 
and  which  increases  linearly  with  distance  from  the  source. 

Using  the  concept  of  independant  zero  and  non-zero 
concentration  eddy  populations,  an  existing  lognormal  concen¬ 
tration  probability  model  was  modified  to  include  intermittent 
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concentrations.  Measured  or  calculated  values  of  concentra¬ 
tion  mean,  variance  and  i n termi tt ency  can  be  used  to  compute 
the  shape  parameters  of  this  intermittent  lognormal  model 
for  the  concentration  probability  density  function.  The 
model  was  compared  with  wind  tunnel  measurements  of  the  pdf 
at  several  points  within  the  plume,  and  predicted  peak-to-mean 

concentration  ratios  within  20%  of  observed  values  so  long  as 

2 

a  dilution  parameter  D/i  was  greater  than  400.  The  highest 
concentration  peaks  were  observed  to  occur  about  two  thirds 
of  the  distance  out  to  the  point  of  maximum  concentration. 

7 . 1  Recommendations  for  Further  Study 

There  are  several  aspects  of  the  fluctuation  model 
that  require  further  investigation.  These  can  be  listed  by 
model  component  as  follows: 

1 .  Plume  I ntermi ttency 

•  Additional  work  is  necessary  to  provide  an  under¬ 
standing  of  the  relationship  between  inter- 
mittency  profile  parameters  yi/2»  zv2an<^ 

plume  spreads  o  9 

•  The  intermittency  model  should  be  refined  to 
account  for  the  constraining  influence  of 
the  ground.  A  truncated  Gaussian  model  of 
plume  motion  near  the  ground  may  be  useful 
in  this  regard. 
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•  More  measurements  are  needed  to  determine  how  the 
centerline  interim*  ttency  depends  on  atmospheric 
turbulence  and  source/sensor  characteristics. 

2 .  Concentration  Variance 

•  The  manner  in  which  the  source  variance 
parameter  sq  depends  on  atmospheric  turbulence 
must  be  established  before  the  dissipated  Gaussian 
model  can  provide  quantitative  estimates  of  fluctua¬ 
tion  intensity  in  full-scale  dispersion  situations. 

•  Detailed  measurements  of  concentration  variance 
at  and  near  the  ground  are  needed  in  order  to 
assess  the  importance  of  incomplete  surface 
dissipation  of  concentration  fluctuations. 

3 .  Intermittent  Lognormal ity 

•  Full  scale  concentration  data  should  be  used 
to  test  the  statistical  model. 

§  The  region  of  model  applicability  within  the 
plume  should  be  thoroughly  investigated. 

•  For  full-scale  measurements,  long-term  average 
pdfs  seem  to  be  exponential,  while  short-term 
average  pdfs  appear  log-normal.  This  apparent 
difference  of  form  requires  more  study. 


204 


The  measurements  needed  to  guide  model  refinements  for 
improved  performance  near  the  ground  can  be  obtained  only  if 
the  fast  response  probe  is  modified.  Specifically,  the  probe 
must  be  redesigned  for  sampling  near  and  flush  to  solid 
boundaries.  Also,  an  alternative  to  the  filter  technique 
used  to  remove  pressure  fluctuations  should  be  found,  in  order 

to  avoid  pro  be -induced  dissipation  of  the  high-frequency  con¬ 
centration  fluctuations  important  to  measurements  of  natural 
dissipation  rates. 

/ 

7 . 2  Concluding  Remarks 

This  study  has  developed  and  tested  a  model  for  con¬ 
centration  fluctuations  in  plumes,  which  can  be  used  to  predict 
the  fraction  of  time  that  a  specified  pollutant  concentration 
is  exceeded  downwind  of  an  industrial  plant.  This  model  can 
aid  legislators  in  the  formulation  of  more  realistic  constraints 
on  the  emission  of  waste  gases  by  industry.  Until  such  new 
regulations  are  put  in  effect,  government  and  industry  will 
continue  to  be  faced  with  unavoidable  violations  of  air  quality 

standards  which  ignore  the  inherently  random  variability  of 
pollutant  concentrations  dispersed  by  the  turbulent  atmosphere. 


. 
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APPENDIX  A 


MEASUREMENT  OF  TURBULENT  SCALAR  FLUCTUATIONS 


Introduction 

The  laboratory  measurement  of  high-frequency  scalar 
fluctuations  in  turbulent  flows  has  attracted  little  attention 
during  the  last  several  decades.  Scalar  properties  of  interest 
have  usually  been  the  fluid  temperature  or  species  concentration 
The  lack  of  data  is  due  in  large  part  to  the  experimental 
difficulties  encountered  in  making  the  measurements. 

Experimental  problems  encountered  generally  include: 
detector  sensitivity  to  extraneous  noise  sources  (environmen¬ 
tal  or  electronic),  insufficient  frequency  response,  inability 
to  obtain  a  dynamic  calibration,  cumbersome  or  expensive  designs 
and  toxic,  explosive  or  messy  tracer  materials.  Of  this  short 
list,  sensitivity  to  extraneous  noise  sources  is  probably 
the  most  often  encountered.  Environmental  noise  refers  to 
factors  in  the  sensor's  environment  which  produce  an  output 
signal,  but  which  are  not  of  immediate  interest.  Depending 
on  what  is  being  measured  and  what  detection  system  is  being 
used,  disturbing  environmental  factors  can  include  fluctuations 
in  laboratory  air  temperature,  moisture  content,  acoustic 
noise  level,  tobacco  smoke  content,  air  turbulence,  air 
velocity,  lighting  systems  and  the  nearby  operation  of  high 

voltage  electronic  equipment. 
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Because  of  the  random  nature  of  turbulent  fluctuations, 
random  electronic  noise  from  the  signal  generation  or  pro¬ 
cessing  equipment  can  sometimes  obscure  the  information 
contained  in  the  signal.  Even  when  the  signal-to-noise  ratio 
is  adequate  for  some  parts  of  an  experiment,  there  can  often 
be  large  regions  of  the  turbulent  flow  field  which  are  closed 
to  investigation  because  the  associated  data  signal  is  too 
weak  to  appear  through  the  noise. 

Light-Scattering  Techniques 

In  spite  of  these  problems,  which  are  after  all  not 
unique  to  the  measurement  of  turbulent  scalar  fluctuations, 
several  measuring  systems  have  been  developed  which  can  give 
satisfactory  results.  Rosensweig  (1959,  1961)  developed  a 
sea ttered -1 i g h t  technique  for  studying  the  concentration 
field  of  suspended  colloidal  particles  marking  one  of  the  fluid 
streams  entering  a  mixing  field.  A  linear  response  was  obtain¬ 
ed  to  both  the  time-mean  and  the  time-variant  components  of 
material  concentration  at  a  point,  in  a  sheet,  or  in  a  volume, 
according  to  the  modes  of  illumination  and  observation.  The 
capabilities  and  limitations  of  the  method  were  reviewed  by 
Becker,  Hottel  &  Williams  (1967).  The  technique  gives 
excellent  results  in  the  convection  region  of  the  turbulence 
spectrum  but  fails  where  molecular  diffusion  is  important. 

This  is  because  the  tracer  particles  have  a  Schmidt  number 
Sc  =  v/D  of  almost  4  x  104  and  therefore  diffuse  much  more 
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slowly  than  the  molecules  of  the  stream  into  which  they  are 
introduced.  A  Schmidt  number  of  0.74  is  appropriate  for 
diffusion  of  air  into  air.  The  bulky  nature  of  the  optical 
system  prevented  extensive  use,  even  in  large  wind  tunnels. 

This  problem  was  overcome  when  Liu  and  Karaki  (1972)  intro¬ 
duced  fibre  optics  as  flexible  light  transmitters  for  the 
high-intensity  mul tichromatic  light  used. 

Yang  and  Meroney  (1974,  1975)  have  further  refined  the 
light-scattering  technique  by  employing  a  5  mW,  He-Ne  gas 
laser  as  the  incident  light  source.  The  great  advantage  of 
using  a  laser  beam  over  an  incandescent  light  source  is  that 
the  monochromatic  characteristics  of  a  laser  beam  permit 
doppler  techniques  to  be  applied  to  the  scattered  light.  From 
such  analysis  the  instantaneous  velocity  fluctuation  components 
can  be  abstracted,  while  the  total  scattered  energy  yields  the 
concentration  fluctuations.  Thus,  there  is  the  possibility 
of  combining  the  two  signals  in  order  to  measure  the  three 
unknown  correlation  terms  u'c',  v'c'  and  w'c1,  where  u',v' 
and  w'  are  the  fluctuating  velocity  components  and  c'  is  the 
concentration  fluctuation. 

Yang  and  Meroney's  probe  was  about  20  cm  long,  8  cm 
high  and  about  2  cm  thick,  whereas  the  volume  sampled  was  on 
the  order  of  only  1  mm^.  In  spite  of  the  small  sample  volume, 
a  probe  this  size  is  too  large  for  use  in  small  wind  tunnels 
or  for  measurements  near  a  finely  detailed  wind  tunnel  model. 
Reductions  in  probe  size  would  result  if  a  smaller  laser  and 
photomultiplier  tube  were  used.  The  sample  volume  has  a  lower 
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size  limited  by  the  signal -to-noi se  ratio  which  decreases 
with  sample  volume. 

Motycka  and  Leutheusser  (1972)  have  constructed  a 
very  small  light  scattering  probe  which  is  only  1  cm  long 

and  which  samples  the  average  concentration  within  a  volume 

3 

of  40  mm  .  The  probe  contains  a  miniature  solid-state  infra¬ 
red  light  source  and  a  silicon  photo-detector;  the  time 
constant  for  the  system  is  stated  to  be  about  2  ms,  which  is 
not  particularly  fast.  Presumably  the  dynamic  calibration 
was  done  by  suddenly  exposing  the  photo-detector  to  the  light 
from  the  IR  source,  and  examining  the  output  signal  on  an 
oscilloscope.  If  such  was  the  case,  the  response  of  the 
probe  is  probably  limited  by  the  response  characteristics  of 
the  photo-receptor.  A  faster  response  time  and  a  smaller 
sample  volume  should  be  possible  with  further  development. 

The  advantages  of  a  light-scattering  system  are  power¬ 
ful.  With  proper  design,  the  measurements  are  obtained 
entirely  remotely  from  an  undisturbed  flow  field.  The  fre¬ 
quency  response  of  most  systems  is  for  all  practical  purposes 
limited  only  by  the  inertia  of  the  aerosol  particles.  Measure¬ 
ments  can  be  obtained  from  sample  points,  sheets  or  volumes, 
and  dual  systems  allow  spatial  correlations  of  concentration 
fluctuation  to  be  measured.  The  size  of  the  probe  limits  the 
space  into  which  the  probe  can  be  inserted,  and  as  noted  above, 
the  system  cannot  observe  the  dissipation  of  fluid  concentra¬ 
tion  fluctuations  by  molecular  diffusion.  Finally,  as 
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mentioned  previously,  a  laser-doppler  system  operated  simult¬ 
aneously  will  permit  calculation  of  the  unknown  velocity- 
concentration  correlation  terms,  required  for  the  study  of 
turbulent  diffusivity. 

Temperature  Fluctuations 

The  measurement  of  turbulent  temperature  fluctuations 
has  met  with  some  success,  and  has  provided  useful  data  on 
turbulent  scalar  fluctuations.  Crum  and  Hanratty  (1965)  used 
hot-wire  anemometers  run  at  low  probe  current  (5  ma)  to  detect 
only  temperature  fluctuations  and  not  velocity  fluctuations. 
These  probes  responded  fully  to  frequencies  in  excess  of  7  kHz. 
For  flow  regions  having  low  signal -to-noi  se  ratios,  a  0.25  mm 
bead  thermistor  was  used  in  place  of  the  hot-wire.  Fiedler 
(1974)  used  a  DISA  55  F05  resistor  probe  which  responded  to 
temperature  fluctuations  up  to  2  kHz.  The  main  advantage 
of  measuring  temperature  fluctuations  rather  than  species 
concentration  is  that  fairly  standard  anemometer  technology 
(i.e.  sensors  and  electronic  equipment)  and  techniques  can  be 
used.  Because  of  the  high  frequency  response  and  small  spatial 
resolution  available,  it  is  possible  to  extend  the  measurements 
well  into  the  dissipative  range  of  eddy  sizes.  However, 
there  are  several  disadvantages  which  must  also  be  considered. 
Acceptable  s ignal -to-noi  se  ratios  require  either  a  sensitive 
probe  or  a  powerful  heat  source.  A  probe  with  sufficient 
sensitivity  may  also  have  unwanted  sensitivity  to  velocity 
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fluctuations,  or  poor  frequency  response.  Sufficiently 
powerful  heat  sources  may  be  impossible  to  provide  without 
creating  undue  flow  disturbance.  Also,  because  of  surface 
heat  transfer,  any  solid  boundaries  in  the  flow  will  act  as 
sinks  for  the  diffusing  scalar,  and  for  a  large  group  of 
dispersion  experiments  this  is  undesireabl e. 

Hot-Wire  Anemometers 

Becker  and  Booth  (1975)  have  stated  that  hot-wire 
anemometry  is  the  only  proven  technique  presently  available 
for  studying  dissipation  in  systems  of  laboratory  scale.  They 
state  further  that  the  only  scalar  quantity  which  is  easy  to 
detect  by  this  means  is  temperature.  There  is  no  doubt  about 
the  truth  of  this  statement,  as  demonstrated  by  the  efforts 
of  Tombach  (1969)  to  apply  hot-wires  to  the  measurement  of 
velocity  fluctuations  in  a  turbulent  concentration  field. 

He  encountered  difficulty  with  the  hot-wire  technique  because 
of  a  "historical"  effect  somewhat  akin  to  hysteresis,  which 
occurred  after  the  sensor  had  been  exposed  to  high  concentra¬ 
tions  of  helium.  This  phenomenon  was  not  a  problem  for  Way 
and  Libby  (1970,  1971)  in  their  work  using  a  dual-probe  combin¬ 
ation  of  hot-wire  and  hot-film  sensors.  They  attributed 
success  to  a  mode  of  calibration  and  data  collection  which  did 
not  require  exposure  of  the  probe  to  high  helium  concentrations 

for  long  periods  of  time. 


216 


As  mentioned  previously,  a  disadvantage  of  hot-wires 
or  films  is  their  sensitivity  to  environmental  factors  such 
as  air  velocity,  temperature  and  species  concentration.  There 
are  several  ways  around  this  problem:  the  sensitivity  to 
unwanted  fluctuations  must  be  eliminated,  the  unwanted 
fluctuations  must  somehow  be  removed  prior  to  reaching  the 
sensor,  or  the  effect  of  the  different  fluctuations  must  be 
explicitly  accounted  for.  This  last  method  has  the  advantage 
of  yielding  the  most  information  about  the  turbulent  flow 
field,  and  was  the  choice  of  Way  and  Libby.  Their  dual-sensor 
probe  separated  velocity  fluctuations  from  concentration 
fluctuations  by  taking  advantage  of  the  different  response 
characteristics  of  hot-films  and  hot-wires.  Briefly,  the 
method  requires  calibration  curves  for  each  sensor,  with  first 
velocity  and  then  concentration  being  held  constant.  With 
proper  probe  geometry,  an  observed  voltage  pair  from  the  two 
sensors  will  correspond  to  a  unique  combination  of  velocity 
and  concentration;  a  set  of  probe  calibration  equations  is 
used  to  calculate  velocity  and  concentration.  Because  the 
equations  for  probe  response  hold  only  approximately,  and 
because  sensor  interactions  strongly  affected  sensitivity, 
determining  optimum  probe  geometry  is  a  major  problem. 

The  main  advantage  of  the  dual-sensor  probe  is  its 
ability  to  measure  the  fluctuating  velocity  and  concentration 
simultaneously,  with  a  frequency  response  and  spatial  resolu¬ 
tion  sufficient  to  observe  the  dissipative  eddies.  It  is  then 
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possible  to  directly  measure  and  study  the  velocity -concen¬ 
tration  correlation  terms  u  1  c  1  ,  v  1  c  1  and  w'c1  for  all  eddy 
sizes  of  interest.  A  triple-sensor  configuration  using  an 
X-wire  and  hot  film  would  permit  measurement  of  two  instan¬ 
taneous  velocity  components  rather  than  one. 

The  most  noticeable  disadvantage  of  the  system  is 
the  large  amount  of  digital  data  reduction  required  to  extract 
the  velocity  and  concentration  fluctuations  from  the  observed 
probe  voltages.  Way  and  Libby  found  continuous  on-line  data 
analysis  was  not  possible  because  it  took  six  seconds  of 
computer  time  to  process  one  second  of  raw  data.  Advances 
in  micro-processor  technology  are  such  that  real-time  data 
analysis  should  be  possible  in  the  near  future.  A  disadvan¬ 
tage  of  more  permanency  is  the  extreme  care  required  during 
calibration,  testing  and  data  reduction  in  order  to  ensure 
satisfactory  accuracy  in  the  results. 

Measurements  in  Liquids 

The  turbulence  fluctuation  sensors  described  so  far 
were  operated  only  in  gaseous  mediums.  Turbulent  scalar 
fluctuations  exist  in  liquid  systems  also,  and  have  been  the 
subject  of  laboratory  investigations.  Gibson  and  Schwartz 
(1963a,  1963b)  describe  the  design,  construction,  calibration 
and  use  of  a  single-electrode  conductivity  sensor  which  detects 
very  small  amplitude  fluctuations  of  concentration  or  temper¬ 
ature  in  aqueous  electrolytic  solutions.  The  spatial 


:* 


218 


resolution  of  the  probe  is  on  a  par  with  that  available  from 
hot-wires  or  films,  and  is  superior  to  that  of  a  light-scatter 
probe.  The  detector  is  capable  of  following  fluctuations  at 
frequencies  of  several  thousand  Hertz.  The  conductivity 
sensor  is  insensitive  to  velocity  fluctuations,  but  is  sen¬ 
sitive  to  fluctuations  of  both  temperature  and  scalar  concen¬ 
tration.  It  is  simple  to  statically  calibrate  the  detector 
for  scalars,  and  to  design  experiments  where  first  temperature 
and  then  concentration  is  the  scalar  being  observed.  This 
easily  controlled  dual  sensitivity  makes  the  sensor  very 
versatile. 

Gibson  and  Schwartz  do  not  describe  their  calibration 
procedure,  but  steady-state  calibration  is  likely  to  be 
straightforward.  A  dynamic  calibration  of  the  conductivity 
probe  appears  to  be  much  more  difficult;  the  theoretical 
response  of  a  spherical  electrode  to  a  sinusoidal  concentration 
fluctuation  was  used  to  estimate  performance  of  the  actual 
hemispherical  probe. 

Salt  concentration  fluctuations  (Sc  -  700)  diffuse 
much  more  slowly  than  temperature  fluctuations  (Sc  -  7).  Thus 
the  conductivity  sensor  has  some  of  the  limitations  of  the 
light-scattering  probe  in  that  the  salt  tracer  (like  the  oil 
fog)  sustains  eddy  sizes  much  smaller  than  could  exist  if  the 
tracer  had  properties  more  like  those  of  the  medium.  This 
problem  is  particularly  important  when  trying  to  measure  the 
modelled  fluctuations  of  air  pollutants. 
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Summary 

Five  techniques  for  laboratory  measurement  of 
turbulent  scalar  fluctuations  have  been  described.  None  of 
the  detectors  are  perfect:  their  relative  strengths  and 
weaknesses  often  make  one  system  perform  better  than  others 
in  a  given  experiment.  The  advantages  and  disadvantages  of 
the  available  detectors  are  compared  in  Table  1.1.  Chapter 
II  describes  the  design,  construction  and  calibration  of 
a  fast  response  helium  concentration  detector  developed  as 
part  of  the  present  work. 


APPENDIX  B 


CONSTRUCTING  THE  HOT-FILM  CONCENTRATION  DETECTOR 

One  of  the  objectives  of  this  study  was  to  construct 
a  fast  response  concentration  detector  that  would  be  small 
enough  for  use  in  a  wind  tunnel  of  cross-section  30  cm  x  30  cm. 
A  hot-wire  or  hot-film  was  proposed  for  the  sensing  element, 
since  their  frequency  response  was  known  to  be  good.  In 
order  to  detect  species  concentration,  the  user  must  exploit 
the  sensitivity  of  a  hot-film  or  wire  to  the  heat  conductance 
of  the  gas  flowing  past  the  element.  Helium  was  chosen  as 
the  tracer  material  because  it  conducts  heat  seven  times 
faster  than  air,  and  was  therefore  easily  detected. 

Also,  it  has  the  desireable  attributes  of  being  clean,  safe, 
inexpensive  and  easy  to  use.  Helium's  natural  bouyancy  is 
expected  to  be  an  advantage  for  the  modelling  of  bouyant 
chimney  plumes. 

The  temperature  of  a  hot-wire  is  sensitive  to  local 
velocity  fluctuations  in  addition  to  fluctuations  in  the 
composition  of  the  surrounding  gas.  For  this  reason  an 
aspirated  probe  was  proposed,  in  which  velocity  over  the  sensor 
would  be  held  constant  regardless  of  the  turbulence  around 
the  probe.  Choked  flow  through  either  an  orifice,  needle 
valve  or  capillary  tube  would  provide  the  necessary  constant 
flow  rate.  The  adjustability  of  a  needle  valve  makes  isokinetic 


220 


221 


sampling  possible,  so  it  was  the  candidate  of  first  choice. 

The  selection  of  position  for  the  choking  device,  either 
upstream  or  downstream  of  the  sensing  element,  is  critical. 

If  the  flow  is  choked  before  it  reaches  the  sensor  (upstream- 
sonic),  there  is  a  risk  that  concentration  fluctuations  will 
be  distorted  or  diffused  as  they  pass  through  the  shock  wave 
at  the  sonic  region.  The  alternative  of  choking  the  flow 
after  it  passes  over  the  sensor  (downstream-sonic)  means  that 
only  the  mean  flow  over  the  sensor  is  held  constant;  turbulence 
external  to  the  probe  may  still  find  its  way  inside  and  affect 
the  probe  output  signal.  It  was  felt  that  sensitivity  to 
external  turbulence  was  the  lesser  of  two  evils,  so  the 
downstream-sonic  probe  configuration  was  chosen. 

A  brass  prototype  was  built  with  a  DISA  miniature 
hot-wire  as  the  sensing  element.  Wind  tunnel  tests  immediately 
revealed  a  serious  problem:  with  a  choked  capillary  downstream 
from  the  hot-wire,  the  background  noise  level  from  the  sensor 
was  much  too  large.  Further  investigation  showed  the  noise 
to  have  three  components:  high  frequency  electronic  noise 
generated  in  the  anemometer  circuitry,  acoustic  noise  produced 
inside  the  probe  due  to  internal  roughness  and  sharp  corners, 
and  wind  tunnel  turbulence  which  found  its  way  into  the  probe. 

A  battery-powered  model  3342  Krohn-Hite  variable  low-pass  active 
filter  was  used  to  eliminate  the  high  frequency  electronic 
noise  (above  1500  Hz).  Removal  of  turbulence  noise  was  more 
difficult.  It  was  found  that  a  short  fibre  plug  (cut  from  a 
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cigarette  filter)  inserted  into  the  probe  inlet  caused  an 
immediate  reduction  in  sensitivity  to  wind  tunnel  turbulence. 
Because  the  filter  material  was  expected  to  distort  or  smear 
the  smallest  concentration  eddies,  a  tri al -and-error  procedure 
was  used  to  obtain  the  shortest  length  and  loosest  packing 
of  filter  material  which  would  adequately  remove  the  wind 
tunnel  turbulence.  Chapter  III  presents  a  method  for  empiri¬ 
cally  correcting  the  probe  signal  to  compensate  for  probe- 
induced  dissipation  of  the  high-frequency  concentration 
fluctuations.  The  last  source  of  noise  was  internally  gener¬ 
ated  turbulence,  and  its  reduction  required  smoothing  of  the 
internal  flow  passages  in  the  probe. 

Along  with  excessive  background  noise,  the  first 
prototype  probe  had  other  problems.  It  was  slow  to  build, 
difficult  to  repair  (i.e.  replacement  of  burnt-out  hot-wires), 
and  it  was  too  large.  These  problems  were  mainly  overcome 
with  the  second  prototype,  which  was  of  modular  construction, 
composed  of  concentric  aluminum  and  brass  tubes  inside  a  1/8" 

ID  plastic  "T".  The  sensing  element  was  changed  to  a  TSI 
1 276-1 OA  hot-film  because  of  its  superior  sturdiness  and  more 
suitable  length.  The  sensor  was  held  firmly  in  the  probe  by 
a  machined  brass  plug  which  fitted  snugly  into  the  plastic  T. 
This  probe  design  was  very  satisfactory  from  many  points  of 
view.  It  was  low  cost  ($75  hot-film  +  $5  materials  +  $10 
machine  shop  labor),  easy  to  put  together  or  repair  (5  hours 
to  build  from  scratch,  2  hours  to  replace  a  burnt-out  hot-film). 
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and  it  was  small  (6  mm  OD,  40  mm  long).  However,  initial 
calibration  attempts  showed  there  were  still  some  problems. 

Up  to  this  point,  the  probe  aspiration  velocity  was  controlled 
by  choking  a  capillary  tube  of  length  chosen  to  obtain  the 
desired  velocity.  This  was  a  temporary  configuration  used 
only  for  initial  simplicity;  the  final  setup  was  expected  to 
use  choked  flow  through  an  adjustable  needle  valve.  However, 
steady  state  calibration  showed  a  "step"  effect  in  the  output 
signal  when  the  probe  was  exposed  to  a  step-change  i n  helium 
concentration: 
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The  first  voltage  change  from  EQ  to  E1  was  expected, 
and  is  caused  by  probe  response  to  the  helium-air  mixture 
as  it  is  drawn  over  the  hot-film.  The  unexpected  second 
voltage  change  from  E-|  to  E^  should  have  been  foreseen,  but 
wasn't.  It  occurs  when  the  helium-air  mixture  reaches  the 
choked  capillary  tube  located  about  1  m  away  from  the  probe 
itself.  The  speed  of  sound  in  the  helium -air  mixture  is 
higher  than  in  air,  and  the  mixture  passes  through  the  choked 
capillary  tube  more  quickly.  This  velocity  change  is  pro¬ 
pagated  back  upstream  to  the  probe  and  enclosed  hot-film, 
where  it  is  detected  as  an  abrupt  increase  in  probe  output 
voltage.  If  the  capillary  is  thought  of  as  part  of  the  probe, 
then  probe  response  to  the  step  change  in  concentration  is 
E^  -  Eq-  However,  the  transient  response  of  the  probe  is 
then  governed  by  the  time  required  for  the  helium-containing 
eddies  to  travel  from  the  probe  to  the  capillary  tube. 

For  the  prototype  configuration,  the  lag  time  of 
several  seconds  was  much  too  long  for  the  probe  to  accurately 
follow  high  frequency  concentration  fluctuations.  One  solution 
is  to  place  the  sensor  and  capillary  tube  so  close  together 
that  the  two  voltage  effects  merge  into  one,  and  probe 
response  is  assumed  to  be  a  single  step  change  of  -  Eq. 

The  alternative  is  to  place  the  capillary  tube  so  far  from  the 
probe  that  a  complete  set  of  measurements  can  be  made  before 
the  c a p i 1 1  a ry  tube  sees  any  fluctuating  concentrations.  A 
large  volume  container  placed  between  the  probe  and  the 
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capillary  tube  would  act  as  a  concentration  fluctuation  damper, 
so  that  the  capillary  would  at  worst  see  only  a  slowly  chang¬ 
ing  helium  concentration.  The  correspondingly  slow  velocity 
change  could  then  be  ignored  or  if  necessary  removed  from 
the  probe  signal  by  high-pass  filtering. 

The  method  of  separating  the  choking  device  from  the 
probe  itself  was  rejected  because  it  limited  the  length  of 
time  the  probe  could  be  used  without  flushing  the  system. 

The  only  choice  left  was  to  incorporate  the  flow  choker  in  the 
probe  itself,  which  meant  that  a  capillary  tube  could  not  be 
used  because  of  its  undesireable  length.  The  adjustable  needle 
valve  was  also  rejected  because  of  size  restrictions;  this 
eliminated  the  possibility  of  isokinetic  sampling.  The  final 
choice  for  a  velocity  limiting  device  was  a  choked  orifice 
formed  by  drilling  a  fine  hole  through  a  thin  brass  plug 
cemented  into  the  probe  body  as  close  as  possible  to  the  hot- 
film  sensor.  This  procedure  re-awoke  the  old  spectre  of  in¬ 
ternally  generated  turbulence  noise-  basically  a  microscopic 
"roaring"  sound  as  the  flow  passed  over  the  sensor  at  about 
14  m/s  and  turned  a  90°  corner  while  rapidly  accelerating  to 
the  speed  of  sound,  all  in  a  length  of  about  5  mm.  This  noise 
was  minimized  by  using  a  tapered  slot  brass  sensor  mount  which 
gently  turned  the  flow  90°  toward  the  sonic  orifice.  Care  was 
taken  to  eliminate  sharp  edges  or  abrupt  cross-sectional  area 
changes  inside  the  probe.  The  orifice  diameter  was  varied  by 
trial  and  error  to  give  an  ingestion  velocity  of  about  7  m/s. 
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which  was  close  to  the  wind  tunnel  velocity  used  in  the 
experiments.  The  final  probe  design  is  shown  as  Figure  2.1 
of  Chapter  II,  which  deals  with  the  methods  used  to  calibrate 
the  sensor  for  steady  and  fluctuating  concentration  levels. 
The  two  figures  at  the  end  of  this  Appendix  describe  electric 
circuits  used  with  the  fast  response  concentration  sensor, 
and  are  discussed  in  Chapter  II. 
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FIG.  B.l  Probe  Resistance  Measurements 
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a)  Offset  Voltage  Supply 
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b)  Voltage  Offset  Box 
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FIG.  B.2  Voltage  Offset  Circuits 


